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RELAXATION METHODS APPLIED TO DETERMINE 
THE MOTION, IN TWO DIMENSIONS, OF A VISCOUS 
FLUID PAST A FIXED CYLINDER 


By D. N. pe G. ALLEN (Imperial College) and 
». V. SOUTHWELL (Cambridge) 


eceived 13 May 1954 


SUMMARY 


In this paper relaxation techniques are applied to the general case of steady 
Jaminar mot of compressible viscous fluid past a stationary cylinder: that is, 
to motion at sper su that neither inertia nor viscosity can be neglected. The 
governing ju O I Ly voley vV2]C 0. 
where V? cy’, and v are the component velocities, v is the kinematic 
viscosity, and ¢ (the v city cv cx—cu/cy. It must be solved in conjunction 
with the 

uCca ovlcy 0, 
which permits the introduction of a stream-function us such that 
u cy " v cus /Cx, C Vu. 

The numerical computations relate to a circular cylinder, but the methods are 
applicable to any shape (an initial conformal transformation changes the independent 
variables from « and y to « and f, the irrotational velocity-potential and stream- 
function for flow past the specified cylinder). The flow-patterns (contours of y and 
{) change as the ‘Reynolds number’ # increases; but an introduction of variables 
involving R makes the change relatively slow, and thereby (e.g.) the accepted 
solution for R 10 is made a good starting assumption for R 100. 

Fig. 6 r 
theoretical estimates. 


lates computed values of the total drag with experimental and other 


1, Introduction. This paper treats two-dimensional (laminar) motion of 
an incompressible viscous fluid in the most general case—viz. when neither 


inertia nor viscosity can be neglected. The governing equations are the 


equation of continuity ow 
== (1) 


ee |p -(X, ¥, Z)+rV%(u, v, w), 
plox oy oz 


[Quart. Journ. Mech. and Applied Math., Vol. VIII, Pt. 2 (1955)] 
5092.30 K 
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in which 

u, v, w are the components of velocity, 

X,Y, Z are the components of body-force, 

p is the ‘mean normal pressure’, 

p is the density and 

v is the ‘kinematic viscosity’ of the fluid (both taken 


here as constant), and 


V2? stands for ¢?/éx?+-6?/éy?+- ©? ez". J 
In two-dimensional motion w Z Oand u,v, X, Y, p are indepen: 
of z, so (1) reduces to 
cu cl 
= 0 
Ca cy 
and (2) to 
C C C l| ¢ C — " 
——+ U + U- (u,v) 5 P- (X, Y)+vV?(u, v), 
ct Cx Cy ploxr cy 
VY? now standing for ©? ca? + 6? /ey?. 


We postulate that the body-forces are conservative so that 


eX oY aa @ 


CY ox 


and on that understanding, having regard to (4), we eliminate p from (i) 


obtain for steady motion (independent of f) 
| w é Cxrtve ey \C vV2C, 
where V2 has the significance stated in (6) and 
¢ (the vorticity) = cv/ex—eu/cy. 

2. In relation to plane two-dimensional motion (4), (8), and (9) | 
to be satisfied in conjunction with appropriate boundary conditions. 1 
paper treats in detail flow past a rigid cylinder (with axis parallel to 
of a stream which otherwise (i.e.. in the absence of the cylinder) w 


have uniform velocity U’. When, as here, the cylinder is kept stationa 


the conditions to be satisfied at its surface are 
u v 0, 
Other conditions (at infinity) will receive attention later. 
3. Even when thus reduced to (4), (8), and (9) the general equat 
have hitherto proved intractable, except by approximate methods an 


relating to fairly small velocities. High-speed solutions have been attail 


only by the introduction—by L. Prandtl in 1904—of assumptions bes 


on the notion of a ‘boundary layer’ and of a ‘wake’ within which ' 


effects of viscosity are confined: and even so it is necessary to posti 
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distribution of the pressure on the cylinder, which strictly should 


merge as a result of computation. The only treatment known to us which 


spenses with this postulate (1) relates to rather slow flow past a circular 


vlinder.t It has some points of similarity with our relaxational treat- 


nt (e.g. its evaluation of the stream-function % at discrete nodal points 


square-mesh net). 


4, Here, too, flow past a circular cylinder is discussed, for the reasons 


that both in theory and in experiment this shape has received more 


ttention than others and (ii) that it is sufficiently ‘bluff’ to exemplify 


ficulties which such shapes may oppose to computation. Actually an 
tial step in the relaxational treatment makes the shape a matter of 
tively small concern: namely, a conformal transformation which 
nges the independent variables from 2 and y to « and 8, the velocity- 
tential and stream-function for irrotational steady flow past the given 
nder. Thereby the field of computation is transformed into an infinite 
ine containing a rectilinear slit, so computation can be effected on 
iare-mesh nets having no ‘irregular stars’: the shape affects them only 
that the transformed equations involve h, the modulus of transformation. 
With a view to the numerical computations, all quantities are made 
on-dimensional’; and in consequence U (the velocity at infinity) appears 


R= UL) (11) 


umerical parameter 


Reynolds number of the motion) conjoined with v and L, a representa- 


dimension of the inserted cylinder. A further transformation which 


replaces 8 by 3R gives to the equations forms such that a solution 


ind for some particular speed U, can be made a starting assumption 


which, with relatively little computation, the solution for another 
eed [’, can be derived. In sections 14-16 yet another transformation 
vised by Allen) is employed to derive equations which, as involv- 
gexponentials, are not closely represented by the customary approxima- 


ns in finite differences 


5. Results are presented in Figs. 1-7, of which their legends provide 
‘wheent explanation. Figs. 1-5 exhibit, by contours of the ‘non-dimen- 
nal stream-function (%) and vorticity (C), the general characteristics 
the flow when R 0, 1, 10, 107, 108. Fig. 6 compares computed values 
total ‘drag’ with experimental and with other theoretical estimates. 


‘ig. 7 shows how the computed pressure-distribution alters with Reynolds 


v flow past a circular cylinder, has recently been 
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Fic. 3 (a). Stream-function (100%”) for R 10. 








Fic. 4 (a). Stream-function (100¢%”) for R 100. 














Fic. 5(a). Stream-function (1006”) for R 1000. 
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Fic. 4 (1 Vorticity (62-5f”) for R 100. 








Fic. 5 Vorticity (62-50”) for R 1000. 
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6, In this connexion a word should be said about stability. It is the 
sence of the relaxational technique that it starts from solutions which, 
-al) ng inexact, would require body-forces for their maintenance, and there- 
systematically corrects the flow-pattern until the body-forces have 
liquidated’ (rendered negligible). Such treatment, though it does 
mtemplate instability, will normally detect any tendency of a solu- 
e, and in this way will indicate the occurrence of instability, 
nnot precise determine the point of transition. It seems 
th while to record that in this problem the flow appeared to be stable 
R 0, but a distinct impression of instability was gained in the 
putation for R 100; for (Goldstein et al. 3, p. 419) ‘the value of R 

t which the unsteady régime commences . . . is probably about 50°.7 
he instability entailed no computational difficulty: that is to say, 
uiy regimes could be computed for R 100 and for R 1000 which 
experiment would not be realizable because the smallest disturbance 


ad upset them. 





at 
7. The governing equations. For two-dimensional flow the governing 
tions are (4)—(9), and the conditions (10) have to be satisfied at the 
: 
dary of the inserted cylinder. The relation (4) permits expression 
und n terms of a stream-function % by 
CY, v Cus Cx, (12) 
i then, according to (¥ C Ved (13) 
ere V* has the significance stated in (6). At the boundary of the inserted 
l¢ t satisfy conditions derived from (10)—viz. 
0, ous /Cy 0), (14) 
ondition t infinity’ we shall assume that 
Uy, sothat (+0. (15) 


8. Reduction of the equations to ‘non-dimensional’ form. Using L 


lenote some representative dimension (e.g. the diameter) of the cylinder, 
U tod te the velocity at infinity, we now write 
] j Ly u Uu', v Ur', ob U Lif’. (16) 
] 


— ( where ‘a V"*b', (17) 








f instability in his calculations for R 40). 
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and according to (8) 


Rl uw’ oo) = V’°C (18 
Ca CY | 
where FR (the Reynolds number of the motion) ULI/v (11)} 
and V"? = 0? /0x'2+ 6? /dy"2. 
Also, by (12) and (16), 
uw” Cus’ oy’, 2)’ Cus’ Ox’ , 


so the boundary conditions (14) and (15) become 


, 


u’ = 0,0 0, at the surface of the inserted cylinder, ) 


u’ 1, v = 0, far away from the inserted cylinder. |} 


9. First change of independent variables (to avoid ‘irregular 
stars’). We can avoid ‘irregular stars’ in the relaxational computation 
by changing the variables from 2’ and y’ to «and £, the velocity-potenti 
and stream-function for irrotational flow past the cylinder under discus 
sion.t Far up-stream 

x 2x’+constant, B>y’, (2 
and on the surface of the cylinder (and in the plane of symmetry) 
B 0. (2 


] 


Here « and £ are conjugate plane-potential functions such that (a 
is a function of (2’+7y’). Accordingly 


Cx op Cx cp 9 
cx’ oY cy’ ear’ 
- . Oa \* Oa \? op \2 op \* ' 
and writing h2 | l | ] il ERE ola 94 
Cx cy Oe cy 
C Cx Cc cB c Ca Cc op ¢ 
we have — : = al = a 
CX CXL CX Cx Cp cy CY Cx CY ¢ B 
therefore Vv’? h?| é? Ca? + 67/68? ] h?V25 (say), | 
re — >| ow’ ¢ os’ ¢ . 
and u - =f. 9) h? sa ' 
C2 cy OB Ox Cx op 
Consequently (17) transforms to 
(= Vay, 


and (18) transforms to 


R Cus’ ¢ i ous’ {" v2 
, _. 73 
cB Cx Cx 4 2 


BS: ? 
This device was employed by Thom (1). Kawaguti (2) employs an iterative prot 
generally similar to Thom’s, but a different transformation (into a finite rectangle). 
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The transformed boundary conditions are 
0 on that part of the a-axis (8 0) which 
corresponds with the inserted cylinder, 
0 (by symmetry) on all other parts of the a-axis, 
R2 


0, > C'—+>Q as (a*+ 8?) > oo. 


10. Second change of variables (to facilitate approximate treat- 
ment). The product terms in (26) and (27) imply that the flow-pattern 


ters with the Reynolds number. But it is evident, even when we envisage 


reak-away’, that the a-gradients of %’ and @’ will not be large, whereas 
heir 8-gradients may attain high values; and on that account a further 
ransformation is advantageous. 
If e 4 ire replaced by B : wn’. ‘ae where 
3 R af’ R-*p” C Ril (29) 
26) is transformed to 
C hav" dbs (30) 
ym pt nee 
) Cl Cw CQ OY ’ 
1 (27 V"2l - — (31) 
CX cB cp Cx 
_— l oe Cc? 
ert \ Vig = 5 (32) 
R R da? © op"* 
Also from (27) and (28) it follows that 
0 on that part of the a-axis (B” = 0) which 


corresponds with the inserted cylinder, 
. ; (33) 
0 on all other parts of the a-axis, 


- 0 _-> C" +0, as (Ra?+-B"") > ow. 


Thus R now disappears from the equations, excepting as it enters into 


ind on the assumption that the derivatives of %” and (” with respect 
und te are comparable, when #& is large a close approximation to 
i V"2 = 02/0B"?, (34) 


the distribution both of %&” and of ©" independent of R, and 
hereby greatly reduces the labour of a relaxational treatment. This aims 


systematic elimination 


f residuals which express the errors of a trial 


SOLUTION na 


now the errors entailed when the solution for R 10 (say) 
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is taken as a trial solution for R 100 are small enough to require by 
little adjustment. Then the solution for R = 100 may be utilized, similarly 
as a trial solution for R 1000; and so on. (R may be either increased 
or decreased in successive solutions. Our sequence was R L000, 10 
0, 1, 10.) 


11. Introduction of the relaxation technique. We now describe th, 
application to our present problem of the relaxational techniques whicl 
have been explained in earlier papers. They are not required for th 
transformation of section 9 in relation to a circular cylinder, since thy 


known solution can be utilized. It is 





when Ox, Oy pass through the axis and when L—the representatiy 
dimension—is the diameter of the cylinder. With a use of these expressions 
the transformation is easy (Fig. 8). 


Contours of a 











12. Hereafter we shall give to B, us, ¢ the meanings which in section 10 wer 
attached to B”, us", (”. On that understanding (30) will now be written as 
CAs L 07s C 9 


R 0a? h>?’ 











and 


It is 


ma 
in 


th 
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quire by nd (31 iter some rearrangement, as 
similar Vaid Cus\ OL 1 (e7C Cus CL) = 
hf pot le Mae ot) on @ (37) 
Nncreas cB" Ca/op R\éa? CB eal 
O00, ] It is not necessary to rewrite the boundary conditions (33). 
2 
fo 
Oo - 
4 
Fic. 9. 

‘or a relaxational treatment, (36) and (37) must be replaced by approxi- 
mations in finite differences. Since (36) is linear, it presents no difficulty: 
in the notation of ‘residuals’ (4) its approximation is 

nt i | a ee 2 — ” 
Ly )o Yo 1 =o T 5 (ih, bs - - 2g) arg h*), 0, (38) 

eSS a 
he s es 0 2, 3, and 4 denoting points so numbered in Fig. 9; 








and from (38) it is easy to deduce the standard ‘pattern’ shown in Fig. 10, 


for the effect on the F,,’s of an increment Ay = 1 at 0, Fig. 9. 
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13.+ Equation (37) presents a harder problem, and being non-linea 
yields a ‘pattern’ which alters as the work proceeds. But it has bee, 
shown in earlier papers that inaccurate patterns can be used to liquidate 
residuals, provided that an exact account of these is kept; and here q 
simple pattern represents sufficiently a relation of rather complicated 
form. 

We contemplate a ‘two-diagram technique’ in which % and ¢ are modi- 
fied alternately. At the end of a stage of ¢-relaxation the F,,’s will have 
altered in accordance with the last term of (38), and the &’s must be 
altered so as (temporarily) to liquidate them: then, with the altered values 
given to uw, the ¢’s must be modified in accordance with a finite-differeng 
approximation to (37). This we proceed to derive. 


14. The two quantities which together make up (37)—namely, 


Cl Ob Of 


(a) TS 
op? Cx ep 
P (39 
70 Cus el 
and (b) co ee 
Cx" op Ca J 


can. be treated similarly in a way we now explain in relation to (a 
Writing ' a 
Cus OG 

’ A K—, 

Ox fas} op 

we have as the solution of the second of (40) when x and A are invariant 


KC AB | P+ Qe —KB 


a9 


C 


um 


(40) 


tw 


P and Q being constants of integration. Then, 2, 0 and 4 denoting 
adjacent nodes on a f-line of the a—8 net (so that 


B, = Bo+a, B, = By—a 


/ 


when a, as is usual, denotes the mesh-length), we have 


K(Co fo) Aa- Ue KPo(¢ Ka__] ) 
K(C, Co) Aa- Ue 3 KBo(¢ KG 1), 
therefore KeX@(Co— Co) + f4— So} Aa(e**— 1). 


We thus have an expression for A—i.e. for (a) of (39); and a like ex 
pression may be formed in the same way for (b). The values of « (= a @ 


and of A ( -Rés/e8), the corresponding quantity in the expression fol 





(b), will of course not be known until the solution has been completed; but 


they may be treated as known when y has been temporarily determined 


+ It should be recorded that the technique described in Sections 13-16 is entirel) 
to Allen (R.V.S.). 
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nd values appropriate to the point 0 may be computed from finite- 


lifference expressions of normal form: viz. from 


2Ka us Ya, 2Aa Rb, y4). 


We then have as an expression for the typical ¢-residual, defined as the 
‘nite-difference approximation to a? x [left-hand side of (37)]: 


Py)y = Ky aer"(la— Lo) +.04—Sos/(e** — 1) 4 


Ag ajpen"(C, Co) T C3 Cos Re — 1). (41) 


It is the latent exponentials in the solution, revealed in (41), which 
ecessitate this special treatment. Except within a very narrow range, 
exponential is not clos ly represented by a polynomial; and on that 
count the customary expression in finite differences for the second 


term in (39), (a) or (6), has insufficient accuracy when x, and/or A, is not 


15. Regard must be paid, in liquidation, to the boundary conditions 


3) section 10. Rewritten in accordance with section 12, and with 





2 
a L 1 
~ + « 
| 
| 
| 
i 
4 
« 
Fie. 1] 
tives replaced by their approximations in finite differences, these 
ecome 
des on the a-axis except the ‘slit’ (i.e. that part which 
corresponds with the inserted cylinder): 
us C ): 
it nodes far away from the ‘slit’ (where «?-+-B? — 0): (42) 
3, C> 0; 


nodes on the ‘slit’ (typified by 0, Fig. 11, in which the node 4 


Is hetitious 


0 U, ibs yb. 
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Only the last of (42) needs special notice. Combined with (38) 
section 12 (which here reduces to 


9 
a’, 
us, us, | = Co 0 
h? 
because 4, = Wy = ys, = 0), it becomes 
: OL? =e 
Go Zhe bo a, (4; 


so permits, at the start of each stage of ‘C-relaxation’, a specification of 
at all nodes on the a-axis. 


16. Except in its use of ‘patterns’ which alter as the work proceed 


and which call for a use of exponential tables, the relaxational procedw 


follows normal lines; ‘residuals’ of the types Fy and F;, being liquidat 


alternately, in stages. Small changes made in ¢-values were found to alt 


y-values largely, and on that account the F;’s were liquidated more 


less completely ; but the stages of ‘ys-relaxation’ were made short becaus 


the F's became smaller almost automatically by reason of the change 


made in the ¢’s. 


17. When acceptable distributions have been found for % and { 


functions of « and 8—that is, in the notation of section 10, for y%” and{ 


as functions of « and £”—the transformations (29) will yield ys’ and ¢' 


terms of « and 8, which have the expressions (35) in terms of 2’ and 


Contours of %” and ¢” can then be plotted on the rectangular (2’, y’) ne 


Figs. 1-5 were thus derived. 


For Figs. 6 and 7, expressions for the boundary tractions were requir 


It can be shown (cf., e.g., (5), sections 325-6) that the normal pressw 


on the circular cylinder is 


Pn (Say) = p—2vpcos20-—¥ + ypsin 2 aii | 


Cxucy Oy? dx? 
and the tangential traction (upstream) on the cylinder is - (4 
9 “9 9 


o? ¢ ‘ c 
Png (SAY) = vp cos 20| — = |-+ 2vp sin 20 — 
CXLCY 


cy" = dox- 
at a point whose angular distance from the upstream stagnation-point is 
/ ‘ 9 9 
‘ C us Cus a a 2h 
At points on the boundary { where yb = = —. = 0, s0 f = h? — 
‘ Cx ox" op 


ous axes C dnda . { (Ex\? ‘Cx\*) Cys 
Oxoy Gu? OB hax dy | 7a) - es Je rep’ 

Tt atin’ dx "i, g Ora Ob 4 0% Co Cp 

dy? da® = ~\ hoa} J Oxdy OB dx Cy CxeB 
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1 


| p, the ‘mean normal pressure’ (section 1), is related with p,, the 


tatic pressure at infinity’, by 


»f,.\ 2 % ¥ 
{Cw ae 2 
p.—) sph? ) + | = dx. (45) 
Cp J OCP 
ict 
p,—pU? | >_ da 
: c B > 
/ | (46) 
pU2t"/ Ri 
1 point on the boundary (where, by (35), « 22’ cos@, 8 = 0). 
m (45) and (46) the ‘drag’ (Fig. 6) and pressure-distribution (Fig. 7) 
computed: Fig. 6 by. integration of 
Png sin 6+-p,,,, cos A). 
18. This investigation, started in 1944, was frequently interrupted on 


ount of pressure of more immediately urgent problems. An interim 
count of it was presented at the International Congress of Applied 
Mechanies (London) in 1948 
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ROTATORY AND LONGITUDINAL OSCILLATIONS OF 
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SUMMARY 


The Stokes stream function is used to obtain the perturbations arising from 
slow rotatory and longitudinal oscillations of axi-symmetrical bodies in an infir 
mass of viscous fluid which is at rest at infinity. The bodies considered in this pay 
are a sphere, an infinite circular cylinder, a prolate spheroid, an oblate spheroid, 
a circular disk. The case of the cireular disk is deduced as a limiting case of | 
oblate spheroid. The paper has been divided into two parts. Part I deals with 


rotatory oscillations and Part II with the longitudinal ones. 


1. Introduction 
THE perturbations arising from the longitudinal vibrations of a sph: 
along a diameter in an infinite mass of viscous fluid at rest have be 


considered by the author (1) by applying the idea of periodic singular point 


and the results are found to agree with those found by Lamb (2). T 
rotatory oscillations of a sphere about a diameter, in an infinite mass 


viscous fluid at rest, have been discussed by Lamb (2). The object of thi 


paper is to utilize the Stokes stream function to obtain the motion oft 
fluid due to the rotatory and longitudinal oscillations of certain axi-sy 
metric bodies in a viscous fluid. The slow motion, in which the inert 


terms in the equations of motion are neglected, is considered in the cas 
when the bodies are in the form of a sphere, an infinite circular cylinde 


a prolate spheroid, an oblate spheroid, and a circular disk. The rotato 
oscillations are about a diameter in the case of a sphere, about its axis 
the case of an infinite circular cylinder, about their axes of symmetry 
the cases of spheroids, and about the normal axis through its centre int 


case of a circular disk. The case of the circular disk has been discussed 


a limiting case of the oblate spheroid. The longitudinal oscillations 

along the same axes about which the rotatory oscillations of respecti 
bodies have been considered. Some of the results obtained agree wi 
results already known, while others appear to be new. For spheroids it 
found that we have to use spheroidal wave functions discussed by Stratt 
Morse, Chu, and Hutner (3). The notation given by them will be us 
throughout this paper. 


[Quart. Journ. Mech. and Applied Math., Vol. VIII, Pt. 2 (1955)] 
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, Equations of motion 


be the general orthogonal coordinates and let the elements of 


> 


sth at the point (a, 8, y) in the directions of a, f, y increasing, respectively, 











oA and € ) . such that 


ds2 e?(da)* | e3(dB)? | es(dy)*. 
be respectively the components of the velocity in the directions 
creasing. For motion symmetrical about an axis (4) we take a 
the general orthogonal coordinates in a meridian plane and y 


| angle ¢, so that e, will be the distance from the axis of 
\ll quantities are supposed to be independent of ¢. The 


ntinuity then takes the form 


€,uU)- (€3 €, V) Q. (1) 


l ous l oud 
Co U C_V - (2) 
Q? 3 
9 UP ec X 
vhether the velocity w round the axis is zero or not so long as 


ent of d. In the case when there is an azimuthal velocity w, 
C4 UW c2* 


3 


he components of vorticity, are given by 





lL 0&2 1 oQ 
) 1) . Par 
5 Cp ; €g€, Cx 
1 | r cus ri ( cus l 
( : 3} D>, 
€1©o €4 oa CP \E3 &2 CP €3 
: Cn | ¢ Co c é, ¢ | 
D2 i— } }}- 
€1 gL CX\Eg €y CX CB \es € OB} | 
ms of motion give 
2X2 = A(Q2, €. l o(us, D7) 2D O(ub, €g) Dd 
2 ‘ ~~ 2 I 2 Q vip, 
ss ASP €,€2€3 C(a,p) €1 €2 €3 C(X, P) 
neglect the second-order terms. becomes 
Dus) vDAd, (3) 
equation for likewise, the equation for Q is 
ok2 l e(ys, Q) 


vD?Q. (4) 
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(a) Rotatory oscillations. In this case 


u = 0, v 0, 3 W = ER 


Thus “4 = 0 and the equation (3) is automatically satisfied. The moti 


is given by the equation (4), which reduces to 
eQ 
ot 


vD?Q. 


(b) Longitudinal oscillations. In this case 


1 eu lL Oy 
34 = — —; é,0 => —— A w Q. 
' AD ‘ . 

€g Cp €, Cx 


Hence Q 0, the equation (4) is satisfied, and the motion is given byt 


equation (3), which now may be written as 


p*(D* i 0. 
v ot 


If 4, and %, be two functions satisfying the equations 


Db, = 0, 


a solution of the equation (6) may be written as 


ob wy | bo. 


; Q, 


(0 


v ot 


Thus we have to solve the equations (7) and (8). The superposition of th 


solutions will give the required stream function. 
To satisfy the conditions at the surface of the body in this case we t 
the origin (2) at the mean position of the centre of the body. 
Part [: Roratory OscrLLaTiIons 
lL. Sphere oscillating about a diameter 
3. Taking spherical polar coordinates R, 6, ¢ for «, 8, y respectiv 
we get . 
a é¢, = 1, e, = R, e,; = Rsiné, 
and “= 0, vy = 0, w = Q/Rsiné. 
The equation of motion for harmonic oscillations of period 2z/c is 


2Q, sind a/1 @aQ,\,,. 
+h?7Q, = 0, 
(a 6) om 


eR? ' R® 60 
where 0 = 0,e%, h? = —iolv. 


The solution satisfying the boundary condition at infinity is 


: ] 
Q, = Ae-*R(1 4. —_|\sin29, 
, ( ih 2) ” 
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nd the angular velocity w is given as 


1¢ hk ] : 
a : 3 iot 12 
pe (!+ aap) Oe) 


\pplying the boundary condition at R a, viz. 
WwW Wo € tot (13) 
tha, ¢ thag3 


nn A : - 
(1+tha) 


(14) 


4, The couple required to maintain the motion 
The tangential stress on the surface R = constant and in the direction 


i d Increasing is 
9 9; 2.9 
= uw,(3-+-3tha—h?a?) . - 
1 -. 5 : sin 6 e?%, (15) 
(1+-tha) 
The couple G is thus given by 

; | — 3tha—h*a*) ted ; 

G TT [LAW om Co, (16) 


| 1-+-tha | 
\ 


ich agrees with the result obtained otherwise by Lamb (2). 


Il. Circular cylinde r oscillating about its axis 


on 


indrical polar coordinates we take « as z and 8 asr. Then the 


yuation of motion for harmonic oscillations reduces to 


17Q, 1 dQ aie 
> —tk*Q, = 0, (17) 
ur@ r dr 
ere h:? a/v and Q ie. 


ippropriate solution is given in terms of modified Bessel function as 


Q), Ark,(vikr), (18) 
d the angular velocity w is given by 
AK, (wikr) . 
1 = ) piot (19) 
, 
[he boundary condition at a, VIZ. w w, ert, gives 
——< (20) 
K,(vika) 
6. The couple required to maintain the motion 
The tangential stress 7’ on the surface r = a and in the direction of d 
reasing 1s given by 
Tr “Hedy Vi kak,(ve ka) ict (21) 


K, (vt ka) 
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So the couple per unit length of the cylinder is 


9 Vi kaK,(vi ka) 


iot ) 


G = —27pw,a* —.———_ _ e™. 9 


K,(vi ka) 
Ill. Prolate spheroid oscillating about its axis of revolution 
7. In this case we use the prolate spheroidal coordinates defined by 


z+ir = ccosh(«+7f), 


€, = €, = C/(cosh*a — cos’), e, = csinhasinf. 
: , ..% eéQ ‘ 
Ihe equation of motion is ; vDPQ, 
C 
where 
— (er a é é \) 
dD? — o "Pee SS eer” a ®@ a as y— > — coth xX . cot B ‘wai . 
\c?(cosh2a — cos?8)\éa? © eB? Ox op} } 


Writing Q, e’ = Q in the equation of motion, we get 


io , 
DQ, = - Q, = ik?Q,, ya 
y 
where k? = a/v, or 
oo 


| at “ _ cotha- = —00t B= ik®¢?(cosh? x — costs) |0, =. (@ 


,* OF é 
Putting : P 
Q, = sinh asin B S(B)R(a), 
cosha = &, cosB = », 


equation (24) separates into the two following equations: 





d | U8). {y_ ipze22 _} L_\s 
—] + {A—tk2c27n?2 — — = 0, 2: 
na” ar ir 
d(x dR\ | . 1 ) 
€*—1) —}+ lA— ee Fe : 2 
de® ) dé} | tk*c2é2 — B1| 0 


where A is the constant of separation. Equations (25) and (26) are identica 
with the equations obtained by Stratton and others (3) in the discussia 


of the prolate spheroidal wave functions. Such equations contain a thir | 


parameter m and have solutions for the suitably related values of thi 


parameter in series of associated Legendre functions or Bessel function 


of half order. In our case m is unity. Thus, for various characteristi 
values of A, there are ‘angular solutions’ S},(vi ke, 7) of (25) which are finit 
throughout the range —1 < » < 1 and are expressed as infinite series 
associated Legendre functions in the form 


S1(vi ke, 7) >” di Pi.,(n), (21 


n=0,1 
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shere the prime indicates summation over even or odd values of n according 


sn is even or odd, and 


a {H(m+1)}! — (4-+1)}! | 
S ind a = “apt? l even, 
1=0 - 

Me L(n+2)'! {4(1+2)}! 

2, tn ta—iyi ~ Geen bat 
& 2 m8 ‘ 


Since Q is to be expressed ultimately as the sum of series of products 


§(8) R(x), the boundary condition 2 — 0 as € + oo determines at once the 


ppropriate ‘radial solution’ of (26). The function which satisfies this 


ndition is R},(vike, €), which is given, for large values of vi kcé, as 





(E2—1)!/(a/2vike&) SS’ i-"dh(n+1)(n+2)K,,,3(vi kc€) 
R38 { V7 / _ n =0,1 . 
>’ d)(n+1)(n+ 2) 
7% =0,1 

(28) 

vm 7 

> é +e 2Y 
Ni kc€ ( ) 


‘he complete solution of equation (24) is 


Q, csinhasinB > A, R})(a)S},(B). 
l1=0 


There fore Ww Ss A R3(x)S! (Bye tot 
hes cae 1l 
t 0 


] 


l the angular velocity is 


x 


l “ 
3 ‘1 Q t 
w = ———_, 5 A, Ria) Sh (B)e. 
csinh asin B 4y 
° [=9 
[he boundary condition requires 
w, ert at a Xge 
- x 
ineretore Wy ( sinh ¥, sin B > A, R?( x9) S1(B). 
l=0 
Using the orthogonal property of S1,(8), we have 


w,csinh a, | sin?8 S},(B) dB Aq, Rix), 





(33) 
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7 1 
To integrate | sin?8 S},(8) dB = | .(1—7?) S3,(n) dn, 
0 71 
where n = cosB, 
we notice that St(7) >" dh P! .,(»), 
n=0, 


and assume the validity of the interchange of the operations of integrati 


and summation. We have 


+1 1 
. : + for n (). 
\ (1 oP. i()) dy | Pi(») PF 1()) dy lo i a 
‘Y . : 
So the equation (34) gives 
4m, c sinh a ' 
A, —3——_ for / even 
od) Rey x9) 
A, = 0 for / odd 
8. The couple required to maintain the motion 
The tangential stress 7’ on the surface a x, and in the direction of 


increasing is given by 





ar eg wW | 
r= ule] (2) 
Mle X\€3 be 
. A, R3(«)S1(B 
bh sinh XQ | > I—0 vit to) a) iot 
- € 
c(cosh?a, — cos)? | ¢ = sinh a | =a 
aie osh op sin B fiw, Zz { R'(a) Oh R? ¥o)}S4,(8) | ms 
| (cosh2a, cos")! 3(cosh?a,) — cos*f)* | 


The required couple G is given as 
G ( all 27, Cc? sinh*a, cosh a, | sin’8 dB 
0 
S77, C8 sinh? ax, > { RF (ax)/% R? (x)} sin?p S1,(B) dp 
l=0 F 
cosh ay+ 4 ¥ {R3i(a9)/q, R3)(xp)} [Sw C3 sinh?a, e', 


oa 
0 


where a dash denotes differentiation with respect to a. 
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Coa Cc Sil 
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7 | 
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lation +5) separates 1 
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) 


j 


oscillating about its axis of revolution 


introduce a system of coordinates defined by 





i? esinh(a+-<f), 
4d 4 
. + cos?6}, €, = ccoshasinf. 
o@2 
yDPFQ 
ol 
a © YO a 22\10 
cotp ik?c?(sinh?« COS“p) | 82, V, 
Q 
cp 
(39) 
ie, i a/v. 
cosh .sin B S(B) R( x), 
€, and cosf = 7, 
» the following two equations, 
nae oe 
| A ik?c? 1? = Ss 0. (40) 
l 4 
IR er l 
lA sk®c*e* 1. — R 0, (41) 
| BL] 
separation. Equation (40) can be obtained 
to ixik. The equation (41) can be obtained 
wing sik to ivik and € to —ié. Therefore the 


10) and (41) can be obtained from those of the 

vy changing vik to —ivik and € to 7g. Thus 

en by 

Vikjn) = ¥ fLPLialn), (42) 
0.1 


n 


ire different in value from the coefficients d!, of 


uined in the same way (5). The radial solutions 


nfl (n+-1)(n+-2)K, (vi ke€) 


(43) 


» 


1)(m ) 
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The required solution in this case is 


ioe 


Q = ccoshasinB > A, R}(a)S},(B)e™, (44 
i=0 
w= > A, Ri(a)Sh(Be™, (45 
[0 
I] — - 
a ee A R3 C 1 erat. f 
OS etek asain B pa 1 Ri,(a)S7(B) (4 


The boundary condition gives 
wyccosh a)sinB = >} A, R},(a9)S},(B). (47 
1=0 


Employing the orthogonal property of S1,(8), we get, as in equation (33 


—4w9¢ cosh ap 





for 1 even 








ingles 34, Ry, Xo) re (48 
A,=0 for lodd - 
Ap 2a 19 
where q = 2 S (n I)(n- 2) Lye, 
ran (2n+3) y 


n=0,1 


9. The couple required to maintain the motion 


T,, the tangential stress on the surface ~« = a, and in the direction of ¢ 
increasing, is 








€a| 0 [w 
1 = ne 2 
res Ou bles 
p.cosh Xp ¢ ™ . A, R3.(x) SLA \ tot ( 
= dale, am? “ae Wo oe A, fiy/(a)s w(B) /¢ a (49 


The required couple is given as 
2c? cosh*ay) } A, R3(a») | sin?B S18) dB— 
i=0 . 


0 


G = 





27pc%w, Cosh*ay sinh ay | sin*8 ape 


<< R3i(a 9) 

a9 P ‘ 5 X ‘ ° » 

| nnoye? cosh?a, p - a 0". — Szrpiw, c? cosh2a, sinh no fe (50 
= N Ry, Xo) 


where a dash denotes differentiation with respect to «. 


V. Circular disk oscillating about its axis 


The motion due to rotatory oscillations of a circular disk about the 
normal through its centre can be obtained from the motion of the oblate 
spheroid as a limiting case. If we take the expression of the function Rj,(§), 
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ROTATORY 





the ‘radial solution’ found in the case of the oblate spheroid, near € = 0, the 
required expression is an infinite series of associated Legendre functions. 
The appropriate form of the function (5), apart from a constant factor, may 


be written as Fi(k, €), where, for even values of I, 


Aik ae ) fo). on(WiE)+fee Gt y(t) 
2. [fis —— a” = 
> ( m) Phy (Vi €)— ie, SH(k, vig) (51) 
—\ J 
a 4)! =|? 
and € Jj ach i Ag): 
mol EY 


while for odd values of J, 


F(k, é) = e'l2] Y fh, Qba(vig) fbr QM vig) 


? l 
~ f! 2n-—1 >1 2 ; v1 ; 59 
( , Pop (Vt £) ve; Sik, vi) (52) 
n=1 
nas {4(U—1)}!)? 
d 1 ac3f! l fiz(— J : 
Having proved that A, = 0 for odd values of 1, we shall be concerned with 


the expression in (51). 


The expression for Q is now given as 


Q2 c cosh xsin B z A gy Fy(a)S} (Be, (53) 
l=0 
vhere As dey cosh ay 
3q, Fy (x9) 
In the limiting case when - 0, Ay, takes the value 
{ 409 € _ 
: 34m Fy(0) 
Proceeding as before, we get the couple on the oblate spheroid as 
~ ( Fox) , i , 
G =uw, c® cosh S/= aA o) | Spm, C? cosh*a, sinh Xo |- 
— (Jor Fy X9)) 
The couple in the case of the circular disk is given by the limit of @ as 
0 wi 


ZT LW C8 S F'ai() ett (54) 
4 Ia F,(0) 
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The values of F,,(0) and F4,(0) are given as 


at ome 4) 
F(0) = ( “Ayn (1 + 2ey/7) > (- yn” : 73 


— n! wy 
a, : =. (n+1)!(—4)! ,, 2 | 
F3(0) = ( uy |2 ar Soa Fin SP ot -. 
n=0 : 


Part II. Loneirvpinat OsciLLaTIoNs 
VI. Sphere oscillating along a diameter 
The sphere is oscillating longitudinally along a diameter. The equations 
governing such a motion are 
Dp, | 
| pe. 
| v et| 
Taking the spherical polar coordinates R, 0, ¢ and writing , = ve! 


wb. = wb, e'", the above equations reduce to 


eb, . sin# =f 1 3) 0 


J (56 
ok? =R* e6\siné e6 
245,  sin® ¢ 1 és, er 
and : a4 : ; : (— We th, = 0, (57 
ok? R? c@\sin@ eé ia 
where h? ia/v. 


The boundary conditions are 


tra lw, a? sind e', 
‘Cus _ 
| ) Wa sin26@ et. 
, R R=a 

u=v=0 at infinity, 


where w,e' is the velocity of the sphere along the diameter. 

With the behaviour of the motion at infinity in view, the solution of th 
-quations (56), (57) is — 
equations (56), (57) 1 sin26 
ub, A iy (od 

R 


, eae . Oe , 
- = Be-Sk)) + sin?é. (99 
iby é | AR in 


1 ; ; 
- Be wa +5 z)|sine@ tot. (00 


Therefore us 
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The boundary conditions give 
= , 9, Ua 
{ ; 3iha tha? me 
2 2 2 J pe 
} w, ae 
> 0 
B=5— 
— tit 
e velocity components ure given by 
l C ] ous 
- v —__—— . w = 0, 
R?sin 6 06 Rsiné eR 


hich when calculated agree with results already obtained by different 


VII. Prolate spheroid oscillating along its axis of revolution 
\s in Part I we introduce the prolate spheroidal coordinates «, 8, y such 
v) c cosh(a iB), 
d 4 
€, = C,/(cosh?a — cos?f), é, = csinhasinf. 
[he equations of motion (7)—(8) are transformed into 
| ¢ - s Cc C . 
coth X -Cot B ly, 0, (61) 
; Q 
\¢ x* 2 CX ep} 
| ,€)  c®(cosh*a — cos?8) ¢ 
coth a cot p—, a —— —i”, = 0. 
' cp} v ot\" * 
(62) 
Writing fb, = ws, e%*, ob 5, e'" and k* = a/v, equations (61)—(62) reduce to 
| ( - oc” Cc am = ’ - 
; coth « cot Pp ly’ 0 (63) 
Ca" * Ca ep} 
nd 
O° c “19-6 9 6 , > 
3 cotho cot B a ik?c?(cosh?a cos?8) Ws 0. (64) 
( CD“ ¥ Cp 


sinh «sin B 6(a)A(B), 


sinh a, n = cosB, 


equation (63) separates into the following two equations with A as the 


separation constant: 
d | dd) | l ) 
Oe 3 Fe) ee” ey (65 
dy|\ dn} n*- 1} ” 
d | d6\ . | L | 
(€2—1) $141) aT 0, 66 
dé\ déj"\" B—1{ ia 
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which are the particular cases of the associated Legendre differenti, 
equation. Thus the appropriate solution of (63), keeping in view th 
boundary condition at infinity, is 
#, = sinhasinB ¥ A, Q}..,(«)Pt..n(B). (67 
n=0 


The solution of (64), as already found in Part I, is 


# = sinhasinB > B, R}(«x)S},(B), (68 
l=0 
where S1,(B) >” d&, Pi..(8) 
n=0,1 


and 


or 
ro 


naa Vt ke€) 


Ly (m/2vikeg) S’ dnd) (n-+1)\(n+2)K 
n=0.1 


RR (vi ke, € 


— 


S’ d'(n+1)(n+ 2) 
n=0,1 
The complete solution of the problem is 


£ 


% = sinhasing| ¥ A, Q!.,(a)P!.,,(8)+ > B, R3(x)S1(B)le. (69 
l=0 


n=0 


10. Determination of the constants of integration A, and B, 


The constants of integration A, and B, can be determined from th 


n 
} 


boundary conditions at the surface of the prolate spheroid, viz. 


eee tt. 1S 


U = miga S ate F 
(cosh?a, - cos”)! (ese op), Xo 


* ene + RQ - 
wy Cosh a sin B viot _ _| 1 ep) 





(cosh?a, — cos)! ler €3 Ga) yay 
which require 


(tb). — a0 —}w,c* sinh?x, sin?8 e', 
} 

CW ‘ . eo ¢ , 

| —wy,c*® sinh a» cosh a, sin? e', 
\¢ x X= 


0 
where wye'” is the velocity of the body along the axis of symmetry. 
Now & and éys/éx are given as 


bs > A, sinha Q}. ,(«)Pt.,(8)sinB+ > B,sinh « R3,( x) Sh,(B)sin 8) 
n=0 {=0 
(70) 
Cus os ‘ , ‘ 
> A,{cosh a Q}. ,(«)+sinh « Q}’, ,,(«)}P} .,,(B)sin B- 
Xx n=0 


> Bicosh « R3(«)+sinh « RY (a)}S}(B)sin Ble tot (71) 
i=0 


where a dash denotes differentiation with respect to «. 
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We rewrite these expressions in forms which are suitable for fitting the 
undary conditions at the surface of the body. For this purpose we 
spand the right-hand sides of (70) and (71) into associated Legendre 
(x) and F,,,(«) such that 


(x) = Q}..,,(x)sinh « 
:, ' (72) 
nd > F,(a)P 3) = sinh a R},(a)S1(B) 


Equating the coefficients of the associated Legendre functions of the 


same order on both sides, we get 


F x) sinh X R?(x)d!. (73) 


mm 


The boundary conditions take the form 


: ; . 
sinh?a, P3 ) > [An tn(%) > B, F,y(%)| P} n(B) 
} 0 ] 0 
<iiiacatia aati 2 [An falco) > B, Fin(a9) | P3 n(P) 
n= 1=0 
(74) 
ich yif 
S BF (3wyc* sinh*a, forn = 0 
_ \9 forn > 0 a 
. ' - (fo3 
{ S BI} , {wy c* sinh ¥, cosh Xp for n = 0 
aia : \9 for n > 0 
liminating | we get 
SB F, )] (Xp) F X_) / (a); 
c? sinh a){} sinh ag f,,(a))—cosh apf,(a))} for n = 0, re 
. - | : (76) 
. for n > 0. 


The B, can be determined by solving this system of simultaneous linear 


gebraic equations. Tomotika and Aoi (6) have made some suggestions 


lor the solution of a similar set of equations. In our case it is not possible 


give similar expansions at this stage as the values of the coefficients d!, 


which have been tabulated to date are of little use to us here. 


Che stream function is thus completely determined. 


VIIL. Oblate spheroid oscillating along its axis of revolution 
For an oblate spheroid the system of coordinates is 
i? esinh(a+7f), 
d = Y: 


c(sinh?« cos’8)*, €3 c cosh ysin f. 
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The equations of motion when transformed into these coordinates becom; 


9 49 


of of C C 
5+ =55 — tanh a —— cot B—]y, 0, " 
Ca* cp* Ox ep)’ 
I i 
. P : ° 
o* o é c ee re 
|. 5 - — tanh « cot B — —— (sinh?a -+ cos*p) - \, — (0, (7 
\Ga? © op? Ox op ov a” 


Writing %, = de’, bo = dge'™, and k? = a/v, we get 


e2 C2 C C , 
—, +—, —tanh «—— cot B- 5 eh 0, 
Ca" op? Cx CP) 
c2 2 C c *l0-6 : 9 2 ’ 
! >+—,—tanha- cot B — — ik*c?(sinh*a 4 cos?p) ly - 0. 
\da2 of? Ca op a 


To solve equation (79) we put 
yo, = cosh asin 8 6(«)4(B), 
sinh « é, cos B = 9. 


The equation separates into 


df dd 
=|? 5] +[2- \¢ 0 


dy dy n?—1 

d dé ; -:] 
and &4-] Se eee g=0 
“ ie|@+ [+P +a 





where A is a constant of separation. 
The required solution of (79) is therefore 
wy cosh asin p> A, Q}..,,(isinh «)P} , ,,(cos B). 
0 


The solution of (80), as already obtained in Part I, is 


hs, = coshasinB ¥ B, R3(a)S4(B). 


l=0 


(8( 


(84 


The superposition of these two solutions gives the complete solution 


u& = cosh sin B| 


n 


11. Determination of the constants of integration A, and B, 


The boundary conditions, in this case, are 


_W Cosh a cos B viot 1 aa 
Aa=X 


_ 





-(sinh2a, + cos?p)! en &g OB 


a _ wo sinh % Sin B hil sn si ( ] 4 
ey X= 





(sinh?x, + cos*8)! 


€3 Oo 





An Qh in(a) PH n(B)+ S By Ria) Shi(p)|e. 
7 1=0 
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shich reduce to 


i. Sw, c? cosh*a, sin?B e* ) 

a ee (86) 
¥ 9 ° 4 9 

| w, Cc? cosh a» sinh a, sin®B e'% 

OX) v= cx 


Proceeding as in the case of the prolate spheroid, we introduce functions 
x) and G,,,(«) such that 


},(x) = cosh a Qt, (a) 
(87) 


> G, ,(«) Pt .,(B) cosh a R},(a)St(B) 


Equating the coefficients of the associated Legendre functions of the 


me ord I we vet 


wa ‘ 


G,,,(x) = cosh a R}(a)f,,. (88) 
boundary conditions take the form 
La, ¢ - cosh? Ly Phi ; > |A, 7,4 Xo) T > B, Gin ( 0) Pt J n(B) 
0 {=0 \ (89) 
— ban . : \! 1 
c? cosh a, sinh Xo P3(B) > |A, Jn(%)4 p? By Giin(a9)| 2 t+n(P) 
ol =0 J 
here dashes denote differentiation with respect to «. 
Equations (89) yield 
lop c2anch2 ar . 
he S BG (a (dw, c? cosh?a, for n = 0 
= 10 for n > 0 
. (90) 
- > BG (x) (w,c?cosha,sinha, for n = 0 | 
= 2 10 forn > 0) 
inatir 1, we get a set of algebraic equations 
N BSG G / Xn)? 
,¢? cosh al $ cosh a9 J,(%)—sinh xpJ,,(%)] for n = 0 (91) 
forn > 0 
letermine the values of the constants B, These equations involve the 
eiicients f',, and those so far tabulated only are of no use to us here. 
[x ¢ ular disk oscillating along the normal axis at its centre 
For deducing the stream function for the motion of the disk as a limiting 
se of the oblate spheroid we take the expression of the function R3,(é) for 
eS O1 wr é 0. The required expression, as in Part I, is K(k, €), 


M 
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which, apart from a constant factor, is given for even values of / as 


F(k, £) = etettl Y fh Qh is (Ni) + fly Qt, (Nié)- 
n=0 


=. /f! 
S ( “a (Vi £)—te, SL(k, vig) 
ay p 


while, for odd values of I. 


Fk, é) = i) fin 1 Qan(ViE)+f0, QU vi é)+ 


n= 


a 7 
S (! 2n-1) PL _ (vi £)—ie, SL(k, vig)| 
head p } 


n=1 


The expression for the stream function for an oblate spheroid is now givens 


% = coshasinf}] > A, Q} 
0 


nt 
A,, and B, are determined by the equations 
s ess. olen : 
3 W, C* cosh? for n = 0 
An Gn(%)+ S B, Ayn (x) | 2"o ? : 
1=0 {0 for n > 0 


Astjeh SB da) an {woc*cosha,sinha, forn = 0 
Jn = v ( 


forn > 0 


where the function H,,,(«), being defined in the same way as G,,,(a), is 


,,,(~) = cosh a F(«)d! (M4 


n* 


In the limiting case when a, > 0, A, and B, are given by 


nt 


. lwic? forn = 0) 
A,g,(9)+ ¥ BA, 10) ax {Wo 
, > ** 10 for n > 0 


x 
A, g,(9)4 > B,H;,,(0) = 0 for all n 
l=0 
Elimination of A,, yields us a set of algebraic equations 


S BiH, ,(0)g;(0)—H,,(0)g,,(0)} = [2%0CIn(0) forn=0 
1=0 : , \0 for n > 0 


for the determination of the B,. They are similar to the equations (91 


In conclusion I wish to thank Professor B. R. Seth for his guidance ! 


the course of the preparation of this paper. My thanks are also due ' 


n\ x) P} -n(B)4 Zz B, Fi x) St,(B) . (9. 
1=0 


Professor L. Rosenhead and Mr. L. Sowerby for their helpful suggestion 
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ON THE DIFFUSION OF LOAD FROM A STIFFENER 
INTO A SHEET 


By W. T. KOITER (Laboratory of Applied Mechanics, 
Technological University, Delft, Netherlands) 


[Received 31 August 1954] 


SUMMARY 
A rigorous solution is obtained for the diffusion of load from a semi-infini 
stiffener into an infinite or a semi-infinite sheet. A Mellin transformation, appli 
to the basic singular integro-differential equation of the problem, yields a differer 
equation in a strip of the complex plane which is solved by means of Lapla 
transforms. The final solution is obtained in the form of an inverse Mellin integn 
which is evaluated by contour integration. 


1. Introduction 

THE diffusion of load from a stiffener into a sheet is a fundamental proble: 
of aircraft stress analysis. It occurs in many forms and it has been discusse 
by many authors, but so far very few if any rigorous solutions have bee! 
obtained. It is the object of the present paper to give the rigorous solutio 
of this problem for an infinite or semi-infinite sheet with a semi-infinil 
stiffener of constant cross-section. 

A semi-infinite sheet with a semi-infinite edge stiffener of constan' 
cross-section has been investigated by Buell (1) by means of a comple: 
stress function for the sheet. Although his solution is entirely adequat 
for all practical purposes, it is not rigorous because the convergence of his 
process of successive approximations to an infinite system of equations has 
not been established. Benscoter (2) has discussed an infinite sheet with 
finite stiffener, starting from the integro-differential equation for the shea! 
stress between stiffener and sheet. He observed that this equation i 
formally identical with Prandtl’s equation for the aerodynamic load distr 
bution over a wing of finite span, and therefore amenable to the sam 
(approximate) methods of solution. Our rigorous solution for a semi-infinit 
stiffener is also based on this equation—which, apart from a different scal 
factor, also applies to a semi-infinite sheet with an edge stiffener—but tht 
method of solution, i.e. the application of Mellin transforms, is quit 
different. It may be observed that our solution is likewise applicable t 
the aerodynamic load distribution over a semi-infinite wing. 


[Quart. Journ. Mech. and Applied Math., Vol. VIII, Pt. 2 (1955)] 
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). The integro-differential equation 


We consider an infinite or semi-infinite sheet (Figs. 1 and 2), loaded by 


| 


eo 














<—_— — 
Po > eel 
To (x) ) 
Fic. 1. Infinite sheet with semi-infinite stiffener. 
y} 














<-_—_ <«— =< a 
Uo (X) 
| 
PO > > > ~~ 
To (x) 
Fic. 2. Semi-infinite sheet with semi-infinite edge stiffener. 


system of continuously distributed forces over the positive x-axis and 


the negative x-axis: these loads are applied by the stiffener 


ind their reactions act on the stiffener together with the end load P,. The 


‘rect strain €, along the x-axis in the sheet in the case of an infinite sheet 
8S given by (2) 


-~ 
x 


(3—v)(1+v) [ 7(€) 
477 E J €—2z 


0 


dé, (1) 
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where the integral is defined as Cauchy’s principal value, and in the cay 
of a semi-infinite sheet by (3) 
=8 
2 T(E) 1é 
(€x)yag = ——z— | = S: 2 
aE J €—2 
0 

where E is Young’s modulus and v is Poisson’s ratio. The axial strain 
the stiffener is given by 
r 


P,—h | r9(é) ig}, 


0 





] 
a a 


8 8s 


where E, is Young’s modulus for the stiffener, A, is the area of the cros. 
section, and h is the sheet thickness. Equating expressions (1) and (3) ( 
(2) and (3) for a semi-infinite sheet) yields the integro-differential equatio 
for the shear stress 7,(x) between stiffener and sheet. 

The integro-differential equation may be written in non-dimension 
form by the following substitutions for an infinite sheet 


EA ee. ae 
a = 1(3—v)(1+v) 2-22. € = 1(3—v)(1+v) 282 
— 8 oF . | 
ee @—wil +) ZA, * 


the corresponding substitutions for a semi-infinite sheet are 
KA, - E,A, P, E 
i - 


Eh a Dee oF Ma 


uty 
ow 


“ 


Dropping the bars over the non-dimensional variables, the basic integr 
differential equation now reads 





- dE + | T(E) dé — 1 0. 


. 


0 0 
The additional requirement that all stresses vanish at infinity is expressed 
by the equation 


x 


| T(€) dé l. 


0 


3. Solution of the integro-differential equation 


We multiply both sides of (6) by x*-! and integrate from x = 0 tox =% 
obtaining 
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Formal inversion of the order of integration in the first term and integration 
hy parts of the second term results in 

sila — ‘ 
“ | gs-1- (£) dé '] 


— n— x7 9(2) dx 0. (9) 


; l—» 8 
0 Uv 0 


Introducing the Mellin transforms (4) 


( x-17,(x) dx, 


0 


(10) 


1 
dx, 
x 


(s ll 11 
H(s) | (11) 
0 


this equation is written in the form of a difference equation with variable 


efficients for the transform 7,(s) 


T\(s+1) = sH(s)Ty(s). (12) 


The foregoing formal analysis may be justified if the following assump- 


tions are made on the behaviour of 7,(): 


a) 7,(a)a4-1 is L(0, 00) for all A in a range b < A < 1, where 0 < b <1; 
b) \r9(¢ 7,(a)| < F(x)\E€—az)| for all |€—ax| < ex, where « < 1 is an 
whitrarily small positive number and F(x)a“ is L(0,0) for all A in 
the range b - A l 
7 ; 
v"| | T(E) dé lj} +0 and 2 [ tT)(€)| dé +0 
0 4 x 
for - oo and b \ l. 
These assumptions may be verified once the solution has been obtained. 
It is now obvious that integration by parts of the second term in (8) is 


justified by assumption in the strip b < res < 1. In order to justify 


the inversion of the order of integration in the first term we writet 


1 da | TAS? g, 
f r 
0 
{ l—é« A ril+e 
"pel B) ae - -t(E) 
im| | a-tde | <0! dg- | asl dx OS! dé 4 
E r 2 J e r 
0 0 r(1—e) 
{ i(l+e A x 
= , <a - — 
y= aa = dé 4 | x*-1 dx _ -- dé}, (13) 
, 1.8 d , gre | 
0 l+e 0 A(1+e) 
The a ; indebted to his colleague Professor 8. C. van Veen for helpful criticism 
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where 0 < « < 1. Bearing in mind that the inner integral in (8) is define 
as a principal value, the inner integral in the second term is also define 
as a principal value, and it follows from assumption (6) that the second 
term in (13) is bounded and of order « in the strip b < res < 1. The last 
term tends to zero for A +o in the strip b < res < 1 on account of the 
inequality 





); A a A o 
Dh a rr, (E 7 i dé 
| x81 dx | ie dé ee | x8-1| dx | T9(€) ya 
0 4(1+€) 0 1(1+€) 
J F — 
< Ac — Ae [ T(E) dé, (14 
A 


where o = res. In the two remaining terms the order of integration may 
be inverted (cf. Fig. 3). Putting x/é = », we write 


- a 

















SONS Ke 
s ae a , 
ee ad 
A(I-E) a(1+€) & 
Fic. 3. 
A x(1- €) A A(l+e) 
iin) | asl dx | Tol8) ae 4 | as] dx | T(E) dé\ 
Ao | cy é —&# . ry é—z J 
0 0 0 x(1+€) 
A(1—e) A A(1+e) E/(1+e) : 
: ° a a ; rast | 
= lim To(€) dé | —— dx 4 | T(E) dé | —— dx} 
4—a P F €—2x J ‘ c J é — 
0 €/(1—e) 0 0 
A(1—e) Alé A(1+e) 1/(1+e) 
: or i am a : = 
— lim £8-17,(£) dé | T dy + | é*-17,() dé | fie 
Avol J - ~ J l—y ; ie eaite l—y 





0 1/(1-e) 0 0 
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The limit of the second term may be written down immediately. In order 
to obtain the limit of the first term, we write, with 0 < B< A, 


i-e i 
oe ° ~t—! 
&8-17,(&) dé | dy 
é S tm 
0 1/(1—« 
Bil—« 1/€ 1(1—e) A/€ 
- get , A 1 
c lr (§ ) dé - dy + i? 17 ,(E) dé __ dn 
| / | 
1/(l—e B(1—e) 1/(1—e) 
Because A/€& > 1/(1—e), the inner integral exists, and it remains bounded 
rd4—>oo in the strip b < res < 1. Therefore 
1(l—e A/é 
; —_— 
, 7° 
lim | &8-1, (€) dé we dy 
io . J i-—7 
0 1/(1—e) 
Bil—e L 
P tg 7° 1 . 
€-175(E) dé | —— dn + remainder term, 
; | es, 
0 1/(1—€) 
re x x 
‘ se ¥ | nf 1 | 
remainder term &8-17 (€)| dé | dn. 
i J |l—n| 
Bil—e 1/(1—e) 


lf we now let B tend to infinity, the remainder term tends to zero. Hence 


ve may write (13) in the form 


TA 
Y 12 ol) ge 
& 2 
ss 1/(1+€) x ’ 
: ' | q 78-1 r 9st 
&8-17,(&) dé = dn + | f dn) +0(6) 
; | J 1l—y J 1-7 
0 0 1/(1—e) 
1—e x 
. : : bs 78 1 y n* 1 
| £*-179(€) dé —dy+ | = ~ dr + Ole). (16) 
| J l—y Ji-y ‘| 
0 0 l-+e 


hi 


is result holds for « arbitrarily small. If we let e tend to zero, we obtain, 


a 
. . . . “ 


s—l 
di &-17,(€) dé | —-— dy = Ti(s)H(s), (17) 


* . : a « . uaa 7) 


0 0 0 0 

vhere the integral H(s) is again defined by Cauchy’s principal value. 
In order to evaluate the integral H(s), we write 

1 € x 


. vf 1 * nf’ 1 7® 1 | 
H(s) ! d lim —dyn- —dn 
to ays fa 


De <0 “7 ‘ ” 
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and change the variable in the second integral into € L/n: 
i—< 1/(1+e) 
aS . — 
P 7 i. y 
H(s) lim! l - dn — | : dtl 
eo| J l—7 g teG 
0 0 
l—e 1/(1+e) 1 , 
i s-—l 8 7 8 s 
: 1 Y = Y —y7 
lim! | ies. dy — | ° * atl - | ain. dy. (18 
€ -o| ‘ | Ul P l —C 7 l Ul 
0 1—e 0 
Applying Legendre’s formula (5, p. 260) 
1 
- ~8—1 Ve 
2 ] dlog T(s 
| = dx = os ) y; (19 
J a—l ds 
0 
where y is the Euler—Mascheroni constant ( )5772157...), we obtair 
finally 4 ‘ 
. dlogT(s) dlogT(1—s 
H(s) — = e+ = 7 cot 7s, (20) 


ds ds 
holding for 0 < res < 1. 
The difference equation (12) may now be written in the explicit form 
T)(s+1) 2s cot zs T,(s), (21 


holding for b < res < 1. The solution of this equation has to satisfy the 


condition that 7)(s) is acceptable as the Mellin transform of 7,(z), ie 


T)(s) > 0 for |s| + 0 in the strip b << res <1. The solution obviously 


contains an arbitrary constant factor which may be determined from the 
‘boundary’ condition (7), viz. 7)(1) l. 
It may now be observed that 7;(s), which is regular in the strip 
b<res<l 
on account of assumption (a), is also regular in the wider strip 
6< ree < 2. 
This statement is proved by considering the inverse transform in the 
second term of (6), 


P . l Pe 
| To(€) dé — 1 seg | Ty(s+ l)a-s ds 
a7 
0 c 1x 
where 6 < ec < 1. On the other hand, we may write 
r r c+ia c+1a 
r t) dé i "Mm (eye-8 de I } lM (9\r-8+1 de 
T(E) AE = a o(S)E" ds = — — Ty(s)a as. 
H J 2m | 27 J s—l 
0 0 c 1x c ix 
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we may obtain the equation 





( c+ia 
j Le l — ' 
[ 7,(s) — 1 |a-8+! ds = T, (s+ l)a-8 ds 4 
dar S | 2m J 8 
€ iz 
c+1+ia 
l r l 
A. r-8+1 ds 
27 =, e—1 
c+1-—ia 
Writing s+- 1] s’ in the first integral in the right-hand member and s 8’ 
in the second integral, this equation becomes 
c 1+iz 
t ; ae 4 l ets 44 Jet 
T;,(s) r-s+l1 ds = | . Ty (8 )\—] Jar s+1 do’, 
S | Z 8 ] 
( l 1x 


From a general theorem? given by Titchmarsh (4, art. 9.9) it now follows 
| : ‘ , ’ . 
that 3 T,(s)—1] is regular in the strip ¢ < res <c+1 and the proof 


that 7,(s) is regular in the strip b < res < 2 has been completed. 

Finally, it may be stated here that the problem may also be attacked 
without recourse to the integro-differentiai equation (6). The state of 
stress in the sheet is governed by a stress function which satisfies the 
biharmonic equation. Applying the Mellin transform to this biharmonic 
equation and its boundary conditions ultimately results in the same differ- 

ce equation (21) for the Mellin transform of the shear stress between 
stiffener and sheet. This alternative approach may also be used in similar 
proble us where the kernel of the integro-differential equation is not known 
beforehand, e.g. for a wedge-shaped sheet with a semi-infinite stiffener. 


We may return to this more general problem in a later paper. 


4. Solution of difference equation (21) 
In order to simplify (21) we write 
~_ aen l 29 
Ty(s) = 2°T'(s) —— y(s). (22) 
sin $78 
Substitution into (21) vields the equation 
COS 778 ‘ 
y(s 1) ————. y(s) (23) 
cos 7s — | 
with the ‘boundary’ condition y(1) 1 The factor 2°T'(s) in (22) has been 
troduced for obvious reasons. The factor (sin }7s)-! serves to make the 
Titchmarsh formulates his theorem for Fourier transforms; its analogue for Mellin 


nsiorms tl nearly obvious 
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cofactor of y(s) in (23) tend to 1 as ims > +o in the strip b < res <] 
Taking the logarithmic derivative of (23), we obtain 

7 Sin 78 


a(s+1)—2(s) = = f(s), (24 


COs 78(COS 78— 1) 

l 8 

where x(s) = — dy(s) (25 
y(s) ds 








Equation (24) is much simpler than (21) because it is a difference equation 
with constant coefficients and it is solved by Laplace transforms.+ Assuming 
that x(o-+-ir)eH” is L(—oo, 00) in the strip b < o = res < 1 for all p in the 


range —7 < —d < yp <d < 7, we may introduce the transforms 
c -— 
X(w) | a(s)es” ds, (26 
c . ix 
c re x c os x 


. ™ sin 78 - 
F(w) = | f (sje ds = a es” ds, 27 


cos 78(cos 7s — 1) 





“. ° 
c ix c hee 


where b << ¢ < land —d <v =imw <d. Equation (24) is now trans. 


formed into (e-”—1)X(w) = F(w) 
: ° e~ 
with the solution X(w) = ——— 7 P(w). (28 
é — 
The solution of (24) is now obtained by means of the inverse transformation 


© oO 





- F(w)e-”8 dw, (29 


: Fw 
| ev—] 


x(s) = a | X(w)e-“*§ dw = — aa 


—C ~= @ 
where 6 < res <1. It may be observed that w= 0 is not a pole of 
(28) because 





c+io 
: d COS 778 
F(0) = | -| log ———_—— | ds = 0. 
J ds| ~coszs—1 
c—io - 


In fact, this is the reason why we introduced the factor (sin }zs)- in (22), 

The general solution of (24) is now obtained by adding to (29) any 
function with period 1. However, such a function cannot satisfy the 
simultaneous requirements that it tends to zero for |s| —> 00 in the strip 
6 < res < 1 and that it is regular in the strip b < res < 2, and it must 
therefore be omitted. 

The evaluation of F(w) (27) and x(s) (29) is performed by standard 
methods of the calculus of residues. The results are 

eter — 1) 
(e+ 1)(e"+ 1) 


+ The author is indebted to his collaborator Mr. J. B. Alblas for the suggestion that 4 


(30 





F(w) = —277i 


homogeneous difference equation may take a more convenient form by applying logarithm! 
differentiation. 
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nagreement with F'(0) 0, and 
7 m(3—2s) ‘ 
XS) rs A ES | (31) 
2 sin 78 sin 278 


[t is easily verified that x(s) is regular in the strip } < res < 2} and that 
t satisfies the assumptions made with the introduction of its Laplace 
transform. A solution of (21) with the ‘boundary’ condition 7)(1) = 1 is now 


given by 


ao 28-1]"(s) ' T m(3—2z 
Ts exp| | [> —+7 = " az), (32) 
sin $78 \2sinzz sin 27z | 
1 
where } < res < 2 and the line of integration lies in this strip. The 


exponential factor is bounded for |s|—> oo in this strip and its cofactor 
tends to zero exponentially. Our result therefore satisfies all requirements 
T,(8). 

[t should be noted that the solution of (21) with ‘boundary’ condition 
T}(1 
wltiplying (32) by any function Q(s) of period 1 which is regular in the 


stripb < res < 2 and equal to unity for s = 1. However, such a function 


s not unique. In fact, alternative solutions are obtained by 


) 


sat least of order e?7'*' for |s| > 00 in this strip, and the alternative solutions 
therefore violate the requirement 7,(s) > 0 for |s| > 0 in this strip. Hence 
32) is the only solution which satisfies all requirements on 7)(s). 

[he exponential factor in (32) may be expressed in Alexeiewsky’s G- 


nection, defined by (5, p. 264), 


G(z+-1) = (27)!e-1 nail [ {(14 2) e-stetan| (33) 
La) ey , 


n 
here y is again the Euler—Mascheroni constant. We write 
| T (3 22 ; 
| a 
2su sin 27z | 


(3 —2z)(cot 72 tan 72), dz 


{Sin 
\_* ) 3 1)\ 
Zt (2 l cot 7z-+-7 cot 72z-+- 27r(z— $ )cot 7( 2 1)\ dz 
{sin 7 7 * 
| ‘ 17 | COS 3 3 | 
2(z— 1)cot 2(z— 1) + 22(z— 3)cot 7(z— 3) dz 
2 sin 7 ” | 
l Ss 


rT oe az’ cot zz’ dz’ 4 nz cot 72” dz”, 
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where in the second integral we have written z— 1 z’ and in the third 
integral z—3 = 2”. Applying the formula (5, p. 264) 


s 


G Dans s 
| zz cot 7z dz = log ne yy 8 log 27, (34 
. G( 1 Ss) 
0 
we obtain after a simple reduction 
i a a(3—2z G(A G(s)G(3 = 
: | lit | dz log! sin lis (2) $) | (35 
\2sinzz  sin2zz | ‘ifs G(3)G(2 _ 8)G(—4 1 ¢)| 


The result (35) remains valid for all s in the complex s-plane except wher 
the line of integration in the left-hand member passes through or ends ir 


one of the poles of x(s) defined in (31), viz.s = Oors l+n(n = 2,3,4.. 
or s = $+k (k 1, 2, 3,...). Using the recurrence relation (5, p. 264) 
7(2+-1) P(z)G(z), (36 


expression (32) may now be simplified into 
2-1 G(s+1)G(§—s) ; 
_ (oi 


- | im 
Ty(s) 1(4) G(s—3)G(2—s)' 


It is now easily verified that (37) indeed satisfies (21) in virtue of (36) and 


the gamma function formula 


7 


sin 7z 


It may be observed that 7(s) is a meromorphic function because G(z+| 


is an integral function with n-tuple zeros z —n (n 1, 2,:3,...). Tm 
poles of 7)(s) are n-tuple poles s = n+1 (n 1, 2, 3,...) and k-tuple poles 
s = 3—k (k = 1, 2, 3....). 


The required non-dimensional shear stress distribution between the 
stiffener and the sheet is now obtained by means of the inverse Mellin 





transform 
c+ix c+ia 
I j l rP «28-1 G(s (3—8 
T(x) = - | T,(s)x 8 ds | ~ ( L)G( ; - Sds. (38 
27. 2m7t J (4) G(s— -1)G (2—s) 
where } <c¢ < 2. The axial load P(x) in the stiffener is given by 
P(a _— L f Ty 
(7) l- o(f) d ] | ol8) 4 Sds 
Fr. P, : 2a Pe l Ss 
0 c—ix 
P »s—1 Gis G 5 rs = 
1 ] | , r(s- r ] ) r( 2 8) yl) 8 ds. (39 
27i J T(4)(1—s) G(s—43)G(2—s) 
( ix 


where } and x is the non-dimensional coordinate. 
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5. Evaluation of solution 
Evaluation of (38) and (39) by numerical integration is cumbersome. 
However, it will be shown that the line of integration may be closed in 


the left half-plane re s c if the contour does not pass through any pole 


97) We take a rectangular contour with sides re s c,res k+ 3, 
chere k is a large integer, and ims tr, where t+ 0. It is easily seen 


fom the bounds on the exponential factor in (32) that the contributions 


fthe sides im s + tend to zero for 7» o. Furthermore, we write 


T 28 1[*(s) Gis)\G(2 8) 
. (4) G(2—s)G(s—}) 
s—1 
2 *T (a) . 
We exp) | 2zcot 7z dz + (s—1)log 27+ 
(3) . 
” 0 
| wzcot 7z dz — (s 3) log 2| 
0 
: s—1 
2°-*}'(a) | 5 | 
ray) °xP | mzcot mz dz-+ 4 log 277), (40) 
(5) \ ‘ = 
where the exponential factor is seen to be O(e*) on the line res k+ 2. 


[tis now obvious from the asymptotic expansion of the gamma function 
that the contribution of the side re s k+-? of the rectangular contour 
tends to zero tor k > OO. 


rhe solution may now be written in the form of the series 





T(x) > residues of [7)(s)a~*] in s k+3, (41) 
K I 
Pia . ° > [Ts ° « 
l SS residues of ol$) 11 8s) ins —k+3, (42) 
A ang l—s a 
where 7)(s) is given by (37), and the series are convergent for all x > 0. 
Because s k+-3isak-tuple pole, these series may be written in the form 
‘ ] dk-1 - 
ry = SY | 4 [st k— 3k T(8)x—]| (43) 
a k | } \ dish 1 7 e k+38 
} l 
P ‘ l [k-1 T,(s 
a fo | (s+ e—3y L008) gre] | . (44) 
i a~ (k—1)!|ds*—! “ I—s = 


The first few terms of the series (43) and (44) may be calculated in ; 
straightforward manner, using the formulae 


f d | f(s 
7 . og] (Ss), 


a? fd. fd, i)? 
P( Ss) i(s) log (s)-—+ log f(s ’ tC. 
f og f(s) \ds og) (8); | etc 


1s? | ds? 
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By means of the recurrence relation (36) the final result for the stiffeng 
load is obtained after some elementary but tedious algebra: 


iy Dy . dy 
Poe) 1/2) 2 /(®):t0g24 40) 
0 \T] on m1 


-{2-+- log 2+-5($)—log #}"] +, 


S| %e 
_— 
a1 
— 
| 
cy 

to 
| 
e 
poor 
+ 


{ ( 2ac\ 
1— /| } [1 —-2(0-25425 0-10610 log a) + 


+-x*{0-0086265 — 0-018888 log x+-0-0033774(log x)?}+-...], (45 


where ¢(z) is the logarithmic derivative of the gamma function, and dasha 
denote differentiations. 

The series in (45) is rapidly convergent when z is not too large, say x < 
as may be seen from Fig. 4, where the curves a, b, c show the result by 
taking consecutively 1, 2, or 3 terms of the series between brackets in (45 
It should be noted that between 2 around 0-5 and x around 1-5 the curv 
c lies below curve b, the maximum difference being approximately 2 per 
cent., and that from a value of x slightly above 1-5 curve ¢ lies above curv 
b, the difference becoming important for 2 > 2. For large x the conver 
gence is slow and the series (43) and (44) are not convenient for numerica 
computations. In fact, it appears that for 2 between, say, 2-5 and 5 th 
approximation by three terms between brackets in (45), curve c, shows 
even larger errors than the approximation by two terms, curve b. 

However, a useful approximation for large x is obtained from the asymp- 
totic expansion, deduced from (39) by contour integration in the half-plaa 
res >c. The first few terms of this expansion are readily written dow 
as the residues in the simple poles s = 1 and s = 2 and in the double pok 
s = 3. The resulting asymptotic approximation for large x is 


P(x) 2 4 
: : : 5) ; (log 2 y—log x) sic (46 
re a as © 


1.0¢ 
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oa 
~< 
—) 


>| 
o 





0.8} 


0.6 


0.4 





This result, retaining one or two terms in (46), is also depicted in Fig. 4 
(curves d and e), and it is obvious that a satisfactory result for all x > 
is obtained by linking the convergent expansion for small and moderate 2 
with the asymptotic expansion for large x. This link has been drawn in 
Fig. 4 (curve f); the maximum possible error is estimated to be at most 4 
few per cent. 

A comparison with Buell’s results shows very satisfactory agreement. 
Buell’s curve for the stiffener load agrees with curve f in Fig. 4 within a 
few per cent. Our solution may serve to judge whether the number of 
unknown coefficients b, /P (notation of (1)) retained by Buell is adequate. 
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10 —-— 
p(y) 
aL! a 1 TERM BETWEEN BRACKETS IN (45) 
\ b : 2 TERMS BETWEEN BRACKETS IN (45) 
C 3 TERMS BETWEEN BRACKETS IN (45) 
d 1 TERM IN (46) 
, e 2 TERMS IN (46) 
0.8) f : ACTUAL CURVE 
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Fic. 4. The stiffener load as a function of the non-dimensional coordinate. 


According to his solution the first two terms corresponding with our 


expansion (45) are (1. equation (42)) 


P(a 2 2~~nzb,] (47) 
VI -- a ’ 
Fs w he 
’ n=1 I 


julvalence of (47) with the first two terms in (45) requires 


> (4) —1 = 0-253314. (48) 


n l 
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Buell has solved his infinite system of equations approximately by retaining 
at most six equations and the first six coefficients. His result yields 


+. nab 
3 a” — 0247108. 
P 


n=1 
indicating an error in the verification of (48) of around 2-5 per cent. 
Finally it may be remarked that the numerical work required fo 
evaluation of the present rigorous solution is negligible in contrast t 
Buell’s approximate solution. Howevei, this reduction in numerical con. 
putations has been bought at the expense of considerably more analytical 
work, involving more advanced methods. 


Note added § March 1955 
A. Pfliiger (‘Halbscheibe mit Randglied’, Zeitschr. f. angew. Math. u. 
Mech. 25/27 (1947), 177) has also observed, independently from Benscoter, 
that the basic integro-differential equation is formally identical with 
Prandtl’s equation. 
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SUMMARY 

Scattering and transmission of sound waves of somewhat more general form than 
nle plane waves are analysed for scatterers in the form of a thin elastic simply- 
ported rectangular plate in a rigid infinite baffle. The plate is assumed to separate 
| different fluid media, and the sound radiation from the plate undergoing forced 
rations under the influence of the incident wave pressure is studied. Unlike the 
| scatterer, whose dynamic characteristics are little different in vacuum and in 
rtain fluid media, the characteristics of elastic plates are profoundly altered by 
uid reaction. By correlating the general solution of the wave equation in the fluid 
lia with that of the Lagrange equations of the elastic plate at the plate—fluid 
lia interfaces, the deflexion of the plate motion is determined. A transmission 
ent 7 is then defined and derived. In the appendixes, expressions for an energy 


y and the reactions at the plate supports are given. 
|, Introduction 
luE scattered and transmitted patterns which result from the incidence of 
pressure wave on an elastic plate separating two fluid media may markedly 
liffer from the patterns obtained from a rigid scatterer of the same form. 
This difference, recently analysed for the case of an elastic shell (1), is due 
the fact that the alternating sound pressure of the incident wave excites 
reed vibrations of the plate. The nature of these vibrations is determined 
y the dynamic characteristics of the plate separating the two fluid media. 
these characteristics, distinct from those of a plate vibrating in vacuum, 
lepend on the density, elastic constants, dimensions, and boundary con- 
tions of the plate, and on the densities of and the sound velocities in the 
two adjacent fluid media. It is desired to find the sound pressures for the 
ittered and transmitted waves in terms of these characteristic parameters 
fthe plate and the fluid media. 
The following assumptions are made: (1) the fluid media satisfy the 
nditions for the validity of the simplified wave equation; (2) the incident 


we pressure is of the form 


p(x, y,2z,t) = Cig exp t(—a, x—a, y)exp tw(t —z/c), 
ere w, d,, and a, are such that w/c, > a, and w/c, > a,; (3) the plate is 
1his pay lraws on research carried out for the Bureau of Ships, Department of the 
, U.S.A., by Reed Research, Inc 
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composed of a material which is isotropic and devoid of damping, ay 
which obeys Hooke’s law; (4) the dynamic deflexions and the thicknes 
of the plate are small compared to its other dimensions; and (5) the plat 
boundaries are simply supported. 

The standard methods of the theory of mechanical vibrations are ug 
in the analysis. The dynamic configuration of the plate is expressed i: 
terms of an infinite series (Fourier) whose coefficients are the generalize 
coordinates of the system. The potential and kinetic energies of th 
deformed plate are then calculated. The sound pressure, in the form ¢ 
the general solution of the wave equation in the fluid media, is integrate/ 
over the surface of the plate to obtain the generalized forces of the plate 
The Lagrange equations for the system are then introduced. Another s 
of equations is obtained from the continuity conditions, equating th: 
normal component of the plate mid-surface displacement to the norma 
component of the fluid particle displacements at the plate interfaces 
Simultaneous solutions of these two sets of equations yield an expressiol 
for the generalized coordinates. With these a ratio between the incident 
and the transmitted rates of energy, called ‘transmission coefficient’, i 
determined. The series expression of this ratio is obtained in terms 0! 
densities and acoustical velocities of both fluid media, the plate pane 
dimensions, area density, and flexural rigidity. Derivations of anothe 
energy ratio yn, and expressions for support reactions are also included. 


2. The total excess pressure (upon the plate) 
The pressure wave equation for an undamped system and no sources iI 
the region under consideration is 
a9 
. 9 o-p 
cV%p = pe (] 
49 
ot" 
where p is the sound pressure, ¢ is the time, and c, is the velocity of sount 
in fluid medium 1. 
In Cartesian coordinates, the complete solution of (1) is given (2) by 
p = 2 > > CryzeXP i( +a, t+, y ta, zt, t), a’ | a; T a? : Wy q 
ds Gy Gs 5 : : 
where for every set of values of a,, a,, and a, there will be a discrete value 
, ; . ‘ : cm 
of angular frequency which is designated w,, in (2). If we assume simp! 
harmonic motion of e+‘! type of dependence on time, only one set of values 
A,, A,, A, = 4, dy, dz remains, with one related frequency w,, = «, so that 
we can write 
9 


a 


C ey ot » I 2.1 g2 3 ,,2/¢2 | 
p = Cygexpt(ta,rta,y+a,z+oat), ay+a5-+-ag = w*/Cy, \ 
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rin abbreviated form 


p p(x, y, z)exp(twt) Cio3 X (x)¥,,(y)Z,(z)exp(twt). (2 b) 


Yewton’s equation of motion gives a relation between the particle 


velocity vector U and the pressure, 


cu 


e grad p, (3) 
Cc 


} 
shere p is the fluid density. 
The particle displacement vector corresponding to the solution for a 


simple harmonic wave of the type (2b) is 


u = u,(z, y, z)exp(iwt), (4) 
and by using (3) the velocity vector may be written as 
' 1 ou - 
u : grad p. (5) 
iw ot twp 
For cases where w/c, > a, and w/c, > dg, u,, and w,, are small compared 


with #, and may be neglected. We are then only interested in the z- 
mponent of u, which for fluid medium 1 is 


X,, (a \¥,(Y) C Z,,,(z) 


u exp(iwt). (6) 
twp, C2 
nm the interface between medium | and the plate z th = 0 we can 
write for the particular velocity 
X,, (x)Y,,.(y) . ,,(0Z,, (2) ” 
uw) 0 - od Tae exp(twt) i ’ (7) 
Ww py \ C2 2=0 
Tr, using (20 
p,(x,y,0)[0Z, (z) /, : 
ij i Pr | Z,,,(2) exp(iwt). (8) 
lwp, C2 / _ 





We limit now the nature of the scattered wave pressure in medium | to 


simple harmonic form 
pi (x, y, z)exp(iwt) X',, (x) ¥;,,(y)Z,,,(zexp(twt), (9) 
vhich according to (5) produces a velocity along the z-axis equal to 
7 ‘ tT > \Y) OBA) ae atind, (10) 
twp, Cz 
vhere dw, which are small for small plate deflexions, are assumed 


7 


e compared with u 
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At the interface z = —}h = 0 we have 
(2, y, 0 ’ ; 
(i, ).-0 — Pil® ¥, | En [Zn exp(iwt). (1) 
—tWp; Cz “nd 


nature to a simple harmonic form 
Po = p(X, y,z)exp(+iwt) = X,,(x)¥,,(y)Z,,(z)exp(iwt), (12 
which according to (5) renders a velocity along the z-axis 


— X,,(x)¥,,(y) € OZ ,,(z ) 


a 1W Ps Oz 





exp(twt). (13 





interface z = +4h = 0 we get 
: o(x, y, 0) [OZ [ , 
(a,)..9 = Pal. y ~: 22) /2,.0)] exp(twt). (14 
ce —tUW Po C2 ; _ 


The condition of continuity requires at the plate interfaces 
(%,,)e-0 + (Uz,)z p= FF = (%.,)2-0 (15 


(w.,). . ne (16 


or 


(%,.)..9 = @ 


where w is the normal displacement of the plate. 
The total excess normal pressure on the plate is 


[ p(x, y, 2, t)]._, = p(x, y, 0, t) 
= p,(x, y, O)exp(iwt)+-p;(x, y, O)exp(iwt)—p,(x, y, O)exp(iwt). (17 


Using (8), (11), and (14), however, (17) reduces to the form 


; ‘ oZ , 
p(x, y,0,t) = =e ‘ps\@ | (t, ).<9— 
5 z=0 


— inps| Zn(2) ) Q Z,(2) _ (u,,)e-0+ lwp, Z,,(2) o2,,{2) Fee” ite 
/ oz 2=0 Co JSe=0 
{ s 


and using (16) we obtain 


: oZ a OZ, (z ‘ 
(2, ¥,0,t) = -iwpy| Z,,(2 2) | nil =) Zip(e) |B) (%,,)e-0— 
/ 2=0 


Oz 


ie os Zp [= a Zp(2) |e) w. (19 
=0 


21 
Cz 


nw 


For the transmitted wave pressure in fluid medium 2 we again confine it 


Here w,, and wz,, are also assumed negligible compared with w,,. At the 
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Using (2a) we find 


, Ww ° ° ) 
Z (2) exp! d (= —a}—a3)| = exp(—i22) 
N \Cy } Cj 
_ Ww . m ; 
Z,,,(2) exp! 2 (S, — aj—ai) = exp(+! —2Z]} }, (20) 
Ni \Cy } Cy 
: w? ° ‘ -W 
Z,,,(2) exp! LZ (= - aj—a; | = exp|—1—2z 
pe J \c2 | é& } | 
; . 1 az ag\-}) 
where a > = - 
Cj w* w* 
(20a) 
ind C (- a | ay , 
s cs w= w* ] 
We then have 
C Z, (Z) @ a Ss 
— = Ly (2) 
¢ Cy 
OZ, (z ee 
Pi ) L i Z', (z) c (20 b) 
re Cy , } 
oZ m4 -W, 
= - U ~ Z (2) 
Cz Co 


which when substituted in (19) give finally the total excess pressure in the 


(21) 


lorm 


p(x, y,0,t) = 2p, €\(U,,)--o—(P1 © + p2€2)e, 


r, using (8) and (20), 


p(x, y, 0, t) 2p,(x, y, O)exp(twt) —(p, €,+ po Fo)w. (22) 


It can be noted that, under the assumptions in this analysis, we always have 


) 
yt) aq 


; and (w/c)? > a} in both fluid media, which results in ¢, = c, 


ind Ee 


; so that (22) may be approximately written as follows: 


p(x, y, 0, t) 2p,(x, y, O)exp(twt) —(p, C+ poCe)w. (23) 


3. The Lagrange equations and their solutions 
Consider the plane 0 at the middle surface of the undeflected flat 
plate of thickness h with its edges at x = 0, y = 0,2 = a, andy = b. For 


small deflexions of a thin plate the effect of the static deflexions upon the 


dynamic deflexions may be neglected. The equation of motion for the trans- 


verse dynamic deflexions of the plate is 


DV4w+-mit (24) 


p(x, y, 9,t), 


where D is the flexural rigidity of the plate. 
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For a simply supported plate. w(x, y,t) can be expressed in the form of 


a double Fourier series 


L x 
a . Mat . nary 
w(a.y.t) = SS dnp(t)sin ——sin—=, (25) 
—— od 


m=1n=1 
where the coefficients g,,,,(f) are to be determined as solutions of the 
Lagrange equations 
Gigs ' oT oV ma 
a T a Qin (<0 
dt CFU mn 


c Ym ” CO mn 


and are the generalized coordinates. In (26) the Q,,,, are the generalized 


‘tmn 


forces. The potential energy or strain energy V for the plate is 


a b 
7 Df f (@w . ew)? 
J | ue oe | dady. (27 
2) hee ey?| 
0 0 
Substituting (25) into (27) we have 
o abm~s m ,ls 
— 2 e ) 
J sets ' S S Wmn YUmn> {<5 
Pe) dowd 
m=1 n=1 


where w,,,,, is the natural frequency of the mnth mode of the plate vibration 


? 


The kinetic energy 7 of the vibrating plate is given by 


aiob 
- Of TT 
7 ua dudy, 29 
y 
0 0 
which on using (25) becomes 
-_ x x 
a. a = — 4 
7 ‘ > © ie (30 
8 4 A Umn> 
« as — 
m=1n=1 


where m is the surface density of the plate. 
The mnth generalized force may be written as 
a b 


Pr . Marx . Nm s 
4 | | p(x, y, O, t)sin —— sin Y dady. (31 
2 4 a b 
0 0 
Substituting (23) and (25) into (31), we have 
a b 
: am i . Mrx . nm 
Onn = Z2exp(iwt) | p(x, y, O)sin sin Y dady — 
‘oe a b 
0 0 


_ a) 
P11 T P22 4 1mn°* 


Expressing now p,(x,¥,0) by a double Fourier series of the form 


py (a, y, 0) >_> d,.8in 
Se ( 


r=1 s=1 





we get 


Hence 


Substi 
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- form a vb } 
. Mnx . nay ab . 
get | px (x, y, 0)sin sin“ dady = (34) 
se a b 4 
0 0 
‘ fence (32) gives 
ab . ab " — 
~ ot the Gm 9 A exp(twt) (py Cy 1 P2C) 4 Gmn- (35) 
Substituting (35), (28), and (30) into (26), we finally obtain the Lagrange 
{ uations in terms of the generalized coordinates 
. . . = 2 2 Dr 
eralize Py Poe . Dna m n <-Ann : ° on 
q - q = 9 | 9 Ymn = exp(twt). (36) 
m m \a2 6b m 
in be seen that the transient solutions will die out and the steady-state 
97 solutions of (36) are as follows: 
exp(iwt) _ 
Gm . r, — (37) 
wn» 
ere Z,,,,, Which is of the nature of a mechanical impedance, is given by 
: m | iw 
Z ; | (we w* —(p, Cy P22) 
ora an tWi i 
- NIG 
1 ; 2 ‘ 
2 py Cy 1 P2Cy) ae (Linn 1). (38) 
9 rresponding to the imaginary part of Pi; i.e. 
im p, Cho3 COS(A, X+-Ay y)siN wt —C}o, sin(a, +A, y)cos wt, 
7) may be written as 
q Imn \COS(wt+ yn + mn), (39) 
Z\A, , {x 1\2]-3 
here | 14 B2.n( mn - (39a) 
W(P1 Cy +> P22) | Ly 
: . 
p n\n 1) ‘ 
tan : (39 b) 
b 4 
mit 
im A 
X tan-! =. (39 c) 
rea 
= mn 
whicl ns 
> ai 
t ») WwW 
7 mi 
) P17 P22 
D | m- n> 
n(—\*( +55}. (40) 
\7 a* b> 
4. The transmission coefficient + 
Le Os 1 } . . . . . 
We define, after Vogel (3), the transmission coefficient 7 as the ratio of 
ean rate of energy transmitted to the second fluid medium at the inner 

















186 J. MARTINEK AND G. C. K. YEH 





interface 2 over mean rate of incident energy in the first medium arriving 
at the outer interface 1: Wee 
°= iT.’ (41) 
1 





— 


According to reference (3), the mean rate of vibratory energy of the incident 
wave arriving at the interface over a period J, can be expressed (under 
the assumption that w,, and w,, are negligible compared with %,,) as 
To ab 
Wie = - | dt [ [3 im[ p,(x, y, 0, t)]im[a, (x, y,0,t)|dady. (42) 
0 Jv @¢ 
0 0 


From (33) we have 
— 1, me 
im[ p,(x, y, 0, t)| = > > |A,,{Sin sin =F" sin(a,,+at), (43) 
and from (8) and (20b) 
er Ls 
im|{%,,(x, y, 0, t)] = —— im[p,(z, y, 0, 4)]. (44) 
Pi“ 
Substituting (33), (43), and (44) in (42) and simplifying, we get 


— l 
V i —_ : : b > ' mn|- as (45 
16p, Cy 


m 


Similarly the mean rate of transmitted energy at the interface 2 over a 


period 7, is T% a 
Wi> — , 4) dt | { sim|[ pp(x, y, 0, t)]im[w,,(x,y,0,t)|dady. (46 
a 5 ny 
00 
Substituting (14) “_ (20) into (46), we obtain 
To a b 
Wi = Z| dt | | ¥ 1 by Cy W?(a, y, t) dady. (47) 


By virtue of (25) and (39) and carrying out the integrations, 


ab “ < 

77 ¢ — 1 9 = 

Vee wan 16 P20” > > IQmn|*- (47a 
m=1n=1 


Substituting (45) and (47a) into (41) gives 


Wee (pyc SS as 
+ ile — (pr WP2 c: 2) 2 > 1Gmn\'- (48 
1 S > Amn |” m=1 n=1 


m Ss n 


_— 


Using the expression (39a) for |q,,,,|, we finally obtain 


_ 4(p1¢1)(p2¢2) >> Ann |? (49 
c oo 9 ¢.)* “ = 1+ mn 1)/x. , 
™ ” ’ = 17 2 * nmi m=" Ft sa 


=1 
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where f,,,, and 2,,,, are defined in (40) and A 


mn 


is given (see (34)) by 


mn 


Mists = re Am n +t im Rue ? 


in which 


4mnn~( 123 


reA,, — > a SEY PM ATs as 
| (mar)? — (a, a)* || (nar)? —(a,, 6)? | 
<{(—1)™+*"cos(a, a+-a, b)—(—1)™cos a, a—(—1)"cosa,b+ 1], (50 a) 
4mn7r*C,,. 
. w 123 
imA,,, --— <= 


[(mm)®— (a, a) ][ (n7)®— (a 6)*] * 


<{(—1)"+"sin(a,a+a,b)—(—1)"sina,a—(—1)"sina,b]. (50b) 


5. Example on evaluation of r 
Let us consider a plane incident wave p,(x, y,0,t) = posinwt. We find, 


from (50), 





a t (16p,/7?2mn for m,n odd, 
4p, ( [{ .. marx. nary 
Ans : | sin sin —~ dady = . . : 
. , ab. - b 0 for other combina- 
0 0 tions of m and n. 
(51) 
Noting that 
136... 1.26 1,3,5... 1,3,5... 1,3,5... 9 9) 9 
a | ey Ee 1 Lar? 2\ (m2) a8 
; Ss S Sle) ee) 7, 
hat ht SM? Law m* Ly n* m*\8 8/\8 64 
m n nm m 
(52) 
we obtain from (49) and (51) 
12%... 135.- 
= 644 (p1C,)(P2 Co) i 9 I (53) 


; . \2 . tint.) OR Siam a 
(pP1¢4 T P2 &2) a oi men | l+ Beent (Teen 1) XLmns ] 
Ifthe two media are the same, p,C, = poCo, and if we use the approximation 

Pi % P22 


7 = 10, we get 


PSO NS 0-64 (54) 
T ’ v0 
—1 hats M*N0* [1+ B ‘(x2 .— 1)/a 32] 


mnt 





mns 


m n 


which is the result Vogel obtained in reference (3) (equation (19)). 
Vogel further showed that, for a medium of low density (e.g. air) and for 
a source of high frequency (e.g. sound source), (54) can be approximated 


by a finite inequality as follows: 


or 
cr 
— 


0-32 ., ‘32 |, 2 , 
. K T = K+ Aprea)". kK |e (i 


abw abw mw" m 
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By similar reasoning it may be shown that, for the case of low density and 
high frequency, (53) can be approximated as follows: 
(=) tee § P ., 
7 (P10, + pole)? abwal m — 


(=) (pr C1)(P2¢2) [P ., 4(p1C1)(P2 2) (56) 
(P1 C1 + P22) 





ter 


abw A m mw? 


2 
7 


6. Energy ratio and reaction forces 
If energy of the vibrating plate rather than the transmitted energy is of 
interest, a new coefficient called ‘energy ratio’» may be defined as 
mean elastic energy absorbed by the plate 


i=— = —— ——_—— 


mean incident energy in the first medium at the interface 1° 





This ratio is derived in Appendix A. It is an infinite series and may also 
be expressed in terms of 7. 

When the dynamic deflexion of the plate is determined, the periodic 
reaction forces at the edge supports R,, R, can be easily evaluated from 
well-known formulae in the theory of plates. They are derived in Ap- 
pendix B. 

While the transmission coefficient 7 is of interest from the point of view 
of acoustics, both energy ratio » and reaction forces R, and R,, are of more 
immediate concern from the point of view of elasticity. 


7. Conclusions 

The analysis gives expressions for the plate deformation w and the trans- 
mission coefficient 7 in terms of various parameters of the plate and the two 
adjacent fluid media. In Appendixes A and B the expressions for the energ) 
ratio 7 and the support reactions R, and R, are derived. 

The results apply to cases where w/c >a, and w/c > a,. This includes 
problems where: (1) the velocities of sound in the fluid media are low and 
the frequency of the incident wave is high, i.e. the wavelengths are so small 
that they may be neglected as compared with 2z/a, and 2z/a,; and (2) the 
incident, reflected, and transmitted waves are so close to a plane wave that 
a, = Oanda, = 0 in both media. 

It is interesting to note in extending the work of Vogel to the more general 
case of two different fluid media adjacent to the plate and to a somewhat 
more general form of the incident wave that such quantities as 7, 7, R,, and 
R,, are still representable in relatively simple algebraic forms. 

The authors believe that the present analysis can be extended to plates 
with other boundary conditions, as well as stiffened plates and plates wit! 
structural damping. Further investigation with still more realistic inc 
dent, reflected, and transmitted wave forms will be rewarding. 
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APPENDIX A 
The Energy Ratio » 
We define a new coefficient 7 called the energy ratio as follows: 


mean elast nergy absorbed by the plate ETo 


an incident ¢« in the first medium at the interface 1 ET" 





(Al) 


wn ineident energy at the interface 1 over one period can be obtained by 


tegrating the mean rate of incident energy given in (42) 


1 b 

E} Ws dt ; | it] it | Lim[p,(a, y, 0, ¢)} im[e.,(2, y, 9, t)] dxdy. 

0 ‘ par (A2) 

I 33 43), a 44 
| =e 7 a 1] sin 20») —- S S | 
Serax’ = mre. So ant 
(for dm» = 9). (A3) 

wder 4 mpute the 1 in elastic energy absorbed by the plate over one period, 
lefir Fi 3 8 


i To , | tim pla, y, 0, t)| dw | | dady, 
rr 0 0 
1 db 
} : | | t] lt | im{ p(2, y, 0, t) w(x, y, t) dxdy. 
00 


5), and (39) in (A5), we find 


x x 
ab ~— — wabr( py Cy + P22) a {~~ P ‘ 
Ln Le , ya a 


S 16 
m=ln=1 
1} 39a) Ww 
I ( ~~ Amn |” iE 5 (= yy] 
i = > f nr 
1+ Beant (22m —1)/?mnd? | mn 
Substitut f (A3) and (A7) into (Al) renders the expression of the energy r 
{ , | Pe 
Ss S \ A 1 Dmni\’mn 1) Ly a 
1 Fee: Linn — 1\*]8 | 
| Jmn — | 
J nn 
\8) with (49), we may express 7 in terms of 7 as 


(A4) 


(A5) 


(A6) 
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(A7) 


atio 7) 


(A8) 


YoY Amn EP a Pee 


Q3 (42 1274 . 
— (1 Pmni("mn 1) Ymns J” P22 
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APPENDIX B 
Reaction Forces at the Plate Edges 


The vertical reactions on the supports of a simply supported plate due to the 
dynamic loading (excess pressures only) can be expressed as 


a3. 3. 

Ow Ow 
R, p| 2% +(2- pL) =| along2x=—0O and x=a (Bl) 

Cx" Cxoy? 

. Ow oew 
an v, ere eet L) a aiong y an y ), (B2) 
and r. D -+(2—; . l { 0 4 U 39 

F cy CYCxr* 


where y is the Poisson’s ratio of the material of the plate. 
By virtue of (25) 


x x 

—, im\3 m marx . niry , 

R, = Dr? > Zz [(”) +(2 »)(™ =)(¢ )" Jamn eos Pa b- (B3) 
m=1n=1 


Hence, after substituting (39) into (B3), 


Re Dr* > 6 (“\[(*)'+0 1)(%2) |sin 2 








b 
m=l1n=1 : 
Ymn Cos(wt - Pmn r Xmn) at x 0, | 
and | (B4) 
m sd | : n\2] . nary 
R. | (4) ‘ (") |sin” ry 
pm > > y(° H) b ‘ b 
m=1n 
YImn|COs(wt + dmn TO&mn) atx a; 
and similarly 
x cw . ‘ Ps 
nN n\* m marx 
y= —ae SF (*)f(%)*4c2—p9(%) Ian? 
> > GG) +e-o (2) ]sin™ 
m=1 n=1 
Ymn|COSs(we + Pinn TAmn) aty 0, 
and . (BS) 
Max 
> Dr’ (Ss F 
Ry = De YY (G+ 2-02) nt 


Amn cos(wt + dmn t Xmn) at Y¥ b. 


The values of R, and R, sa derived represent the periodic forces per unit edge length 
acting on the supports at edges x 0, 2 a and y 0, y b respectively, due 
to the excess pressures. Actual reactions can be obtained by superimposing these 
values on the reactions due to the static loading, if any. 
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A BOUNDARY LAYER PHENOMENON IN THE 
LARGE DEFLEXION OF THIN PLATES 
By Y. C. FUNG? and W. H. WITTRICK{ 
Received 3 February 1954] 


SUMMARY 
In the large deflexion of thin plates with free edges, boundary layers develop along 
se edges under certain conditions. The non-linear effects of membrane stress are 
nfined to the boundary layers and elsewhere the equations are linear. Boundary- 
yer equations are derived and integrated, and are then applied to some specific 


blems. 


|. Introduction 

\slong ago as 1912 H. Reissner (1) pointed out the existence of a boundary- 
yer phenomenon in the bending of thin shells. Suppose, for example, that 
the shell is subjected to a pressure on its surface together with tractions and 
bending moments on its edge. Then, if the thickness is small enough, the 
leformation can be determined approximately by neglecting the flexural 
rigidity of the shell and assuming that the applied loads are resisted entirely 
ythe membrane stresses. This approximation breaks down, however, near 
the edge of the shell where, in a narrow region, or boundary layer, the 
leformation changes rapidly to satisfy the conditions on the edge. 

A different type of edge effect was described (2) by ‘\elvin and Tait in 
1867 in order to give a physical explanation to Kirchhoff’s boundary 
nditions for an edge which is subjected to shear force, twisting moment, 
and bending moment. Using variational methods Kirchhoff had previously 
shown that the two conditions involving shear force and twisting moment 
combine to give only one condition. In explanation Kelvin and Tait pointed 
ut that twisting moments of line intensity /,,, acting on the edge are 
statically equivalent to shear forces of intensity @,,,/és, where s is the 
listance along the edge. The combination of shear forces V,, and twisting 
moments J. on the edge are therefore statically equivalent to shear forces 
1 intensity (V,+0M,,,/és). On applying Saint-Venant’s principle it follows 
that if a solution is obtained which gives the correct values for the bending 
moments and these equivalent shear forces on the edge, it will give an 
iccurate description of the state of stress and deformation in the plate 
except in a very narrow region near the edge. This edge region, or boundary 
iyer, was discussed by Friedrichs (3) in 1949. Friedrichs showed that the 
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width of the boundary layer is of the order of the thickness of the plat, 
a fact which is of importance in the work to be described later. 

The present paper is concerned with yet another type of boundary-laye 
phenomenon. This may best be illustrated by a consideration of an initially 
flat long strip, of width 2b and thickness ¢, bent uniformly in the longituding| 
direction into an are of a circle of radius R. 

It is well known that simple bending theory when applied to such a striy 
involves an anticlastic curvature equal to —/R in the direction of the 
width of the strip, where p is Poisson’s ratio. Experiments show, however 
that this anticlastic curvature does not necessarily occur; in fact with very 
thin plates the cross-section remains practically undistorted. This pheno- 
menon was discussed by Searle (4) and by Case (5) who both concluded 
that the cross-section distortion depends upon the value of a parameter 
b*/ Rt. If b?/ Rt < 1 the anticlastic curvature corresponding to simple bend: 
ing theory develops, but if b?/ Rt becomes large the cross-section remains 
almost undistorted. The transition stage between the small and large 
values of the parameter has recently been investigated by Ashwell, both 
theoretically (6, 8) and experimentally (7). He showed that as b? Ri 
increases the distortion of the cross-section becomes more and more con- 
fined to a narrow region or boundary layer adjacent to each edge of the 
strip. The membrane stresses which exist in the strip are also confined to 
the boundary layers. Thus when the longitudinal curvature is large the 
strip assumes on the whole a cylindrical form, with the generators perpen- 
dicular to the two edges, and in a state of pure bending. Disturbances from 
this state occur in narrow boundary layers along the edges, where the middle 
surface, being no longer developable, becomes stretched and membrane 
stresses appear. The width of the boundary layers is of the order +(Ri). 

It is clear that throughout the cylindrical portion there must exist in 
the plate a uniform lateral bending moment equal to u times the longitudinal 
bending moment, in order that the lateral curvature should be zero. This 
lateral bending moment is obviously zero on each edge and it therefore 
increases through the width of the boundary layer from zero at the edge 
to its required value inside. This is the important function of the boundary 
layer, and the mechanism by which it accomplishes it is as follows. The 
longitudinal membrane stresses which exist in the boundary layer, in 
combination with the longitudinal curvature, can be thought of as giving 
rise to radial forces. If then we consider a narrow lateral strip of the plate 
as a beam, it is subjected to normal forces distributed over the two end 
regions. The resultant of these forces is a pure bending moment about 4 
longitudinal axis which corresponds exactly to the lateral bending moment 
required to prevent the lateral curvature. 
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The purpose of the present paper is to attempt to generalize this concept 
{the boundary layer to apply to certain problems involving large deflexions 
f plates in which a similar phenomenon occurs but which are not as 
menable to a rigorous analysis as the simple case of a uniform strip under 
miform bending discussed above. 

\ particularly interesting case is that of a horizontal flat square plate 
subjected to four equal forces at its corners, two of them upwards at the 
ends of one diagonal, and the other two downwards at the ends of the other 
liagonal. According to the linear, small deflexion theory the plate 
ndergoes a uniform twist, so that any line parallel to an edge remains 
straight. It requires only the simplest of experiments to see that if the 
eflexions are appreciable the surface into which the plate deforms is quite 
lifferent from that given by linear theory. It will be found, in fact, that 
the surface appears to be perfectly cylindrical, the generators being parallel 

one or other of the two diagonals. The plate can be forced from one to 
the other of these two equally possible modes, although once it has started 

deform in one mode it will continue to do so under the action of the 
mer forces alone. (The reader can easily verify this for himself by means 

square piece of card of the type used for filing purposes.) Such a 
ylindrical mode of deformation does not satisfy the boundary conditions 
sociated with linear theory. It can be made to do so, however, by intro- 
lucing a boundary layer, of the type discussed above, along each edge of 
the plate. In this case the generators of the cylindrical surface are inclined 
tothe edge and in order to deal with the situation the von Karman equations 
for large deflexion are used, in conjunction with an order argument, to set 
ip boundary-layer equations which apply in this more general case. 


The problem of the torsion of a square plate described above is a par- 


ticular case of the torsion of a rhombic plate by forces applied at the corners. 


In this case the deformation mode is again cylindrical, although it will be 
ound that the mode which develops is always that which has generators 
parallel to the shorter diagonal. If the deformation is large enough, how- 
ver, it can be forced into the other mode, in which the generators are 
parallel to the longer diagonal. This second mode appears to be stable only 
the curvature is larger than a certain critical value; below this value it 
ecomes unstable and changes suddenly to the first mode. The solution 
ithis problem, which again involves a boundary layer, is given in the paper. 
order to deal with the stability problem involved, however, a considerably 
nore complete solution would be required. 

The other solutions which are given are concerned with the large deflexions 
‘cantilevered plates which are loaded so that they deform into a cylindrical 
de. 
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Since the work described in this paper was completed two interestiny 
facts have been brought to the authors’ notice. First, the transition } 
tween the small and large deflexion modes of a twisted square plate wa 
apparently first noticed by Kelvin and Tait (9), who describe it as 
‘remarkable case ... which deserves particular notice; not only as interesting 
in itself and important in practical application, but as curiously illustratiny 
one of the most difficult points in the general theory’. Secondly, the gener 
concept of thin plates deforming into developable surfaces has be 
used in two papers by Mansfield and Kleeman (10, 11) who are concern 
mainly with the deflexion of swept and triangular cantilever plates. T 
work described in the present paper is somewhat more fundamental tha 
that of Mansfield and Kleeman in that it explains how a developable surfay 
can be made to fit the boundary conditions of the problem, a difficult 


which was not considered by Mansfield and Kleeman. 


2. Uniform bending of a long strip 

Consider an initially flat long strip, of width 2b and uniform thickness! 
under uniform bending as shown in Fig. 1. If the longitudinal curvatur 
is 1/R, the simple theory of bending predicts that the cross-section assume: 
a lateral curvature equal to —/R. It is shown by Ashwell (6) that thisis 
a true picture only if a parameter ab, defined by 


eee 
vb = [3(1—p2)]! 
(Rt)! 
is small, i.e. a < 1. If this parameter is large, however, the cross-sectio! 


remains substantially flat except in a region adjacent to each edge of the 
strip, as shown in Fig. 2. The deviation ¢ of the middle line of the cross 
section from the axis through the centroid of the distorted section is give! 
approximately by 
put 1Vla . ~ . (” 
—- e-*Y(COS ay sin ay) é 
‘ 2(1 1°)! e 
oh Baeaeas f j 
in each of these two edge regions, or boundary layers, where y is the distance 
measured from the edge in question. 
If » = 0-32, the value of ¢ at each edge is given by ¢, 4 0-0975 
whilst the first two values of y at which f = 0 are equal to 0-613.( Rt) and 
3°066.( Rt). Thus the width? 6 of the boundary laver is of order (Rf) and 


Just as in the flow of a viscous fluid over the surface of a body the width 6 of t! 
boundary layer is not a clearly defined quantity. We could define it specifically by sty 
lating that if y 6, 2 ef where e€ is some chosen small quantity. However, this 1s not 
necessary, and for our purposes it is the order and not the absolute magnitude whi 


of importance. 





for large values, the parameter ab can be interpreted as expressing the ratio 
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cresting | ofthe width of the plate to the width of the boundary layer. The horizontal 
. ile to which Fig. 2 has been drawn corresponds to a value ab 10, which 
ves some idea of the value of ab required to produce such a boundary 
eres 
Fic. 2 
It should be noticed that the effect of increasing the longitudinal curva- 
ture, once the boundary layer has formed, is to narrow down the width of 
distal e layer without changing the peak values of ¢. Thus the maximum 
ition from flatness is only about 10 per cent. of the thickness, which 
0-097 practice is difficult to detect. Thus if an experiment of this nature is 
Rt erlormed, the surface appears to be truly cylindrical. 
Rt The membrane tension NV, in the longitudinal direction per unit width 
ne rat nthe y direction is given by 
” Et, > 
N C, (3) 
RS 
here £ is Young’s modulus. Thus the membrane tension is proportional 
the distance of the middle line of the cross-section from the axis through 
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its centroid. Furthermore, the membrane tension is confined to th 
boundary layer and the middle region of the strip is in a state of pur 
bending. 

It is clear that over the cylindrical middle region of the strip there mug 
exist a lateral bending moment J/, equal to —yD/R where D is the flexuy 
rigidity of the plate defined by D = E#®/12(1—p?). But M, is zero on ea¢ 
edge of the strip. The mechanism by which this change in the value of }/ 
in passing through the boundary layer is accomplished is as follows: Th 
membrane tension Nin the boundary layer, combined with the longitudin 
curvature 1/R, provides an effective force N/R per unit area, in a direction 
perpendicular to the surface of the strip. Thus if we consider a lateral ‘bear 
of unit width (as shown dotted in Fig. 1) the membrane tensions produc 
a lateral bending moment M, at a distance y from the edge given by 

y 
My, - * N.(Y)(y—Y) aY. 


« 


0 


For the purposes of this integration the cylindrical region between th 
two boundary layers corresponds to y = 0 so that, on using equation 
(2) and (3), the total lateral bending moment produced by the boundary 
layer is given by 


Mi, = R20) (y—Y )e-*¥ (cosa¥Y — sinalY) dY. 


Je 


0 
On integrating and using equation (1) this becomes 


_ED 


M ; 
, R 


which agrees with the required value. T'hus the important function of th 
boundary layer is to provide a lateral bending moment in the cylindrical 
middle region sufficient to inhibit lateral curvature. 


3. The boundary-layer equations 

Suppose that in a plate undergoing large deflexions a boundary layer 
of the type being considered has developed along a free edge of the plate 
as shown in Fig. 3. Let the direction of the generators of the cylindrical 
surface be inclined at an angle (17—@) to the edge. Choose two sets of 
coordinate axes Ox, Oy and OY, OX; the first pair are parallel and per- 
pendicular to the edge and the second pair parallel and perpendicular t0 
the generators respectively. 
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We denote by w* the deflexion of the cylindrical portion of the plate 


so that C2w* C*w* 
. QO 
Y? oXeY | rm 
; (5) 
o-wr l | 
ox? R 


where Ry, is the radius of curvature of the cylinder. In general Ry is a 
function of X only. 


Ax—DIRECTION OF 
/ GENERATORS 


fa y 
ae ae Seer | 


6 = BOUNDARY 
/ ™ LAYER WIDTH 








Fic. 3 
On transforming from the X, Y to the x, y coordinates by means of the 
p, xcos? y sin @ | 
) x sin 6+ y cos @ 
nd 
cos? 2w*  sin*é ow* sin @ cos @ 6) 
; —-. (6 
R 2 R X CX Yy Ry 
'y hypothesis the tensions and shear force NV,, V,, and N,., per unit width 
e zero in the cylindrical region. In the boundary layer they satisfy the 
itions of equilibrium 
\ \ ON,, c \ il 
U: : cy QO. (4) 
( cy Ca 
We now define / as a lengt haracteristic of the length of the edge of the 
( e width of the boundary layer. Also let 
x Y 
(8) 


re) 
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so that in the boundary layer both € and y are O(1). Equations (7) nog 


become ' yo é ee 
oN ry _ fa) ON, C N, Dy OC - 
on l < 's C n , 3 € 


It will be seen therefore that 


9 


Po 0(°), N, o(5). : 
N, iL | N, [? 

Again by hypothesis the width 6 of the boundary layer is small compared 
with the typical dimension of the edge /, so that in the boundary layer ¥: 
compared with \V.. Thus the membrane forces ir 


r 


can neglect N, and N,, 
the boundary layer consist primarily of tensions in the direction of the edg 
Now let the deflexion of the plate in the boundary layer be given by 

w wrt Cl, (10 


where ¢ represents the departure from the cylindrical form. The equatio 


of equilibrium (12) of the plate in a direction perpendicular to the sy- 


plane is 2 ne 2 
C-w O-w o-w 


DV4w q+ N, — | 2N N. (1] 


02x? "Voxoy | ” dy?’ 
where q is the applied pressure on the surface, V‘ is the biharmonic operator 
and D is the flexural rigidity of the plate defined as 
Ee 
12(1—p2)" 
In the cylindrical portion, where the membrane stresses are zero, equatio1 
(11) reduces to DV4w* = q. (19 
Hence, on using equations (5), we see that q is a function of X only, whicl 


means that a cylindrical surface is possible only if the surface pressure is 


constant along each generator. 
It will be seen from equations (6) that in general the values of @w éz* 


6?w/oy", and é2w/éxéy are of the same order of magnitude. Hence, on 
substituting equations (10) and (12) into equation (11) and neglecting the 


terms involving N,,, and N, we obtain the equation 


ao 
, O-w 
pvt = n,2”. [ 
” Gaon 
From equation (8) we have 

» 1a, 2 al lag 

V*C a al xan eee ; 

0* Cn 0°l? 0€70n? (U4 ag4 
and since 6/1 < 1, we need include only the first term of this expression, 


so that approximately 
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Hence on using the first of equations (6), equation (13) becomes 
ot cos"6é ., 
AN. (14) 
cy" DR, 
The second of the large deflexion equations of von Karman (12) can be 
written in the form 
> y . c Th . C 2 4 2p Ps 
v2(N.+A Et} | etre! | (15) 
CxXCY Ox" oy* 
Now V, can be neglected compared with V.. Furthermore, we have 
1 oN 1 o2N 
o- on /2 C&* 


nd only the first term of this expression need be retained. Thus the left- 
ind side of equation (15) is approximately equal to c?.N,,/éy?. 
Using equations (6) and (10) the term in square brackets on the right- 


hand side of equation (15) is equal to 





[2Psnbeosé atl sin?6 670 cos*@ 6 ( real ) eC aC 
R, Cxcdy Ry cx?” Ry Oy? excy Cx oy" 








,4 o'sS org o(7) (16) 





OXCY oy” l 
»y 9¢ (9 
oC Ca & o* ~" 
wna a = ol 3} (17) 
cx} Cy” [? 
so that the first two terms in the above expression may be neglected 
comparison with the third. Hence equation (15) now reduces to 
l &N cos*@ oC OC \? al aC , 
ee 5+ Hee — a5 aa (18) 
Et ey" R x cy” CXCY ox" oy* 
In the case of the uniform strip under pure bending discussed earlier we 
found that O(t) and that the width of the boundary layer 6 = O(. Rx t). 


We shall now assume that these orders still apply and confirm the assump- 

tion later. The ratio of the term in square brackets to the first term on the 

right-hand side of equation (18) is then of order 5?//? and hence only the first 
term need be considered. Equation (18) then reduces to 

2N Et cos?6 e27 

= (19) 

y” Ry dy? 

Equations (14) and (19) are the boundary-layer equations which take the 

place of von Karman’s large deflexion equations. If we multiply equations 
14) and (19) we obtain the equation 

O4C o2N 12(1—,?)cos*4 vy! C 


C y* oy* t? R*. "s y2- 
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Hence on using the second of equations (8) we see immediately that 

8 = O(N(R x t)). (20 
Furthermore, it will be seen from the discussion of the uniform str 

under pure bending that the membrane stresses in the boundary layer ax 

responsible for developing a bending moment M, in the cylindrical portion 

of the plate which is of order D/R,. Hence 


iD 
so that N, (3). 
‘ o*, 
; ‘ Et? 
that is N, o| =} (21 
\Ry, 


Finally, from equation (17) we have 


; IN, Rx 
ol- _ ioat_9 Yi 
. Et 


and on using equation (21) this becomes 


t = Olt). (22 


Equations (20) and (22) confirm the assumption which was made in deriving 
equation (19). 

We can integrate equations (19) by observing that if the boundary layer 
is thin we may assume that Rx is independent of y for the purposes of this 
integration. Thus 

Et cos*é 


N, 
Ry 


C+ y f(x) +fo(2), 
where f, and f, are arbitrary functions of integration. But as y increases 
indefinitely both N, and ¢ approach zero so that f, = f. = 0. Hence we 
have . P 
. Et cos?@ 
NX, af 
Ry 


és (23 


so that the membrane stress in the boundary layer is proportional to ¢ 
the deviation of the middle surface from the cylindrical form. 

On substituting equation (23) into equation (14) we now obtain the 
equation 


“y 
C 
> 447 0 (24 
oy’ 
. 9 
’ i ein One on 
where a? = [3(1—,p?)}} (25) 








Aga 
equati 


‘ 
But as 


to det 
edge. 


The fi 
M, is 
condit 

On 


(Consi 
(22) t 
order 
is of 


dition 


On s} 


and | 


If th 
edge 
the | 
ence 
surf; 

bent 











LARGE DEFLEXION 





OF THIN PLATES 201 


Again, assuming that in the boundary layer Ry is independent of y 
equation (24) can be integrated to give 

C e~*Y| A(x)sin ay+- B(x)cos ay|+-e%[ P(x)sin ayt Q(x)cosay]. (26) 

But as y increases indefinitely ¢ tends to zero so that P = @ = 0. In order 

to determine A and B we must now consider the conditions on the free 


edge. These conditions are that 


o~2 o-u - 
0 aty G. (27) 
cyu- Cars F 
oi Cu 
2—p)— 0 aty 0. (28) 
oy? Ox=Cy ; 


The first of these two equations expresses the fact that the bending moment 
VM, is zero on the edge; the second equation is the Kirchhoff boundary 
mdition which says that the equivalent shear force is zero on the edge. 


On using equations (6), (10), and (17) these two conditions become 


‘. ] 
(sin“0+- 4 cos“0) — at y U, 
¢ um ly 
; ann @i 4 
sin 6} sin*t 2— )cos*6@] ‘| at y V0. 
dX\R, : 


Consider the second of these two equations. On using equations (20) and 
22) the left-hand side is of order 1/(R, 6) whilst the right-hand side is of 
rder 1/(R Hence we may replace the right-hand side by zero since it 

fa smaller order of magnitude than the left-hand side. The two con- 


tions for determining A and B now become 





(Sint Lcos-0) Feb 
; = Oo, aty V0. (29) 
R oy” ; 
| On substituting equation (26) we finally obtain 
: t(u+-tan*é) 
fo yt? 
| J 120 pL ) 
| ( | tO? I nes 
tan°@) , 2 
: — e—*4(COs ay — SIN ay). (30) 
| \ A és | [L™ } i 
| If the generators of th cylindrical surface are perpendicular to the free 
| edge (/ = 0) equation (30) becomes identical with equation (2) which defines 
} 1 . ° . ° ° ° 
| boundary layer of a uniformly bent strip. There is one essential differ- 
} €nce, however; even if 6 0 the radius of curvature Ry of the cylindrical 
Suriace, and therefore a, is in general a function of 2, whereas in the uniformly 


bent strip both Ry and « are constant. 
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4. Bending of cantilevered plates 

Consider a thin plate ABCD of the shape shown in Fig. 4, fixed alon 
the edge AD and perfectly free otherwise. Suppose that it is subjected 
a normal pressure on its surface which is a function of X only, so that alon 
any line such as PQ the pressure is constant. In addition it may be subjected 
to a shear force and a bending moment 1, uniformly distributed over th 


GENERATORS OF 
’ GYLINDRICAL SURFACE 
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tip section BC. Then, if the deflexion is large the plate will bend into a 
cylindrical surface, whose generators are parallel to the root AD, except 
for a boundary layer along each of the two edges AB and CD. The curva- 
ture of the cylinder can easily be calculated. For, if we neglect the contribu- 
tion of the membrane stresses in the boundary layer to the total bending 
moment at a section such as PQ, we have 


| M 


> (3] 
Ry 2bD 


where M is the total bending moment at the section PQ due to the applied 
loads and 20 is the width of the section. 

The deviation of the plate from this cylindrical surface in the boundary 
layers is given by equation (30) where « is defined by equation (25) and Ry 
by equation (31). The membrane stress exists only in the boundary layers 
and is given by equation (23). 
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Clearly this is the solution only if the width of the boundary layer at 
y section is much less than the width of the plate, i.e. if \(tRy) <b. It 
1y be that the loading is such that it is impossible to satisfy this condition 


ill sections; for example in the case of uniform tip shear force the 


ending moment is zero at the tip, and however large the applied shear 
ree, there will always be a region near the tip in which the condition for 
boundary layer to form is not satisfied. The plate will then assume a 
viindrical surface only inside a line such as that shown dotted in 
Fig, 4. 

Furthermore, there can be no boundary layer at the root section AD 


since the clamp will ensure that there is no distortion of that section. Hence 
, small region near the root section must also be excluded. 
The special case of a rectangular cantilevered plate is given by @ 0, 
| that of a triangular cantilevered plate is obtained by making B 
id C coincide. These are both included in the more general solution of 
the plate shown in Fig. 4 
It will be seen from equation (30) that the value of ¢ at each free edge 
the plate (y 0) is given by 


t(u+-tan6@) (32) 
Ja 


9 


V1 2(1—p*); 


[his is independent of Ry so that once the boundary layer has formed ¢, 


remains constant, both along the edge and with further increase of load. 


It will also be seen that ¢, increases numerically as @ increases. If up = 0-32, 
for example, the numerical value of ¢, when @ 45° is over four times as 
great as when 6 0. Even then, however, it is only about 40 per cent. of 


the plate thickness so that the maximum deviation from the cylindrical 


form is not large 


5. Torsion of a rhombic plate 

Suppose that a plate in the shape of a rhombus (Fig. 5) is loaded by four 
equal forces P at the corners, in a direction perpendicular to the plane of 
the plate. The forces at A and C are applied in the opposite direction to 


those at Band D. On the assumption that the deflexions are small enough 


rt the linear small deflexion theory to apply, this problem has been solved 
by E. Reissner (13). The expression for the deflexion w of the plate obtained 
irom Reissner’s work is 
P [ v2 — y? (x2 + y2)eos 20 
u 0 Yo rl 0 Yo) (33) 


${)sin 26} | pb l+p 
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where the x, and y coordinates are as shown in Fig. 5. On differentiatiny 
this expression to obtain the curvatures of the surface we find that 


C2w P ] 
l—p 


cos | 


ex 2D sin 20 ltp 
Cw P ] cos 26 
oyr 2Dsin26}1—p =1+p] 


Thus the curvature is constant in both the x, and y, directions. Further. 


GENERATORS OF” 
CYLINDRICAL SURFACE 


more, the curvature in the x, direction is positive and that in the y, direction 
negative. Thus the plate bends into an anticlastic surface with constant 
principal curvatures. 

Such a surface, however, is not developable and therefore the plate 
stretches, and membrane stresses are introduced. The presence of these 
membrane stresses modifies the deflexion surface and by the time the 
deflexions become large the shape bears no resemblance to that predicted 
by the linear theory. In fact a simple experiment showed that the plate 
bends into a cylindrical shape, the generators being parallel to the shorter 
diagonal of the plate. Such a deflexion mode does not satisfy the boundary 
conditions, however, but this can be explained by the presence of a boundary 
layer of the type under discussion which forms along each edge of the plate. 
This was easily visible in the experiment which was performed; also it was 
noticeable that near the corners, especially the acute ones, where the 
boundary layers interact there was an appreciable departure from thi 
cylindrical shape. On the whole, however, the plate deflexions are ver’ 
well represented by the cylindrical shape. 
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[It is easy to calculate the curvature of the cylinder. If we again neglect 


the contribution of membrane stresses in the boundary layers to the total 


bending moment carried by a section parallel to BD at a distance X from 


where X bcos @) we have 
' 2X tan@ D 
PX aa 
Rx 
| , 
¢ cot é. (54) 
R, 2D 


Thus the curvature of the cylinder is constant, so that the width of the 

undary layer is also constant. The deviation from the cylindrical form 

the boundary layer is given by equation (30) and the membrane stresses 
by equation (23), where a is defined by equation (25) and Ry by equa- 
ion (04) 

Using either equation (33) or (34) it is a simple matter to calculate the 
elative deflexion of the corners A and C from the corners B and D from 
either the small deflexion theory or from the large deflexion theory given 


ere. If these are denoted by A, and A, respectively we have 


Pb? ; 
J —— sin 26(1-+--+- 2 cot*2¢é), 
1(1—p*)D 
Pb? 
\ cot 6 cos*6é. 
1p) 
” A 2 tan*6 . 
hus —(1+-y+ 2 cot*26). 
\ l—p? 
In the case of a square plate (@ = 45°) with p 0-32 this gives A, = 2-94A,, 
whilst if @ 0°, A |-47A,. Hence the small deflexion theory gives a 


considerable overestimate of the deflexion of the plate. 

It is interesting to observe that a square plate can bend with the genera- 
tors of the cylindrical surface parallel to either diagonal. Once it has started 
to deform in one mode it will continue to do so unless forced into the other 


mode. In the case of the rhombus, however, the plate always starts to 


deform in the mode with the generators parallel to the shorter diagonal. 
If the curvature is not very large and an attempt is made to force it into 
the other mode, with the generators parallel to the longer diagonal, it springs 


back into the first mode as soon as it is released. If the curvature is suffi- 
ciently large, however, the second mode appears to be stable though the 


plate will never deform into it without being forced. 
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6. Comparison with other well-known boundary layers 

The novel character of the boundary layer discussed above can be broug} 
out by a comparison with other well-known boundary-layer phenomen 
such as Prandtl’s boundary layer in a flow at large Reynolds number, 
Kelvin, Tait, and Friedrich’s boundary layer in the bending of plates wit 
regard to Kirchhoff’s boundary conditions. In these classical examples 
the effect of boundary layer on the field outside is such that the terns 
involving derivatives of the highest order in the complete differenti: 
equations are negligible. In contrast to this, the effect of the membran 
stress boundary layer discussed above is to make the terms involving th 
highest power of the dependent variables negligible in the field outside th 
boundary layer. The equations of motion for the flow outside Prandtl’s 
boundary layer remain non-linear, while the equations for the bending of 
the plate outside the membrane stress boundary layer become linearized 
The highest derivatives in Prandtl’s case represent the effect of viscosit 
of the fluid: when they are neglected the fluid behaves like a non-viscous 
one. The highest derivatives in the plate-bending equations represent th 
effect of bending rigidity: they are fully effective outside the membrane 
stress boundary layer. 

It may be pointed out that although the derivation of the boundary 
layer presented in § 3 fulfils the purpose of the present paper, the basi 
assumptions can be made less restrictive to cope with more general loading 
conditions. For example, the assumption that the deflexion surface outsid 
the boundary layer is cylindrical can be replaced by simply being develop- 


able: the tractions on the unsupported edges may be finite. 


7. Some general observations regarding ‘applicable surfaces’ in 

the theory of plates and shells 

From equation (15) it can be seen that if membrane stresses vanish in a 
plate the deflexion surface is a developable surface. This is related to a 
number of important phenomena concerning bending and buckling of thin 
shells. In differential geometry, two surfaces are said to be ‘applicable’ to 
each other if one can be deformed continuously into the other without 
stretching or tearing. A developable surface is applicable to a plane. Hence 
the phenomenon discussed in this paper can be stated as that under certail 
conditions the plate deflects in such a manner that, except for a boundary 
layer, it remains applicable to the original surface. 

It is well known that when a cylindrical shell is subjected to a sufficient] 
large end compression it buckles into a diamond-shaped deflexion patter 
(see sketch in Fig. 6). An ideal diamond pattern consists of a series 0! 


triangular flat surfaces and is applicable to the original cylinder. The actua! 
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ickling pattern deviates slightly from the ideal applicable surface for the 
lowing reason. The bending rigidity of a shell varies as Et?, while the 


lity resisting the stretching of the mid-surface of the shell varies as Et. 


Hence as tne thickness t becomes very 


ull, it becomes far easier to bend the 







ell than to stretch it Hence it is 
ti] t expect that the buckling 
le of the shell should involve little 
ching. However, if the buckled sur 
is exactly applicable to the original Me 


\ 


' 
\ 


ce, the stretching of the mid-surface —}\ 


\ 
ishes identically and there is no strain tS 
f% 


nergy due to membrane stresses. The 


g | ? 

nding strain-energy nowever, be ! | 
mes large in spite of the small vending } | 
gidity because of the high curvatures | | 
r the idges. The actual physical | | 
henomenon strikes a balance to require | : 
small stretching strain energy to re- L y, 
ice the bending strain energy. The cl 
esult is a mode that is almost applicable Fic. 6 


the cylinder. Although it would be 

ong to assume that the deflexion mode is exactly applicable to the 

inder, the recognition of the diamond pattern enables an approximate 
tation of the buckling load of cylinders to be obtained, as demonstrated 
m Karman and Tsien (14) and Leggett and Jones (15), ete. That 


e diamond pattern is almost applicable to the cylinder has been discussed 


recently by Yoshimaru (16 


Similar situations arise in a number of plate-bending problems. A case 


l interest to aeronautical engineers is that of the deformation of a swept 


ving und id. An idealization of this would be the case of a thin flat 


rip clamped obliquely at the root cross-section and subjected to some 


of loading at the tip. A simple experiment with such a configuration 


indicate that the plate deformed into a surface which contained 
light generators. These generators were perpendicular to the edges of 
strip at some distance from the root, but in the root region they rotated 
| right at the support the generator coincided with the root section. 
e question therefore arises as to the possibility of obtaining a solution 


1 


rthis pi em in which the mode of deformation is a developable surface 


wit 1 bor lary layel along each edge in order to deal with the edge 


If the loads are applied at the tip the deflexion function would 
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also have to be biharmonic in order to satisfy equilibrium requiremen 
However, it is shown in the appendix that it is impossible in general ; 
satisfy the biharmonic equation and the developable condition simuly, 


neously. The only surfaces which do so are cubic cylinders (and possib)) 


SSE CRETRARERIARY 
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conical surfaces although this has not been fully investigated). If the plate 
is loaded by means of a pressure on its surface the right-hand side of the 
biharmonic equation is no longer zero. However, in general, the distribu- 
tion of pressure required to obtain a solution in the form of a developable 
surface is so unrealistic that further investigation is unwarranted. Thus 
the deflexion surface of a swept plate will not be an exact developable 
surface with a boundary layer along each edge. However, when the plate 
is very thin the deviation of the deflexion surface from a developable surface 
is small. 
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In Fig. 7 the lines of zero curvature for a swept plate under a load at the 
re sketched. Similar sketches are given for the deflexion pattern of a 
ry thin rectangular plate which is twisted beyond the range of validity 
he linear torsion theory, and that of the torsion of a rectangular plate 
ted to forces not acting at the corners. These patterns can be demon- 
ted by very simple experiments. These deflexion patterns would be 


ful tor deriving approxin ite solutions by energy methods. 


APPENDIX 


Surfaces that are both developable and biharmonic 


ypal ur r than a cylinder ora cone there is associated a 
This LS} curve whose tangents sweep out the developable 
17 18 
t r(s) be t position vect { a point on the line of regression, where s is a 
ng the posit f the point on the line, i.e 
r f(s), g(s), A(s)). 
cto! Tay { n the de velopable surface is 
ris) ur(s). 
3) + u f(s) 
ats ud(s) (35) 
h(s)—uh(s) 
| the distance of the point along the tangent from the point of tangency 
ne of regression, and lot refers to differentiation with respect to s. Hence 
u f* 
ws " (36) 
J u 
1—f9)—h(fa—gf)+hfa—fa 
rrespona between the vr,y and uw,s coordinate systems, the 
! e above ex nm must not vanish. Writing 


q* als S$) 
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If the developable surface given by equation (35) is biharmonic, i.e. if V4w 


the coefficients of u~*, u-4, and u-> must vanish independently. From the last ten 


in equation (37) it is seen that the only possibility for this to be true is dé = 0, B 
¢ = 0 implies % = 0, which in turn implies é?w/éx? = é*w/éy? = 0, and the surf 
must be a plane. 

Developable surfaces that are not represented by equation (35) are cylinders ay 
cones. It is easy to see that any cubic cylinder is both developable and biharmo, 
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SUMMARY 


[It is shown that the boundary conditions and differential equations governing the 


brium of an anisotropic homogeneous elastic solid, whose 21 elastic constants 
tisfy 14 specially chosen conditions, can be expressed by means of a linear trans- 
mation in a form identical with the corresponding equations for an isotropic solid. 
s gives 1 to a new method of treating equilibrium problems in the classical 
y of elasticity for solids having anisotropy of a special kind. In the case of a 

| which has three mutually orthogonal planes of elastic symmetry, 9 of the 14 
nditions are automatically satisfied, and the transformation involves simple 
inges of scale in the symmetry directions; in the particular case of transverse 
ropy, the number of conditions automatically satisfied is 12. The transformation 
applied to problems involving two or more elastic solids if they have the same 

stie constants and are similarly oriented with regard to their directions of elastic 
try. As an illustration, the transformation is used to derive, from Hertz’s 

tion for the isotropic case, expressions giving the shape and size of the contact 


gion formed when an elastic solid having a special kind of elastic anisotropy is 


ssed light gainst a rigid plane. 
In the case of a transversely isotropic solid, it is shown (without restriction on the 
mstants) that the class of solutions of the differential equations of equilibrium 
able by a method due to Elliott may be enlarged by using three stress functions 
mit pe instead of two. This leads to an interesting alternative derivation 
ms on tl istic constants for the case of transverse isotropy. 


1. Introduction 
Ix considering the solution of a problem involving the contact of trans- 
ersely isotropic elastic solids it was found that a considerable simplification 
esulted when the five elastic constants involved were taken to satisfy two 
nditions. When these conditions were introduced in the general equili- 
rum equations of elasticity theory, it was found that these simplified also, 
nd were expressible in terms of new variables in a form identical with 
it of the corresponding equations for isotropic solids. Any given problem 
nvolving homogeneous solids with the right kind of anisotropy could 
ereby be transformed into an isotropic problem, in which form it should 
less difficult to solve. A brief note of this result has been published 
lsewhere (1 lhe object of the present paper is to extend these trans 
mations to include as wide as possible a class of anisotropic solids. 


Quart. Journ. Mech. and Applied Math., Vol. VIII, Pt. 2 (1955)] 
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It is found that even with the most general linear transformation (th; 
variables being referred to a rectangular Cartesian coordinate system), th; 
21 elastic constants for the general solid must satisfy 14 conditions, so thy 
only 7 constants are independent. For solids having orthorhombic elast 
symmetry all but 5 of the conditions are automatically satisfied, and ther 
are 4 independent constants; for transversely isotropic solids all but 2 ¢j 
the conditions are automatically satisfied, and there are 3 independen 
constants, as mentioned above. Since non-linear transformations are w 
likely to be of interest, the extension given here probably represents t! 
best that can be achieved by the present method. 

The following assumptions and restrictions are made throughout: 


(i) The solids are ideally elastic, and the strains are infinitesimally sma 
(as in the classical theory of elasticity). 

(ii) The elastic constants satisfy equations (11) below [or (16) for orth 
rhombic symmetry, (26) for transverse isotropy ]. 

(iii) Only equilibrium problems can be dealt with; motion and vibratio: 
are excluded. 

(iv) The elastic solids are homogeneous. When two or more such solid: 
are involved in one problem, they must have the same elastic co 
stants and the same directions of elastic symmetry. The hom 
geneity restriction can be relaxed with regard to two of the constants 
but this is not likely to be of interest in practice. 


The method is based on the use of new dependent and independent 
variables. The ones used in the present paper involve the choice of a speci 
rectangular Cartesian coordinate system (for the case of an orthorhombi 
solid, the coordinate planes must be parallel to the planes of elastic syn- 
metry), to which the equations of any given problem must be referred 
before the transformations given in the paper can be applied. The use of 
a special coordinate system is not, however, essential to the method; the 
actual transformations to be applied in any other coordinate system could, 
if required, be derived from those given in this paper. 

In section 2 the basic transformations of the differential equations 
equilibrium are derived, together with the restrictions which have to be 
imposed on the elastic constants. In section 3 the special case of solids 
with orthorhombic elastic symmetry is treated, a summary of the results 
being given in section 4 in a form convenient for application to actual 
problems. In section 5 the special case of transversely isotropic solids Is 
treated, and a digression is made in order to give a generalization of Elliott's 
transformation. In section 6 the transformation of boundary conditions '!s 
considered briefly, and in section 7 the methods are applied to obtain @ 
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eneralization of Hertz’s solution to a problem involving the contact of 
elastic solids. The paper ends (section 8) with a discussion of the applica- 


bilitv of the methods to actual materials. 


2. Transformation of the differential equations 
In the classical theory of elasticity the partial differential equations which 
ive to be solved to determine the equilibrium state of an anisotropic elastic 


slid subjected to given forces and constraints may be written in the form 


CY * 
Pir 4 pxX ; 0, (1) 
CX, 
i ov i WW (2) 
=} j : “ 
C€;; Ce;; 
ou ou ‘ 
% Fan J (3) 
J ° y 
c uj CX; 
J . 9 . > 4 ! 
NW C11 O71 FC y2 C11 Cog HC1g C11 € 3g H 241 (C14 C23 + C15 31 1 C16 © 21) 
1 PA ) » 
dC a0 C39 + Cog Cno €gg + 2lg0(Co4 Cog + Cos C31 + C26 C21) + 
1 2 9 ’ , , | 
3 C33 C33 + 2l33(Caq Cog +Ca5 €g1 + Ca¢ a1) 
9 t. 9 a. 
2 93(C4q Cog + =C45 €g1 + 2C gg Co.) 4 
) 2 o%% 2 
2655 €31 + 456 €31 C12 + 2Ceg Cie: (4) 
where p;;, €;;, %;, U;, X,;, p, W denote the stress components, infinitesimal 


strain components, coordinates, displacement components, components of 
body force per unit mass, density, and strain-energy per unit mass. The 
coordinate system is rectangular Cartesian fixed in space. Latin suffixes 
take on the values 1, 2, and 3, and each term containing a repeated Latin 
suffix is to be understood, here and throughout, as summed for the three 
values of that suffix, unless the contrary is stated. 

The coefficients c.g (a, B ae 6), which are constant for a homo- 
geneous material, characterize the elastic properties of the material; for 
the most general type of anisotropy, 21 of them are independent (cf. (2), 
p. 159, allowing for differences in notation, namely e,, = €1;,€yz = 2€gg, ete.; 
our ¢,g are the same as Love's). 

The differential equations are linear; in certain circumstances there exists 
1 linear change of variables of the form 

li. >2Z:, Uz:>U,, ¢ > Cis, Pig > Pijs Xi > Xj, (5) 


where 
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the a,; being constants satisfying the condition 
det(a;;) - 0 


for linear independence of the x 


;, Which transforms the differential equ 
tions into a form identical with that of the corresponding equations for a 
isotropic solid: equations of the same form as equations (1), (2), and (3) ay 
obtained, while the strain-energy function (4) takes the form 


Q , , 


W= (308; 55,54 19 j, O;5)€;; Eps 


appropriate for an isotropic solid of Lamé elastic constants A, pu. 5;; equal 
1 or 0 according as 7, j are equal or unequal. 

In order to verify this statement, perhaps the simplest procedure is { 
start with the isotropic equations and work backwards. We need onl 
outline the proof, which is straightforward, though lengthy. The thre 
equations (1), (2), and (3) are found to be invariant in form under thy 
transformation (5), whatever the values of the transformation coefficients 
a;;may be. Application of the inverse of transformation (5) to the isotropi 
form of the strain-energy function (8) leads to an expression which may bx 
written in the form 


é 


13 \ 


Y = (3AA,;; As I vA; A;,) 


Ci; 


where A 5; = Ay; Ay; (10 


This expression, which involves only 7 independent constants (6 indepen- 
dent A;; and the ratio A/u, for example), has to be identified with the 
strain-energy function (4) for the general anisotropic material. This is 
possible when (and only when) the 21 constants c.g in the latter satisfj 


§ 


certain conditions, 14 in number, which are found to be expressible : 


follows: 


on | a . . } 2 la 
2C441+-Cog = (Coo C33)? + 4C34/Coo, 
op . » } n2 lp 
“C551 C3, (C33 €y,)*? +4035 /C33, 
ee eee 2 
2C 6+ Cy2 = (C41 Cog)? + 476 /C11, 


In ‘Se me oe — 
256 t-Cyg = Coq(Cq1/Cog)? + 2Cg5 Cy (C33 C41) 


on » - > $ | 9 ~ + 4 
<C64 1 Co5 €35(Ca9/C33)? + 2Cy6 Cog(Cy1 Coo) *, 
rn » > } » » > 4 
2C45 1 Cgq = Cy6(Cg3/C11)° + 2Cog Cy5(Coo C33) *, 
a p la 4 
€15/¢35 (€31/€33)° 
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2 C33)" — 4C24/ Coo (C33 C11) 4035 C33 (Ci Co9)}- -4ci¢ Cn 
” 4C54/Coo C31 — 4C35/C3g Cyg—4C56/C) 
Co4(C11/Co2)° — 2Cg5 Cy6(Ca3 C11) : °33)* - -2C16 Coa(C11 Coe) 
Cyq— C35 Cig (Cg C11) 246 C24(C11Co2) 
C16 (C33/C11)° — 2Coq Cg5 (Con Cag)” (11) 
€3g— 2C 94 C35(Co2 Cgg)~* 
1+ 2u/A (12) 


Equations (11) are therefore a necessary and sufficient set of conditions 


for the transformation to an isotropic form by means of a transformation 


fthe form (5), (6) to be possible. This form can be shown to be the most 


general linear one which leaves the differential equations (1), (2), and 


3) invariant. The positive signs have to be taken for the radicals in 
11) (and throughout the paper) in order to give a real transformation (6). 
In order to make use of the transformation in an actual problem, we 


have to be able to express the transformation coefficients a,;; in terms of 


the given elastic constants c,g. The relations required for this may be 


btained by identifying the two forms (4), (9) for the strain-energy func- 


tion. and are found to be expressible as follows: 


Ay [¢4,/(A- 2) |* 


A 56 Coo (A T 2) |! 
A 33 C33 (A+ 2u)|! 


(13) 

Ay, 2CyqlCoo(A 2) | 

Az, 2¢35| C33(A+ 2) | 

As 2€ 461 C4,(A-+ 2) | 

The elastic constants A, », of what may be called the ‘equivalent isotropic 
iterial’, are determined only in so far as their ratio is given by equation 
\2), and they may be given any values consistent with this equation. 
There is also a certain freedom of choice in the values to be given to the 
} transformation coefficients a;;, because they are determined only by 
juations (13), where they occur in only 6 independent combinations 4 ,, 
The transformation as given above has been restricted to materials which 
homogeneous (i.e. the c.g are constants). This restriction is likely to 


x 


most felt in those problems which involve more than one elastic body, 


fan 3 


rit means that the bodies must not only have the same elastic constants 
ut must also be similarly oriented with regard to their directions of elastic 
symmetry. A slight, but probably unimportant, relaxation of this homo- 


zeneity requirement can be made, because the c,g can be regarded as 


‘unctions of a;;, A, and yw, and it is the transformation coefficients a;; which 
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have to be constant, not A and yz, which can be functions of 2,, ice, thy 
equivalent isotropic material need not be homogeneous. 

It has been stated that the transformation to isotropic form is possibj 
only for equilibrium problems, but in view of the importance of elastic way; 
propagation in anisotropic materials, it is worth noting what happens unde 
the above transformation when the inertial term for small motion 
p0*u,/et?, is added to the right-hand side of equation (1). The transforme 
equation is found to be 
Dir 1 pX’ Prt: 


, — 
OX, ct 


A,,. 4 
When the inertial term is not zero, it is seen that the form of this equatia 
is invariant only when A,; = 6,;;, which corresponds to the degenerate cas 
in which the original material is isotropic and the transformation can lx 
interpreted simply as a change from one rectangular Cartesian coordinat 
system to another. 

Some (but not all) of the 14 conditions (11) on the elastic constants ¢,, 
are automatically satisfied when the material possesses lines or planes oj 
elastic symmetry. Two such cases will now be considered. 


3. Solids with orthorhombic elastic symmetry 
For a material with orthorhombic symmetry or with three mutual 
orthogonal planes of symmetry, 14 of the 21 elastic constants c.g are zer 
and the strain-energy function may be written in the form 
y l(e 2 . 2 2 . 7 7 i, 
HW 2(Cy1 €71 + Co9 €29+ Cgg €33) + Cog Cae €gg + Cay €gg C11 + C12 C11 Cn2t 


22 »9 


+ 2(C4q C33 +55 €31+ Coe C2), (15 
the coordinate planes being taken parallel to the symmetry planes (2, p. 154 

In this case it is found that 9 of our 14 special conditions (11) are satisfied 
leaving 5, which may be expressed in the form 


2C 44+ C23 (Cop C33) 
2¢55+C31 (€33 C44)’ (16 
266+ C12 (C11 Co9)* 

(Cop Cgg)?/Cog (C33 C11)" /Cgy (C11 C2)" /Cy2 | 1+ 2y/A] 


in terms of the c,g, or, alternatively, in the form 


893+ 3844 = (So9 S33) 
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nterms of the elastic moduli s,2, which are related to the c,g by equations 


f the type 
A = Cy9C33—Co3; 893 A = Cy Cyg— Cy Ca; Syq = Cqq’; ete., (18) 
vhere A = det(e;;)  (t,) 1, 2, 3), 


nd to the Young’s moduli #;, shear moduli G;, and Poisson’s ratios o;;, 


by equations of the type 
— _ -_—* oe ( 
BE, s,', G, a;,*, Org * 812/81); ete. (19) 


ind o,. denote Young’s modulus and a Poisson’s ratio for a tensile stress 


n the direction Ox,, and G, denotes the rigidity modulus for a shear in 


l 


which planes normal to Ox, move in the direction Ox,; the other symbols 


have similar definitions and o = 3A/(A+p) is the Poisson’s ratio of the 
equivalent isotropic material 

Equations (16) or (17) are thus the conditions to be satisfied in the case 
f an orthorhombic solid if the transformation to isotropic form is to be 
possible. These conditions have some direct physical significance, for the 
first three of equations (17) are the conditions for the polar diagram of the 
reciprocal of the fourth root of Young’s modulus to be ellipsoidal (instead 
fa more complicated quartic surface (2, p. 161)), while the two remaining 


equations can be expressed in terms of Poisson’s ratios in the form 


, . - 2 9 
23 %32 531 [13 912 9%) S |. (20) 


Another form of the conditions which is perhaps worth noting is the 
following ya . 7D 
te G,/Gs (Ey E;) 
¥21/ F293 G3/G, (E, E,)} ' (21) 
. , » Gy’ \t 
%32/%31 oe (Ey E,)* 


These are necessary but not sufficient, only four of them being independent. 
In the case of orthorhombic symmetry the transformation coefficients a;; 


can be given the simple set of values 
a,—0 (i +), | 
a E./E)t (i not summed), ! 


where £ is an arbitrary positive constant, equal to the Young’s modulus 
of the equivalent isotropic material. This follows from equations (10), (13), 
16), (18), (19), and the known relation E = p(3A+-2u)/(A+p). 

In the next section we collect together all the results, for the case of 
orthorhombic symmetry, which will be needed in applying the present 


transformation methods to any actual problem. 
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4. Summary of transformations for the case of orthorhombic 
symmetry 
In the differential equations of equilibrium, referred to a rectanguly 
Cartesian coordinate system whose coordinate planes are parallel to the 
planes of elastic symmetry, let the variables 


x; Uu; é 


4 —_ . 99 
ih a? ij? Piz, x i —_ :. a 3), (23 
representing respectively the coordinates, displacement components, in- 
finitesimal strain components, stress components, and components of body 


force per unit mass, be replaced by the quantities 

E\it, /(E\i, E? \}, E,E\} , E\iw |. ; 

ji 4 -)* a’., =1"€, oi" Bis, ')*_X* (7, ) not summed 
EB i E i EB E ij D) ij ’ i P 

) v; i; E; 


k* E 
where £; are the Young’s moduli in the symmetry directions, and £ is a 
arbitrary positive constant. If the elastic constants of the solid satisfy 
equations (16) (or (17)), the resulting equations in the new variables z 
u;, ete., will have the same form as the corresponding equations for ar 
isotropic elastic solid expressed in terms of the usual variables, x,, u 


(24 


+ Ct 
The Young’s modulus of this ‘equivalent isotropic solid’ is equal to the 
arbitrary constant H; the Poisson’s ratio o is determined by the Poisson's 


ratios o,;; of the given anisotropic solid according to the equations 


o = O52 Cao Oa, O12 = 019051. (20 
23 “32 31°13 12 “21 


5. Transversely isotropic solids 

A particular case of the class of orthorhombic solids worth special mention 
is the class of transversely isotropic solids, that is, solids which have an axis 
of elastic symmetry. The strain-energy function contains five independent 
constants and may be obtained from the strain-energy function (15) for 
an orthorhombic solid by putting 


, —T , , . n l/p . D5 
C11 = Cea, Coz = Cg, C55 = Cam» Cog = 211 — Cyn), (28 


when the x-axis is taken parallel to the axis of elastic symmetry (cf. (2) 
p. 160, groups C%, etc.). 

In order that the transformation to isotropic form may be carried out 
for a transversely isotropic solid, its elastic constants must satisfy two 
conditions, viz. —. ee ee 26) 

Cyg + 2Cgq = (C41 C33)° = 31 C11/Cr2, \ 
which may be obtained from equations (16) and (25). Other sets of necessary 
and sufficient conditions may be derived from those given in section 3 for 


orthorhombic solids; we mention here only one set, namely that the polar 
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iagram of the reciprocal of the fourth root of Young’s modulus is an 
lipsoid of revolution about the axis of elastic symmetry, and 
O31 13 = Ofe- (27) 
This follows from equation (20) and the fact that, owing to the transverse 
sotropy, 19 O91; F320 O31 and Oo3 O73: 
The changes of variables which enable one to effect the transformation 
isotropic form in the case of transverse isotropy may be readily obtained 
by putting LZ, FE, in those given in section 4 for the case of orthorhombic 
symmetry, and need not be given separately. 
We now make a brief digression to state a new result in the general theory 
f transversely isotropic elastic solids which has some connexion with the 
main topic of the present paper. For a transversely isotropic solid of general 
type (whose elastic constants do not necessarily satisfy our special con- 
litions (26)), it is a straightforward matter to verify that the differential 
equations of equilibrium are satisfied identically when the displacement 


mponents w, are expressed in the form 


c Cd. 
Ue (d, ; dbo) t $ 
CX, * CXs 
F 2d. 
Us (d, + dp) - 4 (28) 
i COX>5 " CX, 
4 
Us (hy Jy +key by) 
= 
terms of three functions ¢; which satisfy the equations 
; aati y 3)? QO (a 1, 2, 3, not summed) (29) 
vy Cxs Ons] 
It are otherwise arbitrary 
Here, v, and v, denote the roots in v of the equation 
C14 Ogg ¥° + (CT 3 —C yy Cgg + 2C 13 Cag) + C39 C44 U, (30) 
2C44 (C11 —Co) (31) 
Bal I. . » » 9 29 
nd Ky (C11 %qy—Cqg)/(Cig + Cqa) (oe I, 2). (32) 


This is an extension of a result due to Elliott (3), who gave the particular 

ise 6, = 0 appropriate to problems in which the component of rotation 

1 / CX, »/€x,) about the direction of elastic symmetry is everywhere 
zero. The present result will be applicable to a wider class of problems, but 
whether it is completely general or not is not known. 

It is instructive to relate this result to the transformation method of the 
present paper which applies to transversely isotropic solids whose elastic 
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constants satisfy equations (26). These equations are equivalent to the 
equations ; 
, V1 Vo = Vg; (33 
as can readily be verified using (30), (31). When these equations are satis. 
fied, it is obvious that the differential equations (29) can be put into 
isotropic form by the transformation 
tL, = 2, hy = Be, Sy == ES, (34 
We now return to the main topic of the paper. 


6. Transformation of boundary conditions 

We have seen that for suitable types of anisotropic elastic material the 
differential equations of equilibrium can be put into isotropic form, for 
which many solutions are known. In any specific problem it is necessary 
to find a solution which fits boundary conditions characteristic of the 
problem which may be expressed in or derived from one of two general 


forms, namely ' = 
‘ . u; = u on s(x) = const.t (35 


when the displacement components have prescribed values on the surface 
of the body s(x) = const.. (36 
or pin; = Y; on s(x) = const. (37 
when the components of surface force per unit area have prescribed values) 
Let n; be the components of the unit normal to the surface; these are given 


by the equations — , 
F OS OS |—-%3 ; 
- — a (38 
OX, OX, 


It is easy to verify that, if we apply to these equations the transformation 
(6) which we applied to the differential equations, equations of the same 
form in the primed variables are obtained, provided we define the primed 
function s’ by the identity 


nN; 


i 





a (se } = oie), (39) 
and the primed variables Y, by the equations 


Y, = 10,:¥;, 40) 
where n = (A,,n,n,)}. (41) 


[t follows that the equivalent isotropic problem is a problem of the same 
kind as the given anisotropic problem, the possible differences being in the 
shape of the bounding surface, the distribution of surface forces or displace- 
ments, and the distribution of body forces. General methods of solution 
available in isotropic elasticity theory may therefore be used; the isotropic 
solution may even be known already, in which case the only labour 


+ xis an abbreviation for x, 2, %3. 
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OF 





EQT 





nvolved in solving the given anisotropic problem arises from making the 
ppropriate transformation in the isotropic solution, An example of this 


will now be given. 


7, Example: contact of an elastic solid with a rigid plane 

\s an illustration of the use of the above method in the solution of actual 
roblems, we shall now use it to determine the shape and size of the contact 
egion formed when an elastic solid with a curved surface is pressed lightly 
gainst a rigid plane with a given normal force P (say). Hertz’s solution 
rthe isotropic case is well known, and we use it to derive the solution for 
, homogeneous anisotropic solid having orthorhombic symmetry of the 
special type considered above, i.e. the elastic constants are taken to satisfy 
equations (16) or (17) in addition to the conditions of orthorhombic sym- 
metry. For simplicity, we assume that the planes of elastic symmetry are 
parallel to the planes of principal curvature and the tangent plane at the 
point at which contact is to be made. 

\s in Hertz’s theory, we assume that the surface of the elastic solid is 
convex to the plane and sufficiently regular near the point at which contact 


sto be made to be represented there by an equation of the form 


x3 /ro+23/15 22); (42) 
where ? ure the principal radii of curvature; the coordinate system is 


axis is normal to the rigid plane at the point of contact 


chosen so that the a, 
ind the other two axes lie in the rigid plane and are parallel to the principal 
lirections of curvature. 

When the force is applied we shall assume that the contact region formed 


is bounded by an ellipse 


a2 /b?+-22/c? l. z, = 0, (43) 
whose semi-axes b, c are to be determined. As in the isotropic theory, this 
wssumption will prove to be justified because it enables us to find a solution 
satisfying all the conditions of the problem. 


As in Hertz’s theory, we assume that the contact region is so small com- 
pared with the rest of the elastic body that this may be treated as if it were 
infinite in extent in formulating the boundary conditions. These conditions 


nay then be written as follows 


ve lr a2 |) ' when 22/52 ts Cc? i. xy 0). 
O, wu, 1 (a constant), when x = o. 
Al] 7 ) whe n x2/b? x2 /c? > l, 2, 0. 
All except 7 0 when 23/b?+-a23/c? < 1, %, = 0. (44) 
dx, dx. P, the integration being taken over the 


9/9 
| ae/p2 « . 0 
Le | ¢ l = 1 . 
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These are the boundary conditions for the given anisotropic problem, 
To obtain from them the boundary conditions for the equivalent isotropi 
problem we replace the variables (23), wherever they occur in (44), by the 
variables (24). We then find that equations are obtained having the same 
form as the given boundary conditions (44) but with the constants P. Moy I 


, 


b, c, d replaced by P’, r,, ry. b’. c’, d’. where 
; a> 13 


E\3/ E? \} 
pr — P| -)*( a ae 
E,} \E, z) 


From this it follows that the equivalent isotropic problem is a contact 
problem of the same kind as the given one, but with different values 
(P’, rs, 73) for the applied force and the radii of curvature. The reason for 
this specially simple result is that the surface on which the boundary 
conditions are specified is a plane parallel to a plane of elastic symmetry. 

The solution of the equivalent isotropic problem may therefore be taken 
directly from Hertz’s solution. For convenience we use the form given by 
Timoshenko (4). Allowing for the fact that in our case one body is a rigid 
plane, and making the appropriate changes in notation (a, b, A, B, P, k,, ky 
in Timoshenko’s notation have to be replaced by c’, 6’, (2r,)-!, (2r,)-*, P, 
(1—o*)/7E, 0 in our notation), it follows that the contact region is an 
ellipse whose semi-axes b’, c’ are given by the equations 


bo oc’ (3P" 


(46 
n m | 2 





1 2\1 

| | 1 —o*)3 
E | 

where m, n are functions, tabulated by Timoshenko, of a variable 6 defined 


by the equation 


cos 6 = (r,—1r4)/(73 +75). (47) 


To obtain the solution of the anisotropic problem we use (45) to express 


(46) and (47) in terms of P, ry, 3, 6, c, and we use equation (20) to express ¢ 
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n terms of the Poisson’s ratios of the anisotropic material. After a slight 
rearrangement we finally obtain the result that the contact area is bounded 
yan ellipse as assumed, and that the semi-axes 5, c of this ellipse are given 


\ the equat 10ons 





(#, ! b (=)! c (3P BS ‘ ki * 1—o93 O59 \} (48) 
\E,} n E,) m \|2 (|r, Ht rg Ht| (#2, E,)*} ’ 


cos 0 (7 ; Es . r 4) ("3 E; + lo Es ty. (49) 
his completes the required solution, which I believe to be new. 
\s a check on these results, we may consider the particular case in which 
the elastic solid is transversely isotropic in planes parallel to the plane of 
ntact (x, = 0). In this case we have £, = E, and og, = og, so that our 


esults (48). (49) give 


b ¢ | £,)3/3P[1 L | ~*1—o33)4 (50) 
n ” | E, | 2 if. rg E, } 
cos 0 (7%3—1)/(%g-1e)- (51) 


Qn the other hand, the solution of the same problem for a transversely 
sotropic solid of general type may be obtained independently by a straight- 
rward calculation of Love’s type (2, p. 193), based on Michell’s solution 


5) for the corresponding “problem of the plane’; when our special restric- 


ons (similar to (26), but with the 2,-axis as axis of symmetry) on the 
lastic constants are introduced, the results obtained agree with (50). 
For a transversely isotropic solid in which the direction of elastic sym- 


etry is not perpendicular to the contact plane (in contrast to the last 
example) but parallel to it, the contact problem loses its symmetry and 
‘nown methods do not apply; the present methods were in fact developed 
nan attempt to deal with this case. Taking the X_-AXI8 to be parallel to 
E, and O93 O39 Ore (from 


e direction of elastic symmetry, we have £, 


27)),and the solution is found from that forthe orthorhombic case ((48), (49)) 
» be riven bv the equ wwiONS 
b ( (3P l v | l OF: )4 


n in | 
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where v = (E£,/E,)' may be taken as a measure of the anisotropy of the 
material. 


8. Discussion 

The limitations of the present method arise from the special restrictions 
on the elastic constants which are imposed by the mathematical simplifica- 
tion of the differential equations of the theory and may not be satisfied by 
any actual materials. If there are in fact any materials which satisfy the con. 
ditions, it is likely that they will be orthorhombic or transversely isotropic, 
because such materials automatically satisfy most of the required condi- 
tions. However, none of the varieties of wood or the metallic crystals in 
these classes whose elastic constants are given in Hearmon’s comprehensiy 
list (6) satisfy all the required conditions. 

The conditions consist of relations between physical constants, and it is 
natural to ask whether they have any special physical significance that 
indicates the kind of materials which might be expected to satisfy the 
conditions. For orthorhombic materials, three out of the five necessary 
conditions are equivalent to the condition that the polar diagram of the 
inverse fourth root of Young’s modulus shall be an ellipsoid (whose principal 
axes are normal to the planes of elastic symmetry); in an isotropic material 
this polar diagram is a sphere. One might expect such a condition to b 
satisfied in materials like synthetic textile fibres whose elastic anisotropy 
is due to their having been permanently deformed by a pure, homogeneous 
deformation from an initial, elastically isotropic, state in which their 
constituent molecules are randomly oriented and distributed. There are 
however, not enough experimental data available to test this rather tenta- 
tive suggestion. Even if this condition is satisfied, there remain two other 
conditions (between the Poisson’s ratios), and we have not found any 
reason for expecting them to be satisfied. It is a curious fact that on 
of the zine samples in Hearmon’s list satisfies the ellipsoidal polar diagran 
condition very accurately. 


For materials which do not satisfy the required conditions the present 


methods might give a rough picture of the stress distribution and displace- 
ment field in problems which have not been solved exactly by other means. 


Such a picture might be of use for highly anisotropic materials which maj 
come nearer to satisfying our conditions than the conditions of isotropy 
In particular problems a knowledge of this rough solution might give 

useful lead to the form of the exact solution. 


This work forms part of a programme of fundamental research undertaket 
by the British Rayon Research Association. I am indebted to Dr. H. ©. 


Howell for suggesting the problem of generalizing Hertz’s contact theor 
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where t 
SUMMARY 
General direct methods of solution are provided for both the stress and displacemer y 
. ° . — nere ut 
in a wedge under force or couple nuclei acting at any point. The method is to ad; } 
Tranter’s Mellin transform treatment of the stresses to the complex variable approa 0 
and to extend it to deal with the displacement. Solutions are then given for a wedg 
under an internal couple nucleus with stress-free or rigid boundaries. 3 For 
An i 
1. Introduction the orl: 


A TREATMENT of generalized plane stress expressing the mean stresses and 
displacements in terms of complex potentials has been given by Stevensor 
(1) and here we shall follow the notation of that paper. In particular} “eT 


complex potentials giving rise to force and couple nuclei are also givei 


by Stevenson (2), but it appears impossible to find complex potential _— 
in algebraic form which would remove the resulting flank stresses over 
the wedge boundaries. Tranter (3) has developed a Mellin transfor 
treatment to find the Airy stress function for wedge problems with For 


distributed flank loading. This, with some modifications to suit the com-| :, = 2 
plex potential method, has been extended to include a treatment for the} tise to 
displacements. It will be shown here that the usual displacement functio 





can be side-tracked and the displacements derived directly via the Melli 
inversion integral, so that they involve only the Airy stress function trans-| 
form. Solutions for normal and tangential point loading on the flanks wer 
given by Shepherd (4) by Airy stress-function methods, but no attentiol 
appears to have been given hitherto to the problem of force or couple nuclei) 4, TI 
in the material of the wedge. The object of this paper is therefore to put) — str 


forward a definite method of solution for these problems. The 


2. Notation 


We shall use the mean stress combinations in polar coordinates 7, 4 givel 
where 
ié 


by z= ax-+1y re 


integ: 
9] © 


0’ = +00, ®’ = *—60 +278, 


[Quart. Journ. Mech. and Applied Math., Vo’. VIII, Pt. 2 (1955)] 
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1 can be expressed in terms of complex potentials Q(z), w(z) in the 


20) (2’(z) + Q’(Z), (2.2) 
20’ — 20"(3) 4+ ~0"(z) (2.3) 


ibsence ot body force 
‘he complex mean displacement D’ is given by 
D U.+iU,) = e-”D, (2.4) 


t 


e the cartesian mean displacement D is given by 


SuD K42(z) 20)"(z) @ (2). (2.5) 
is the rigidity and + 3—4o, where a is related to Poisson’s ratio v 
1—o)(1+-v) - (2.6) 





Force and couple nuclei 
\n isolated force singularity F F.+-iF, acting at a point C, which is 


origin of a complex variable z,, is given by the complex potentials 


SSeS al bate.) P log z,, W(21) xk Pz, log z,, (3.1) 
! oP 
r¢ r ; (3.2) 
m(1-+-k) 
‘imilarly, an isolated couple singularity G may be represented by 
¢ 21G “— 
= 0) .(2; 0 w,(Z)) log z,. (3.3) 


r completeness in connexion with the problems detailed later, in which 


may be mentioned that these latter complex potentials give 


stress and displacement combinations 
= iG ] 
; (-)’ 0) (p’ ; (3.4) 
Me = (7 de®?) 
va 
dD : (3.5) 


tru r—de'? 
‘. The Mellin transform and its application to generalized plane 
Stress 
The Mellin transform f() f the function f(r) is defined by 
| rP 1A() ) dr, (4.1) 
0 
te p must be such that the integral converges. Provided that the 


eora ly is bounded for some | , an inversion integral 







































x 
since it follows from the definition of y that it is real if p is real. 
The stress combinations may be similarly treated since 


o* ofl 
: )x 21 . | Dy). 
ore, or\r 
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may be expressed in the form 
c+ia giving 
7 7 -—p { 
i (7) S(p)r-” dp, $9 
( “ia nd 
where c > k > 0. 
It is proposed to apply this transform to the equations which express th pom (4 
stress components and displacements in terms of polar coordinates, 
The mean stresses. The components 77, 60, 76 are given in terms of Airy'd Supposi 
T'} t rl t 60, 76 t fA Dy 
stress function y, assuming no body forces, by 
l dy l ex — o*y oy C l ex ie 
imi = —— » © / < 66 = = P ré | 4 . (4 Tran: 
r= Ob* ror or or\r oo ‘ 
jucing 
where Vey = @ | 
= oe ila.180 eres 
and V2 = —+——., 4.5] Equati 
Ore ror = r* 0? 
The mean displacements. Here U. and Up are given (see Coker and Fil 
(5)) in terms of y and a displacement function % satisfying 
: : On tra 
— C Os ‘ 
Vu QO: : (; V?x (4 
cr\ o@ 
as 
, CY Cw and 
2ul e x (l—o)r “ 
i ' 47} Elimir 
_ l ey , OW 
2ul f) - x i (] o)r- 
r 00 ob 
Integration of (4.1) by parts gives the Mellin transform of derivatives 
the form said 
na ( "1 (p+n) ; 
[ fore tn-t dr = fp) is 
J * ) 
0 (P) 
provided that 
. Inv 
pp+m-lf(m-l(r) +0 asr>Oandasr—>oo form ® ae nm. (404 
7 reac 
Writing D = d/d@ we may transform (4.4) and obtain 
(D?-+- p?){ D?+- (p+ 2)" 0, (4.10) Q’ 
the solution of which is 
Acir® + A4e-iv) 1 Beil +200 1 Be—iip+20 (4,11 
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g | & l dr (D2 p")x (4.13) 
| @’rPtl D-+-ip){ D+-i(p+-2)ty. (4.14) 
19) when » | it will be seen that we have here assumed that 
2)’ +0 asr—>0Oand asr—>oo. 
nosing © O ) at infinity, we have, since it is finite at the origin, 
2 p< n—Z. (4.15) 
[ransformation of the displacements may be accomplished by intro- 
9 liary displacement function %, given by 
“us yp. (4.16) 
ions t.6) ther take the torm 
—: ( ] Cus 9 - 
v(t) <0, &(1 4A) _ gay (47) 
‘ cor\r co 
ransforming these equations, we have 
/) (p- 2)* hb, ) (4.18) 
} | Dis, (D?4 p)x. (4.19) 
nating 1.7) now become 
(D224 97) y 
| / (1) DY (] a) } “ 
(p 1} ) 
(4.20) 
D( D+ py 
| 7 Dy—(1—o Px 
(p+ 1)(p+2) 
ilae. Application of (4.2) to the stresses and displacements 
L\ b l¢ ids Lo 
| 
= ax Dp \ Dp 
+ [ (D+ipD+i(p+2yiRr-»-* dp (4.21) 
l r a)(D—in\yi D ,( 1 2)! 
| ; D ip)| | Py AP ‘ly? I 1 dp 


(p +-1)(p |. 2) 
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The boundary conditions. These may be conditions of displacement 
of stress. Using the complex boundary stress r@-+i@ we may write 
(ro ibO)re +} dr (p+ 1)(D+ip)x 4) 
0 
so that if | (70+-i00)re 1 dr = t(p) 
0 
then (p+1)(D+ip)y = t,(p) on é x 
t,(p) oné Ms 4.2 
where 6 a are the flanks of a wedge with origin at the vertex. Similar 
if the displacements are given on the flanks and 
} (U,+-iUp)r? dr d(p), 
0 
then 
i , (1—o)(D—ip){ D—i(p+2)} 
(D-+ip)|1 a) p) (p+2)}]3 Kio it 
2h (p+1)(p+2) 
d,(p) oné \. 


We may of course have displacements given on one boundary and stresses 
on the other. 


5. Isolated couple nucleus at z — d (real) with stress-free boundaries 


Using the vertex of the wedge as origin, the complex potentials (3.3) be 


come 





Further 


415). 


eading 


and 


where 


The st 


in wh 
and v 
the li 
the in 
integ 


take 
W 


whe 





_ , aes 
iz) = © w(z) logz, where z, -d. 5) 
These give rise to flank stresses 
i" ons 7 | 
r6-+-200 —, ond X. ).2 
2a (r—de+t)2 
To remove this stress we accordingly take 
t.( 7 { rPHdr 
)) ———= 5 
i\F ? : (7 de i m2 
0 


which may be evaluated by a contour integration giving 


Gd ipa (p | 1) 


t(p) . 
NE 2 sin pr 


t,(p). 


also 


t,(p) 


| wel 


and 
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1 


ement rther, for convergence of the integral —2 < p < 0 which agrees with 


rite t,10 


ie boundary condition (4.23) accordingly gives 


iGd”. ip(m—a 





{e'P* 1 B(p+ 1ljee+20— Be-ty+2a . 
4sin pr 
: o iGdretPa-a 
At tt B( p+ 1)e-tp+20e_ Betly+2)a ; 
4 sin pr 
ne te 
iGd K(p. x«)sin prt G(p, «)cos px p 
{ : u u : I I (0.0) 
tpG(p Sin pir 
iGdv 
tu - 
B . ; (5.6) 
4G(p, x) 
G(p, a p+ 1)sin 2a — sin 2(p+1)a ) me 
° ( 4) 
K(p, x) p+ 1)cos 2a cos 2(p- L)a J 
stresses may now be evaluated from (4.21); we find 
G *  (d\” (p+1)sin(p+2)6 ~ 
(-) | i - I dp, (0.3) 
qT ; ? G(p, x) 
vhich the complex p in replaces the real p of the Mellin transform 
| } ) re} } 
| we must take 2 0. If we attempt to press the value of c to 
imit zero, the integrand of ©’ must be examined for possible poles on 
imaginary axis. The zeros of the denominator G(p,«) in the range of 
ndaries| Integration occur at p = 0 |, —2 but only p = Oisa pole. We therefore 
ke ¢ 0 and put p in (5.8). 
Writing now G(p, x) E,+10,, 
here 
E, n 2a(1—cosh 2na), 0, 7 SIN 2a cos 2asinh 2yx, (5.9) 
EB FE, (sin 26 cosh 76 7 COS 24 sinh nf) } 
O,(7sin 26 coshné + cos 26 sinh 7), 
U / sin 26 cosh 76 cos 26 sinh n@) 
O, (sin 20 cosh 74 — 7 cos 26 sinh 78), 
G | . *jeimn( EB | iO) 
ly y 
¢ ”) } 11m | — = dy 7R 
Tr | . | H+ 0} 
R is the residue of the integrand at p = 0, so that finally 
G * Ecosm Osin myn asin 20 - 
Q Jo vy) dy site wee, (5.10) 


E;+ OF sin 2~— 2a cos 2a 
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Similarly ®’ may be reduced to 
; G a d\P p+l)etr? — 
D’ — — . | _ _ Je — Z(p, 8) dp, (5.1) 
2nr* | r} G(p,«)sin pr 

where Z(p,@) G(p, «)cos pr+-[K(p, «)—(p 2)e2?|sin pr. (5.12 
Since Z(—2, 0) = 0, p = 0 is again the only pole of the integrand. Writing 

E, cos 2a(1-+-cosh 2x), 

O, = 7 cos 2a — sin 2a sinh 2yq, 

E, = —O,sinhzy-+ E, cosh 7y+ ne?sinh 7, 


O, = E,sinh y+ O, cosh 7n—2e?sinh zn, 
E, E, E,+ 0, O:+74(O, E3— E, Os), 
O, = E,0,—0, E3+-7(E, E;+ 0, 03), 
we have finally 
A) i* 508 KN y si iv 
@' = — Ja [ Pscoetnt Ree dn —oR,|, (68 
2 sinh 77( £}+- OF) 


: 27(cos 2a —e?"?) 
where k=m+i0 and 7k, = 1+———. 
sin 2~— 2a cos 2a 
The displacement is given by (4.21) as follows: 


(D—ip)(D+ip){ D—i(p+-2)x 8i(p+1)(p+2)Be-t+, 





° y — G , ‘d v Y @ ip0 
leading to U,+-iUy = | (“) (p, O)er™ 


r} sin px G(p, x) 
where 


Y(p, 0) sin px|{K(p, xjJ—(p-4 1 )e2? 1 (3—4oa)e—2 | 1 Gi(p, x) COs pr, 


Eval 
of « al 
involv 
the 01 
from (: 


G 


stress 


‘ 


Total 


6. S 

AY 
a ser 
may 
intes 


have 





(5.16 
so that, writing 
E, = E,coshzy + sinh wy] ne?” — O,4-(3—4e)ae sinh 278], 
O, = O, cosh ry + sinh z7| £,—e??+ (3—40)e-2"eosh 28], 
E, = E, E;+0, 9;, O, = O; £,—E;0,, 
we have 
U,+i0y = SH |2 ( Resin bat Oecetn gy eR}, (67 
Smpur| , sinh 7n( £7+ OF) " 
0 
al 2 cos 2a —e29+- (3—40)e-2? 
where 7k, 14 "| saeco a = 


sin 2a— 2a cos 2a 
The complete solution is now obtained formally by adding to the stresses 
and displacements of (5.10), (5.13), (5.17) those due to the original complex 


potentials as given in (3.4) and (3.5). 


see | 


EF; 


whe 
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Evaluation of these integrals does not appear possible for a general value 


fy and recourse must be made to numerical computation. The work 


nvolved is of a very elaborate character, but calculation of the shear stress 


the only non-zero component) on the axis of symmetry has been made 


from (5.13) for the case « 


re 
' 
n 
io 








lor; the values are given in the following table. 


aR 
y* 
2 


2Np $S2J%sS 


/ 3 

rs / 
0 i 4 
a i a Cc 
a — - ___ —— ooo | 
2 Wd 1 ¢8 
~ Pal _ Total stress I< 
, ° 
Oe ——-—-—Stress due to boundary load t2 6 
a 


Shear stress r@. 


6. Semi-infinite plate under couple nucleus 





An algebraic solution of the above problem can be found for the case of 


semi-infinite plate corresponding to the value « 


‘7 since the integrals 


iy then be integrated directly. This accounts for the fact that explicit 


ntegrated forms for the complex potentials cannot be found in general but 


have been found for this special case by Stevenson (2). It is of interest to 


see how the 


From (5.10) we have 


ibove results reduce to and agree with these. 


(~) isin 26 = I, sin 20-- J, cos 28, 


| dn, 


C cosh nA si n cosh 76 cos 
where I. | cosh 7@ sin myn -+- 7 cosh 76 ¢ OS M7 
J sinh 77 
0 
* sinh 70 COSMy 7 sinh 7 sin mn 
I, = | a eet he. 


sinh 77 
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The values of these integrals are given in (6, p. 277) and lead to 


4Gr sin 6(d-+-r cos @) 





0’ se (6.1 
a(r2+-d?-+ 2rd cos 6)2 
’ Gi de 210 r(e 310 1 De i) 
a - 6.2 
(d-+_-re i)3 (9 
: Gi [rile 210) Ke-200 
D (6.3 


diz (d—+-re i0)2 (r+de i?) ; 
Stevenson shows that the boundary stresses may be removed in this cage 
by the complex potentials 


Q2(z) ; w(2) Nog Zo+—}, Za = 2+d, (64 


the stresses and displacements from which agree with the above results, 


7. Isolated couple nucleus at z = d (real) with rigid boundaries 
Using the complex potentials (5.1) the displacement of the flanks is 
Ean ali iG l 
Uti, on @ +a. (7.1 


4 r —de*'* 
To remove this displacement we accordingly take 
. eo 
iG r? dr 
4a | r—de ta’ 


0 


=) 


14.2 


d,(p) 


which may be evaluated by a contour integration giving 
iGd”s ip(a—a) 
d,(p) . =~» 1.3 
4x. 81n pa 
also we must have d,(p) d,(p), (7.4 
and for convergence of the integral 
I p=. (7.5 
The boundary condition now becomes 


iG dPe-'v7 


pA + B( p+ 1)e?!*+- « Be-t2p+2)a : : 
4 sinpz 
; LG dPeip7 
and pA B(p | L)e Zia” Bei2p+2)0 d = 
4 sin pr 
leading to 
Gd 
A H,(p, «)cos pxr+-J,(p, «)sin pr}, (7.6 
4pH,(p, eo pet nas ee NPs Pa} 
B iGdP 77 


iH (p.a)’ 
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H,(p, x) Pp 1 )sin 2« t kK Sin 2(p + L)x ) 
; (7.3) 
J (p, ’ 1 1 )cos 2~—x cos 2(p- l)a } 
3 may now be obt ined from (4.21) as follows: 
G *  (d\” (p+-1)sin(p+2)6 * 
(-) Mt Bhs eso! ia dp, (7.9) 
. ? H,(p, x) 
(). 
ise the integrand has no poles in the range since x + 1 and so 
0 and put p . Writing 
H », MX) XxX, iY; 
in 2a K cosh 2na y; 7 Sin 2x | K COS 2a sinh 2a, (7.10) 
XY, (sin 26 cosh 76 7 COS 24 sinh n0)-+- 
Y,(7 sin 26 cosh 7@ +- cos 26 sinh nf), 
X sin 26 cosh 76 cos 24 sinh 6) — 


Y, (sin 2@ cosh 74 — 7 cos 26 sinh 78), 


2G [ E,cosm7 O. sin m7 Pa 
fe) k& al ‘dy. (7.11) 
. Xi+YF 
] pP x) X- i) ; 
2a(1 K COSN 27a); iF 7 COS 2a -« sin 2a sinh 2na, 
Xo Y, sinh T X , cosh nz ; ne?’ sinh 7, 
if X, sinh n7-+Y, cosh nw—2e?” sinh yz, 
ae ee Se De 
O, = X,¥,—-¥, X40 RAM 
Gi | * O,coskn+F, sin ky ~ 
@ ake dyn —1). (7.12) 


2a \~ | sinh n7( XF t Y?) 
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Writing 


X, = X, cosh nz + sinh yal ne??— Y,+ nie? sinh 276], 


Y, = Y, cosh nz + sinh ynx[ X,—e??+ xe? cosh 28], 


Ey = X,X,4Y,Y,; 0,=Y,X,—Y,X, 


is ae Gi * E,sin kyn+ 0, cos ky 
then U, Tt il A > 4 9° 7 ue os 7 ( u] ] ° (7.13 
Sar sinh ry(X}+ Y}) 


0 
Again, to complete the solution, the stresses and displacements due 
the original complex potentials must be added. 


8. Isolated force problems 

The above theory and problems forms part of a thesis entitled ‘General 
ized plane stress in an elastic wedge (under distributed and isolated loads) 
for the Ph.D. degree of the University of London (1953), in which isolated 
forces are also applied to points of the wedge. 

Forces acting on the flanks at z, = de*'* may be treated using complex 
potentials which include terms of type 


Q(z) = Elogz,, w(z) = —Ez, logz,—z, Elogz,, (8.1 
whilst a load at a point on the axis requires 
Q(z) = —Plogz,, w(z) = «Pz, logz,+dP logz,. (8.2 


A fundamental difference between the form of solution in these two cases 
arises from the intrusive « in the w(z) of (8.2). This has the effect of doubling 
the labour involved. Since the solutions to these problems are so elaborate 
they have not been detailed in this paper but may be referred to in the 
thesis mentioned above. 


9. Conclusion 

The facilities offered by the Mellin transform have allowed mathematical 
solutions to be obtained for nuclei of force or couple acting at any point of 
an infinite wedge, so that the usual tentative approach offered by an Airy 
stress-function method has been avoided. At the same time the use of the 
stress combinations has shortened the amount of work involved, and often, 
particularly in the general theory of section 3, has led to a neater expres- 
sion of the results. 
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PROBLEMS OF THIN CLAMPED ELASTIC 
PLATES 


By R. TIFFEN (Birbeck College, London) 


SOMI 





R ed 9 April 1954} 


SUMMARY 
with prol f transverse displacements of thin elastic plates 
pakiuaien | fl 1) half-planes, (6) those which can be mapped 
nfinite strips, the boundaries being clamped. 


Solutions t te problen e been given by several writers, Stevenson (1), 


~ 
Ww 


The present investigation aims at general methods 


lar solut Proofs of uniqueness are given, with emphasis on 
the cas f terials extending to infinity in some, but not all, 
t ited interior loading are considered for all the above 


1. Fundamental equations 

THRoUGHOUT this paper rectangular Cartesian axes O(x, y) are chosen in 
the mid-planes of the plates. The notation is that due to Stevenson (1). 
\ll relevant mean stresses and mean stress couples may be expressed in 


rms of w, the transverse displacement of the mid-plane, where 


Viw = f(z, y). (1.1) 


The ven lution to (1.1) may be written 
Q(z) +-2O(Z) + w(z)+@(2Z)+ F(z, 2), (1.2) 
ere | is a particular integral of (1.1) and Q(z), w(z) are functions 
which are analytic in the region R of the z-plane which is 
cupied elastic material. Problems of given surface loading and 
specified conditions of boundary clamping are considered. Thus, w, éw/en, 
where lenotes differentiation along the normal to the boundary C of 
the mate ive assumed given along C. Such information yields imme- 
boundary values of éw,cz. The existence of solutions of such 
problems y be shown by methods similar to those used by Sherman (4). 

2. Uniqueness of solutions 

It i imed that the region R is bounded by contours Cj, C,,..., C. 
nd a finite surrounding itour C,.,. Cuts are made to make RF simply 
nnected d the contours (C),..., C,,., and the cuts, traversed once in 


ach direction, are referred to as the contour C. Let W,, W, denote two 
systems of transverse displacements corresponding to the same surface 
loading of the plate, including isolated forces, which give singularities in 


(Quart. Journ. Mech and Applied Math., Vol. VIII, Pt. 2 (1955)] 
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w, and to the same boundary conditions. Then the function y = w 
satisfies the equation Viy—0 9) 

1X -. 
in FR, and, on the boundary of the material, we have 


fo 
CZ 
Further, since w,, w, have the same specified isolated singularities in B 
are uniform, and have uniform partial derivatives, except at these singu- 
larities, the function y is free from singularities, is uniform, and has wi. 
form derivatives throughout the open region defined by R. 

Taking 


X = 20 9(z)+20,(Z)+-w,(z) +a, (2), (2.3 
CX -f\! - , 
then ‘ 20 (z) + 225(2) + wo(z), (24 
Oz 
Cc 2y t= B- 
—— = 2,(2) +2, (2) (2.5 
Cze2 
O34, Ev. 
Ps - 0),(2), (2.6 
C202~ 


where accents denote derivatives. 
The boundary conditions (2.2), together with the uniformity of the 
above functions, make the integral 


- Tov e2v Re. 
| | | x < x dz x : Xx as| 2.1 
P CZ C202 Zz 
( 


(a. 
c 


equal to zero. Assuming continuity of Q5, Q5, QF, wo, wo in the closed 


region defined by the material, 7 may be transformed to a surface integral 


i 2,,\2 3 : 3 1 
° Cc 4 va € d CY C § Cc 4 ’ 
r= 2i f |(2X)', x Axo @x _, Ox.) as 
p C202 CZ 0202" 02 Oz0z" ~ dz*0z* 
2 
Thus | [Q6(z)+-2,(2) 2 dS = 0, (2.8 
R 
so that 5 (2) +04(2Z) = 0 (2.9 
throughout R. Thus (5) Q),(z) wcz-+ a, (2.10 
where c is a real constant and a is a complex constant. 
From (2.4), (2.10) it is evident that w)(z) is uniform. Also dy/éz = 0 
OM Gh. C41, 80 that 


w(2) 
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nthe boundary. If z is any point in the interior of R, 


l f —adZ . 
Wal ~ me - Xe (2.52) 
zm C¢—z 


use of the Cauchy integral being justified by continuity and uniformity 


hus vo(Z) 12 +B, (2.13) 
ere B is onstant, throughout R. However, on (C),..., Crs X 0. so 
t her edi c8-+5) +B-+B-+(—ae—az) = 0 
, 3 ik (2.14) 
here /: is eal constant Finally 
2) xyz+tk. (2.15) 
From (2.10), (2.15) it follows that 
\ 0 1.¢ wy uw, (2.16) 
ougcnout R 


3, Uniqueness of solution when (’,_, tends to infinity in all direc- 
tions 
\ssume initially that the conditions specified at infinity imply no 
restrictions on (2, w which satisfy the conditions 
Q o(z-*), Q o(z), Q,w° = O11), w= Oz). (3.1) 


n, at infinity, y satisfies the conditions 


2 3 
C7 o"-y c 
O(z) : O(1), X — o(:-), x 


CZ C202 CZ2C 


o(z-*). (3.2) 


-9 


Hence the integrands in (2.7) are o(z~1) at infinity, and, since C’,,, may 


be chosen to have a length which is O(z), the integral takes the value 


ro as ¢ tends to infinity. The equations of the preceding section hold 
s far as (2.12). which is replaced by 
, » (+adl ; P ht? xc = 
2ziw(z) = S | + lim = ot. (3.3) 
— 9 Soe J $4 
l ( Cn 1 
rhe first integrals are zero, whilst the last integrand is O(f-!) at most. 
\ssuming that (7) 
CW(C = ‘ 
in ——— “f (3.4) 
mmplex constant, the same in all directions, then 
w(z) x. (3.5) 


, 


shows that « vx. The remainder of 
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The demands made in (3.1) may be justified as follows. Let @ denot 
any large contour enclosing (C; C,, and a large circle @4 which itse 


genes n 


surrounds (},..., (',. Since @ lies entirely within the elastic material 
Cy¢g — = Cyg [O’(2) + 0'(Z)| = 0. (3 


Hence CygQ'(z) 277 


where y is a real constant. 
It can be seen that 
Cy¢Q"(z) = Cygw"(z) = 0. (3.8 
Thus outside @, Q’(z)—ylogz, w”(z) may be expanded in Laurent series 
and it can be shown that the most general forms of Q, w at infinity are 


Q(z) yzlogz+-A logz+B+ } (a,2"+b,2-"), (3.9 


n=1 . i 
w(z) = Azlogz+-8logz+Cz+D-4 > (c,2"+1+d,, 2-"), (3.11 


where y, 5 are real constants and A, B, C, D, a,, b,, ¢,, d,, are comple 


constants. It is evident that, from (3.1), the conditions specified at infinit; 


enable the determination of the constants y, A, a,,, ¢c,, but no others. Th 
terms specified at infinity give all finite stress couples and infinitesimak 
other than those which are 0(z-1). Since C’,,, has length O(z), the abov 
infinitesimal terms have the unique values required for equilibrium of the 
material, as in the case of generalized plane stress (5). The smaller terms 


automatically take unique values, so that the solution is unique. If th 
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material has boundary of length O(z) at infinity but does not tend to infinity 
in all directions, then it is sufficient to specify all the terms at infinity othe | 


than those given in (3.1). That these conditions are also necessary is see! 
from an example given in section 9. 


4. Plates in the form of half-planes 

The plates are chosen to occupy the region y > 0 and two methods 0! 
solution for w under specified conditions of clamping along y = 0 at 
given briefly. They lead directly to methods of solution of problems 6! 
plates which can be mapped conformally on to the half-plane and t 
solutions of problems involving infinite strips. 


where 


ind £ 
lirect 


youn 
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{part from a constant, the boundary conditions along y 0 are deter- 
ined by (€w/éz),-9. Taking w in the form 

w =()(z 20)(2)+-w(2z)+-@(2), (4.1) 
a —— — 30)'(z)+-.0(2)+w'(z). (4.2) 


w,. where 


w3(2 2004 (2) +0, (2), (4.3) 
¢ displ icements Ws; where 
f cow, 
treQ,(x) = (“") f(x), | } 0. (4.4) 
i 0 Cy y=0 
\oain, the notentials 42,, ; where 
(2, (2) —Q,(2), (4.5) 
e displace ments w,, wher 
{OW>s CWs > 
*) 0 #1m £2,(x) ) Q\x). (4.6) 
; CA y=0 \ oy y-0 
Hence the problem of the half-plane with 
Cw ° = 
2 =| 1q(x) h(x), say, (4.7) 
s been reduced to the determination of functions Q,, Q, which are 
nalytic in y > 0 and have specified real or imaginary parts when y = 0. 
If 
Q) 2 S2, w W1-+Ws, (4.8) 
w 20)" (2) +O, (2) —Q, (2). (4.9) 
if n be expressed as Fourier integrals, then (6) 
1O(z) = A(z fin’ (2) 2H’ (z)+H*(z), (4.10) 
wi - 7 [{ h(t)dt 
here H | | A(t)e-" dt (4.11) 
T ; 7 2—t 
| H*(z) is the transform of h(t). The Cauchy integral is the simpler for 
lirect evaluation, but the Fourier form is more suitable when a second 


undary exists in the finite part of the half-plane. 


5. Isolated interior loading of the half-plane 
The potentials 


>A 2, log (z) = } —z, A,(z—z,)log(z z,) | 
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where D is Love’s plate constant and W, are real constants, correspond t ' 
isolated loads W, at z,. Let z, be in the region y > 0. If the plate is clamped 6, So 
so that let e 
cw . 
| - Q, 5.9 ymntoul 
Ce | y=0 n to tl 
it is clear that additional potentials free from singularities in y >*0 are | to the 
required to counteract the effect of (5.1) along y = 0. The logarithmic | from a ‘ 
terms given by (5.1) in the boundary condition may be removed by th 
addition of the ‘image potentials’, 
; n n 
Q,(z) = > {—A,(z—z,)log(z—Z,)},  w,(z) = > {4+2, A,(z—z,)log(z—z,)}._ | Let the 
r=0 r=0 
(5.3 
-_ < zy) : 
Thus QO, = 2,42, = > A,(z—z,)log ; (5.4 
r—0 o— 2) 
If 
$ {= Ate—2)log?— mo 
Wo WaT Wp hes p<“ zr) Og — =. \v.0 hen, i 
r ol ~ Z,J 
If wy, is the displacement given by Qo, wo, 
u - (2 \(z—Z,) 
cw “ ~~ 0 — Lin NS Sy » = 
‘ S A ,(—i,)log ——~_— +. I(z, 2), (5.6) } If 
CZ hon 2—Z2Z,)(Z—2Z,) 
- mn - - - __ | then 
where I(z,Z) SY A (Z—Z,)(z,—Z,)/(z—Z,). (5.7 
r=0 : 
pee a Thus 
re CW x , 
Thus xt — I(x, x) J (2) 4 (z,—Z,) (5.8) | and w 
Cn /y=0 r=( 
and 
Cc Uo . : CW9 >: ; ind 
= 2re J (x) f(x), | 2im J (x) g(a). 
ox y=0 cy }y=0 Te fy 
(5.9 AT TS), 
- . solution 
The solution may be completed by the method of the previous sectio! 
Equations (4.4), (4.6) are satisfied by 
20, (z) J (2) 20,(z). (5.10) | Where 
Hence, using (4.3), (4.5), the potentials additional to Q,, Q,, w,, wy, are 7. Iso 
‘ | 
0,=0,42,=0, w,=a,+o, | J(z) dz. (5.11) Lett 
: Fin (5.] 
HA,= A, A, = A,= A, 1 0, the complete solution is give | displac 
by the potentials 
Q. = A(z—zp,)log|[(z—2)/(z—%p)], (5.12) 
w Az, (z- z»)log| (z 29)/ (%— 29) | (Zo 29) Az. (5.13) 7 Also 
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, Solution by conformal mapping on to a half-plane 
im pe ; : 
Let elastic material occupy a region R of the z-plane, bounded by a 
tour C, and let the function z(Z), where ¢ = €+7%7, map R conformally 
the region 0 of the ¢-plane, the boundary C corresponding 
60 a the axis 0. The boundary conditions along C are given, apart 
it] a constant, by (cw/cz From (4.2) 
Q'(£)/2’(f) +. Q(f)+-w' (Z)/2’ (2). (6.1) 
the req ed boundary conditions be 
| Cw \ ‘ ° a . 
a—) f(€)—ig(€) h(€). (6.2) 
0 ‘ Y |} 0 
€)Q,(f) —2(2)Q4 (0), (6.3) 
intl 
u CW, x 
i treQg(é), | } 0 (6.4) 
} A CY }n=0 
Wo 2'(C)Q,(C) —2(C)Q4(Z), (6.5) 
rT CW, ° ‘. 2 
: | 0 | *) $1m 0),(€). (6.6) 
Cr |} 0 cy 0 m 
hus the boundary conditions (6.2) are satisfied by choosing Q = Q,+Q, 
where 
freQ, é tim 22,(€) g(é), (6.7) 
(¢ €)(Q, —Q,) —2(f)(Q4 +Q4). (6.8) 


fy the conditions of Fourier’s integral theorem, then the 


1() H(C), ti’ (C) 2 (C)H*(f)—2(C)H'(C), (6.9) 


I’ H Dt uined by 1’é placing 2 by C in (4.11). 


Isolated interior loading 


Let the r isolated loading be given by the potentials 2, w, given 
.1), the points z, lying in the region » > 0. For the corresponding 
5 Sp emel 
n 
| a “ > = ] ¥ m4 s ar 
9 >A log| (z—z,)(7—Z,)|+ > A,(Z—i,). (7.1) 
0 r=0 


C.)b,(C), 
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where ¢,(¢) are analytic and non-zero in » > 0 since 2’(¢) does not yanis 


in 7 0. 
Thus ‘ 
Q.(6) = > A,(z—z,)log| (C—C,)6,(0)], (7: 
r=0 
nr oe ae Gs = 
w(f) = > {—2(¢,)A,(z—z,)log| (C—L,)4,(¢) ]}, (74 
r=0 
whilst (7.1) gives 
w = - - , , > > ¥\ T > ik - = 
pied > A,(Z—Z,)log| (¢—C,)(C—¢,)b,(0)6,(0) | + >} AMZ—2Z,). (75 
CZ r—0 r=0 


The most convenient image potentials in this case are 


n 


2,(0) = ¥ {—A,(2—2,)logl (C—£,)4,(2)]}, 1 
r=0 
n _ . _ . 

w,(C) z 4 +A .2(¢,)(2 z,)log] (¢ l)b,(C) I}, :: 
r=0 


which are free from singularities in » > 0 and give 


€ Wy, 


n = . 7 a . 
> A,(Z,—2)log| (¢—%,)(€ —f,)b,(0)4,(2) | 
Cz r=0 
S Alem) die Egg, as 
a (C,—C)4,(C) dz 
For the combined potentials 
n ; 
QO, = 2,40, > {+-A,(z—z,)log[ (¢—Z,)/(C—¢,) }} 7 
r=0 
n 2 a - . 
Wy WW, Wp >4 | =(C..)(2 z, log] (Z BS) (C Celits 7.1 
? 0 
uw 9 CMF - 
im 1 (2—2 log § See — Sr) | NGO (7.1 
Cz Ler (6—G,6—4G,) 
: (C—£,) $,(2) 
ae 3 is >\.T,F 4) Gere 3 - 
where Wé,2)= S Ae ¢.)o,(¢) ——_+- — (7.12 
ya (C & 2(C) 
r 0 
1 C Wy’ bo "1 
Thus | l(é,€) = J(é) (7.1 
CZ J y=0 
ne _ - ~~ 
where J(é) = 3 {4+ A,(6,—G,)b,(E)b,(E)/2'(E)} (7.1 
r=-0 


ig(E)| 2. (7.18 


The logarithmic terms required to give the isolated loading 
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near in the boundary function J(€). To complete the solution the 
thods of the previous section may be followed, analytic functions being 


sen to satisfy (6.7). Completion of the solution is exceptionally easy 


J(0), ie. [(f, 0), is free from singularities in y > 0; for, in such cases, 
] re the real and imaginary parts, when 7» = 0, of the same 
ytic function, viz 2J(C). Hence (6.7) are satisfied by 
202, (¢ J(C) 202,(C). (7.16) 
is the required potentials are 


() () u 


2'(C)J(C) dZ. (7.17) 


~ 


general, however, the functions ¢,(Z), and therefore J(¢), are not 
ible in 0 and f(é g(€) are not usually the real and imaginary 


rts of the same function of ¢, analytic in y > 0. 
8. Isolated load in plate with clamped parabolic notched boundary 
If rd 7a)" 202, (8.1) 


0 corresponds to the parabola as 2a*y x? 0 (7). Let there 


single isolated load V W, at z 2) corresponding to 
.+iB (B > 0) 
s Co)(C +, + 20a)/22 (8.2) 
Dy C+-¢,-+ 22a) /2e. (8.3) 
potent (),. wy in (7.9), (7.10) become, replacing A, by A, 
Y) { Cy) (¢ 2ia)log = i Wo $a (Co 1a)?Q.. 
> SO 
(8.4) 
\ls J { €,—21a)(C+-¢,+ 22a)/2(C+-2a), (8.5) 
nd is free from singularities in 7 0. Hence the additional potentials 
given by (7.17 In the present case 
3al?+-30(C,+ 2ia)(l,—2ia) ]/6. (8.6) 
9. Uniqueness of solution 
Consider the transformation 
id) Piz, (9.1) 
er Ss a positive integ The above method of solution for a plate 


unded 1 0 may be applied without difficulty. If, however, in (6.5), 


() 
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so that, apart from a possible constant, the clamping conditions aloy; 


» = O are not affected by the potential in (9.2) and w, given by (6.5), yj 


Wo if(C+ia)"/2—12(C-+-1a)"+1/2(mn+-1)—2(C—ta)"*"/2(n+1). (94 
Thus Wy if" (nm —1)/2(n+1)+O(C") at infinity. (93 
Hence, at infinity, 

Q2, O(C) O(zZr), Wo O(E"+) O(2e stdin), (Q 


Since n may be made indefinitely large, it is clear that the requirements 
Q) O(1), w O(z) 97 
made in (3.1) cannot, in general, be reduced. The solution given in 


previous section is not unique, for, at infinity, the potentials Q — Q,—( 


Ww = wy+w, may be written 
Q ABL+O0(1), (ABO /64-0(2), 3 
and when » 2, (9.2), (9.4) give 
ih, = é, w, = ((C—bal?+ 3ia2Z) 6, 9.9 


which are the same orders of magnitude as (9.8). If (9.9) is multiplied | 


AB and added to (9.8), the unique solution 


Q(C) -LiA(E—Ly)(F4+- Cp +-2ia)logl (f—L_)/(C—Z,)|4+ APL, (9.1 
w(C) 1A (Z, va)*(C— Co )(S+Lo 2ia)log| (¢ Co)/(¢ Co) 


C?AB(ia+ 2a)/2—i ABE (C,+ 2ia)(Cy—2ia)—a?]/2, (9.1 


is obtained, where Q = O(1), w = O(2) at infinity. 


10. Modified treatment of the half-plane 

Throughout this paper the principal aim has been to consider inter 
forces. This has not hitherto been done for infinite strips, and, in order t 
consider problems involving these, a modified treatment of the half-plane 
is necessary. The previous method, which was the most suitable for exter 


sion to conformal mappings, leads to difficulties with regard to convergent 


of integrals when applied to infinite strips. The transverse displacement 


wis given by x 





w 20)(2) +-2022(Z)+- w(z)+- @(2), (10. 
Cou l cu -Cu pone _ i 
so that 2 202’ (z) +-O(Z)+w’(z). (10.2 
C 21 O2 cy 
If w,(2) 20,(2)+-2 | Q,(2) dz, (10.5 


‘ of 
then Ws o=9 tre Q), (x) dx. (‘ ) 0. (10.4 
™ y-0 
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[f w,(2) 202,(z), 
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then (10.6) shows that w, satisfies the conditions given for ws, in (11.3). 
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OW. . , 
then ees (), | "| 41m Q2,(2). (10.6) 
: cy 0 7 
Hence the boundary conditions 
cw on 
we A(x). | = Giz) (10.7) 
CY | y-0 
\ be obt ined by choosing 
tre | Q,(a) da h(x), fim Q,(x) = g(x). (10.8) 
With the notation and restrictions of (4.11), the transverse displacements 
piat e given Dy 
K (z)—iy| K'(z)+-0G"(z) |}. (10.9) 
is assumed, henceforth, that specified conditions of clamping along 
0 have been obtained by the above method for the boundary y = 0 
finite strip 0 < y 1, the appropriate potentials being denoted 
0, Q \dditional potentials are required to give specified 
y 1, without altering the conditions along y 0: ef. 
t] treatment of generalized plane stress (8). 
> e e e cu 
11. Potentials which give w, , | 0 
cy 0 
If Q,+2 | Q, dz, (11.1) 


Qod v1 iyo) Ut(y, iys)e sul du, (11.2) 
where » re real fun ms of the real variable u, then, from (10.4), 


‘ (3 0. (11.3) 


2Q,, (11.4) 


Q, ; 105) é (6, —20,)¢ ul du, (11.5) 


The ctions y,, yo, 5;, 5, may now be determined to give specified 
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conditions along y = a. Using these two pairs of potentials in (10,] 
(10.2) and combining the results. the following conditions are obtain 


along y = a. 


s . ‘ 
Wy a | feos xul —8, as+y.(s—Ac)|-+sin rul —8, as+-y,(s—Ac)}\ du, 
0 ai 
fOw oe , , . . : 
| | {cos ru| —sudy,—96,(s Ac) | f sin xu V1 SuA—0,(s+-Ac)]|\ d 
OY | y=a 7. 
0 / | }- 
where A au, ¢ cosha, s sinh A. 
Let the conditions required along y = a be 
ae ; 
Wen | (a, COS TU o.s8in vu) du. Ls 
0 
cu \ af : 
} | (7, cos. ru+7,sin ru) du, (11.9 
cy rr = 


ya . 


where it should be noted that the functions ¢,, 04, 7,, T. allow for the effect 


of the potentials Q,, Q,, 


w,, w, as well as the clamping conditions along 


y =a. From (11.6)-(11.9), we find 


C7 05a8 Y(s Ac) 
Oo 0, a8 (8s Ac), 
(11.10 
Ty Yo SUA— 0,(8-4 Ac). 
T> y, suA—4,(s+Ac), 
whence 
(s Ac)o> AST. (Ac 8)T> suros 
2 2 5 ie 2 
me 2)? st) 
(11.11 
(s Ac)o, AST, (Ac S)T, sudo, 
Yo — ; O» es . 
si #—® 2» 


The derivation of (11.6), (11.7) required two differentiations with respect 
to z under the integral sign of the integral for Q, and one for Q,. Validity 


of the operation requires that 


0}; Fo o (u-4), Tas To o(u-*) at wu 0. (11.12 
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7, are bounded at u 0, it is evident that 
2A7/5)—ar,(1-+-A*/30)}+-O(1), 
-A“/ 3D) az,(1 A2/30! O(1), 
($2.3) 
O(1) 
O(1). 
1.5) at u 0 are 
u~) iy.(20zu) O(1). 
7 (12.2) 
*u*) +-16,(27zu)+ O(1), 
\difications are required if the solution is to 
potentials Q., w, and Q,, we, where 
(12.3) 
{ “u*/2)+-(y,—tyo) l—izu—z?2u?/2)| du, 
(12.4) 
10 du. 
where. Hence (11.2), (11.5) may be modified 


1] —ozu-+z*u*/2) 
iyo )(e l--izut-22u?/2)|du, (12.5) 
1) + (6, —16,)(¢ 1)| du. (12.6) 
6) at u 0 are bounded. It can easily be 





il signs are also bounded at u 


ned by differentiations of (12.5), (12.6) with 


0. 
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13. Isolated load in the interior of a clamped infinite strip 


From (5.12), (5.13), the potentials 


Q),(z) A(z—z,)log(z—z,)/(z—2Z,), w(Z) 2 §29(z) + (Z)—2) Az, 

(13 

correspond to an isolated load at z 2%) in a semi-infinite plate y > 0 wit 
Ww,» (Cw/0z),, 9 0. In the following work z, r+718, 0O<8 


From (10.1), (10.2), the displacements along y = a, given by (13.1), are 


2) 


, Bs x—a)*-+-(a—B)? 
Wo y=0 Al (wx v)*+(a B)2}log | Biot A 1aBpA (2). 
(x v)“+(a- B)? 
(13.2 
Also 
Vo x )21-(a—B)? SAaB(a-+-B 
2 "| 2Ala B)log \ a ( <a — ) x(2 
cy}, (v—a)?+(a-+B)? © (a—a)?+(a+p)? 
(13.; 


If (x), x(a) are expressed in the forms of (11.8), (11.9), it is found that 


o, we" Am cos au a, wel" Amsin au 
Bu cosh Bu(aw-+-1)—sinh Bu(l+-au-+-aBu2), (134 

T, ue Az COS xu T, wee" Azsin au 
aBu cosh Bu-+-sinh Bu(—B--a-+-aBu). (13.5 
Since o,, o, 7, 7, are O(1) at u 0 and are O(e%-®") at u 0, tl 


solution is valid. 


14. Conclusion 

Whilst these problems show much in common with the corresponding 
problem of generalized plane stress, their interest lies in the fact thai 
each presents its own particular difficulty. This does not seem to hav 
been appreciated by writers such as Muskhelishvili (9, 10). 
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VOTE ON THE ACCELERATION OF LIN’S PROCESS 
OF ITERATED PENULTIMATE REMAINDER 
By A. ¢ AITKEN (Edinburgh) 
R ed 14 October 1954] 


SUMMARY 


terated p timate remainder, introduced by S. N. Lin for 

ymial, is examined from the point of view of 

g or ace ting convergence. It is shown how this can be done 
proposed nomial by a suitably determined catalytic poly 

rl rocedure is nu I v illustrated for a linear and a quadratic factor. 


1. Introduction 
[HE present note arose from perusal of the note by Morris and Head (1) 
n Lin's method of iterated remainder (2) for approximating to a polynomial 


tor of a given polynomial. We had elsewhere (3, 4) already answered 


f 
7 


raised by Morris and Head on criteria of convergence, and 
ud even ide the suggestions we now develop. These suggestions seem 
ikely to be valuable in computation, and therefore deserve more than the 
immary treatment we gave them. They had reference to a device, which 

shall call the method of the catalytic multiplier, for removing divergence 


| greatly accelerating the convergence of Lin’s iteration. 
2. General criteria of convergence 


It is necessary to summarize the main result of our former investigations. 


t the polynomial to be factorized be 


? 1 1,4 1 Ay x! : ( ya 
Let L,,(2 b, z™-1+-b,a"-*—...+-(— )™,, 
be a polynomial factor of , and let an approximation to d,,(x) be 
b, (by +-€9)a™-?—...+(—)™(On+€m)- 
Let the quotient when /{ is divided by d(x) be 
; . , wie: Co Xt” ieee. \ one m* 
Let € denote the vector fe,, €5....€,,. of small errors in the coefficients of t,,,(2). 


Lin’s p cedure consists in dividing f,,(x) by t,,(a) as far as the penultimate 


m 


stage, at which the remainder, which we proposed to name the penultimate 


narnd r.), is of the degree m of the divisor. Dividing the p.r. by a 
nstant » that the leading term becomes 2”, we obtain the reduced 


penuttima mainder (r.p.r.). The r.p.r. is now taken as a fresh divisor, 


[Quart. Journ. Mech. and Applied Math., Vol. VIII, Pt. 2 (1955)] 
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and in this way penultimate remaindering is iterated, line by line, on th; 
work-sheet. The question is whether the resulting sequence of successiy; 
divisors, really a sequence of vectors having their coefficients for elements 
will converge to d,,(x) and, for practical purposes, converge well. 

We treated (3, 4) this problem from the point of view of the line 
transformation of the error-vector e. the elements e; being assumed s 
small that their powers and products of degree higher than the first wer: 


negligible. With d,(2) we associated a canonical matrix 








m . ol 


familiar in matrix theory as the companion matrix of d,, (a), its characteristi 
* I m 

polynomial, though here transposed from the more usual form. With thi 

r.p.r. we associated the reduced penultimate quotient (r.p.q.), namely 


n—Mp—1 S»n—m » »n—m-lj_p w—m—2 | n—m—lp “l 
) ‘ a Cyt Cot Sts ( ) Cn—m-1%}- 


The fundamental theorem was that the matrix which transforms th 
error-vector € at one stage of r.p.r. iteration into the error-vector at thi 


next Is 


R = (—)"-™e7 1, { Be-™—c, Bum (—)nmte BY. 


1—mi non 
The condition of convergence is that if p be the latent root of R of largest 
modulus, then p 1; the smaller |p| is, the more rapid the convergence 
By standard matrix theory the latent roots of R are the values 
assumed by the r.p.q. when x takes the values €,, &5...., &,, these &; being 
the roots of d,,(~) = 0. We thus have a criterion of convergence based 
jointly on d,,(2) and on the r.p.q. Naturally such a criterion in practic 
can be used only with tentative approximation, since both divisor and 
r.p.q. emerge only de proche en proche; but experience has shown that it 
serves well, and even (3, 4) provides further useful information concerning 
] 


a 


»(t). The results of Morris and Head are equivalent to cases m lm=2 
of the general theorem. 


3. Simple linear illustration 
Morris and Head, illustrating the linear case m = 1 by the simple exampl 
f(x) (w—1)(a—2)(~a—3) xv—6x?+-lla—6, 


point out that r.p.r. iteration converges in the case of the factors 2—1 ant 
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3, but not in the case of x—2. This is because the respective values of 


1 p.q ure aS follows 


r—l1;r.p r?—5x)/6; for x l, p s 
xr—2 “_47)/3; I 2,2) +. 
r—3 ry x) /2 r 3, p 0 
» case of 3 discloses a feature not explicitly commented upon by 
Morris and Head, namely that since x—3 is also a factor of the r.p.q. we 
> 0; the effect being that the r.p.r. process obliterates the linear 
t of error, leaving as residuum that higher type of convergence called 
ulratic convergence, characteristic, for example, of the Newton—Raphson 
teration 
In this instance the feature is accidental. We can however—and this 


s in fact (4, 333) our suggestion—deliberately produce it as follows: 


Let a be the trial divisor and f, (2) be the original dividend. Now 
eplace this dividend by (a—h)f, (x), where A has been so determined that 
her.p.q if h )y (7) with respect to the divisor x x iS zero when x a 
\s we have said, this can be done only approximately; but it is open to us, 

d often advantageous, as improved values of the desired factor emerge, 


recalculate / and thus further accelerate the convergence. If the factor 


t interests us is not linear but quadratic, d,(x), we shall replace f,,(2) by 

h,x+h,)f,,(x), where h, and h, have been so determined that the r.p.q., 

the p.q. simply, of (a h,x+h,)f, (a) shall contain d,(x) as a factor. 

Similarly the case of a divisor of degree m would require the introduction 

ta catalytic multiplier of degree m: but for practical reasons m is unlikely 
exceed 2 


4. Theory of the catalytic multiplier 


We consider the general case. Let q,,_,,(x), the quotient arising when 
is divided by d,,(2), be again divided by d,,(x), the division being 
ntinued into negative powers as far as may be required by what follows. 
Let the re he 
/ . ky yn—2m—-1 i fh, an—2m—2 
We require that the p.gq 
| fA, x™-*—...+-(— A, 3, (2) 
with respect to d, (x) shall be exactly divisible by d,,(v). This p.q., being 
penultimate, is a poly ial having zero for constant term. It follows 


expansion of the formal product 


1 > 


yi? l > » . ‘ » an 
? a2 h (a? 2m -ky x’ 2m—1 kh, x’? 
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the constant term and all later terms shall vanish. Taking the first ; 
conditions for such a property, we have 


hin k, 2m+j Tt h,, kn -2m+j+1 Tee hk, mt+j-1 k, m+j 0 
(j = 0, 1, 2,..., m—1) 
m simultaneous equations for h,, hs,..., h,, with a persymmetric matrix 


This is the simplest practical way of expressing the conditions determining 
h,, hg,..., h,,; for the quotient 


Sf, (%) (d,(a)}* Tn m(t) (d,,(2)} 


will be run off on the machine (always approximately, as we insist all 
through) by the usual synthetic division, and the coefticients 1. Ky, Bye 
will be shown in a single row. 


In the linear case m | there is a single equation, yielding 
h a 
In the quadratic case m = 2 there are two equations 
hak, sth, k, _3+k, . = 0 
hak, 3 h, k, ot k,, 1 0 


and it is hardly necessary to point out how, for the solutions, three isobari 
determinants of second order will at once be computed from k,,_ 4, /,,3.k 
k,,_,. four consecutive coefficients visible in a row of working. 

These two cases are likely to be the most valuable in practice, and will 
now be illustrated. 


5. Numerical illustrations, linear and quadratic 
The Laguerre polynomial of fourth degree is 
L(x) x*— 162°+ 727?— 9627-24. 
One of its zeros is near to x = 4-5. We shall find by r.p.r. iteration, with 


suitable catalytic multiplier, a linear factor near to x—4-5. The work-sheet 


with the convention of sign we personally favour (indicating elementary 


symmetric functions of roots rather than alternately-signed coefficients), 
appears as follows: 


$°5 I Ties 20°25 13875 


11-25 §5°5125 
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First, / 15°5125/11-25 = 5 nearly. We therefore multiply L(x) by 


and proceed with r.p.r. 


p-! 
7 106°125 4375 120 
77 100 IOICc $514 120 
I 105 575 O*45140 120 
Vith rapid convergence we obtain the factor 2«—4-53662. Two things 


y be noticed; first, that r.p.r. iteration on L,(a) unmodified would have 
ngly diverged in this case; second, that though the catalytic multiplier 


5 is quite close to the factor «—4-53662, no ‘fouling’ has been caused. 


However, this in its turn provokes the cautionary remark, that the device 


fthe catalytic multiplier will not succeed if the first approximation to a 


tor is too crude. There is a certain radius of effectiveness. 

\s a second example, given that two roots of L(x) 0 are not far from 
£5, 2 |-75, let us combine these in the approximate quadratic 

tor x7—6°3 7-9. Unmodified r.p.r. iteration will again be found to 


erge We proceed to determine R,, he: 


e second row shows | ky, ky, ko, kg. The equations for h,, h, are 
h,+-3-4h,—26-32 = 0, 
3 4/ 26-32h, + 138-956 = 0, 
elding / 6:03, he 5°82 
raking 62-+-5-8 as approximate catalytic multiplier, we multiply 
by it and proceed with rpr.: 
700°8 I 
7 F fg I10*°4093 
5 74°47 17 10°4182 I 
onvergence is rapid. The factor in question is actually 
628242 77-9198 (a—4-5366)(x— 1-7458). 
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CORRIGENDUM 
to the paper on ‘A source in a rotating fluid’ by 8. N. BARUA 


Che published figure is inconsistent with the table in the paper. It should 
be contracted in the axial direction by a factor 2. When corrected for this 


discrepancy, the bulge will appear to be less spread out along the axial 


direction. 
Ref. Vol. VIII, Part 1, p. 22. 


CORRIGENDUM 


to ‘A note on a paper by Davies and Walters on “The effect of 
finite width of area on the rate of evaporation into a turbulent 
atmosphere’ > by L. E. PAYNE 


The 7 in equation (12) of this paper (Vol. VII, Part 3, pp. 283-6) should 


be changed to a - 
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